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ABSTRACT 


This  report  presents  the  results  of  an  analytical  and  experi- 
mental research  program  to  (1)  improve  the  understanding  of  the 
physical  p-echanlsns  that  control  the  occurrence  of  pressure 
fluctuations  in  lonr  and  shallow  rectangular  cavities,  exposed 
to  hirh-spced  external  flow,  and  (2)  devise  and  evaluate  devices 
that  would  el  t tier  substantially  reduce  the  amplitude  of  such 
pressure  fluctuations,  or  totally  suppress  the  occurrence  o'' 


pressure  fluctuations, 


ft 


Durinp.  this  invest  lration,  the  understanding  of  the  complex 
interaction  of  the  external  shear  layer  and  the  cavity  Internal 
fluid  medium,  which  constitutes  the  generating  mechanism  of 
high-intensity  periodic  pressure  fluctuations,  was  substantially 
furthered.  hath  the  analysis  and  Investigation  of  the  physical 
mec nanisms  were  aided  through  extensive  shallow-water  flow 
simulation  experiments. 


Large-scale  experiments,  using  the  NASA  Lewis  Research 
Center  8 x 6 ft  Supersonic  Wind  Tunnel,  •.•■ere  conducted  to  sub- 
stantiate some  of  the  analytical  predictions  and  to  study  in 
detail  the  aeroacoustic  behavior  of  cavities  in  the  length-to- 
depth  ratio  ranee  of  2.3  to  5-5,  and  the  Mach  number  ran.p;e 
of  0.8  tc  2.0.  Detailed  information  was  obtained  on  the 
normalized  levels  of  the  first  three  resonant  modes  in  the  cavity 
for  a range  of  cavity  length-to-depth  ratios  and  freestream  Mach 
numbers . 


r 


Several  concepts  for  pressure  oscillation  suppression  were 
developed  and  evaluated  in  wind  tunnel  experiments.  The  most 
promising  concept  utilizes  a slanted  trailing  edge,  which 
stabilizes  the  "free  shear  flow  above  the  cavity,  thus  effectively 
suppressing  discrete-tone  generation. 
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SECTION  1 
INTRODUCTION 


f Cavities  and  cut-outs  in  structural  surfaces  of  aircraft 

: that  are  exposed  to  high-speed  external  flow  can  Rive  rise  to 

intense  pressure  fluctuations.  The  amplitude  of  these  fluctua- 
tions can  be  of  such  magnitude  that  they  affect  the  structural 
integrity  of  nearby  aircraft  components,  of  sensitive  instrumenta- 
tion, or,  at  the  very  least,  interfere  with  crew  performance, 
communication,  and  comfort.  In  fact,  levels  close  to  180  dB 
at  resonant  frequencies  can  easily  occur  at  transonic  flight 
speeds. 

Substantial  effort  has  already  gone  into  the  study  of 
cavity  pressure  oscillations,*  resulting,  in  some  understanding 
of  the  physical  mechanisms  and  the  complex  interaction  of  the 
cavity  external  shear  layer  and  the  cavity  internal  fluid 
medium,  which  is  responsible  for  the  generation  of  the  highly 
periodic  pressure  fluctuations. 

In  1970,  the  aeroacoustic  characteristics  or  narrow  and 
shallow  rectangular  cavities  in  the  Mach  number  range  from 
0.8  to  3.0  were  investigated  by  Heller,  Holmes,  and  Covert 
(1970).  The  main  result  of  this  study  was  an  improvement  of 
an  analytical  expression  (first  developed  by  J.K.  Rossiter 
in  1966)  that  relates  nondimens ional  resonance  frequencies 
and  freestream  Mach  number.  However,  only  an  upcer  bound  for 
resonance  amplitudes  could  be  derived.  At  that  time,  the. 
phenomena  were  still  understood  too  poorly  to  predict  reliably 
the  occurrence  of  discrete  pressure  oscillations  and  their 
energy  distribution  within  the  cavity. 


In  1972,  the  U.S.  Air  Force  conducted  extensive  flight 
tests  using  an  R.F-AC  aircraft  with  a modified  SUU  Al  test 
pod.  Results  became  available  on  the  aeroacoustic  environ- 
ment of  shallow  cavities  in  the  Mach  number  range  from  0.6  to 
1.3  for  realistic  flight  environments  (Smith  e.t  al . , 197A). 

It  seemed  that  t’ne  previously  assumed  upper-bound  level  for 
resonant  modes  was  generally  too  high.  Furthermore,  it  was 
realized  that  each  resonant  mode  must  be  considered  individually 


* An  Air  Force  Technical  Report  (Heller  et  al . , 1970)  describes 
Drevious  relevant  researc.i  and  provides  some  50  references. 


for  prediction  purposes , since  each  node  shews  a different 
level  dependence  or.  trie  freestream  7iaeh  number. 

Tnis  lack  of  thorough  physical  understanding  has  also 
hampered  the  development  and  implementation  of  devices  that 
would  either  substantially  reduce  the  oscillation  amplitudes, 
or,  ideally,  suppress  the  occurrence  of  resonant  oscillations 
altogether . 

In  view  of  these  facts,  an  analytical  and  experimental 
research  program  was  initiated  to  (1)  further  improve  the 
understanding  of  the  mechanisms  that  control  pressure  oscilla- 
tions in  f lew-exposed  cavities,  (2)  investigate  onset  and 
feedback  mechanisms,  and  (3)  devise  techninues  to  control  the 
onset/feedback  mechanisms  in  order  to  eliminate  the  occurrence 
of  oscillations  or  to  reduce  their  intensity. 

This  report  discusses  the  research  program  and  its  results. 


SECTION  2 


PAST  RESEARCH  EFFORTS  AND  PRORLEM  DEFINITION 

2.1  Overview 

On*.1  of  the  first  cxper inental  studies  investigating  flow- 
induced  oscillations  in  a simulated  weapon  bay  environment  was 
performed  by  Karanacheti  (1955).  Discrete-frequency  acoustic 
radiation  from  the  cavity  was  observed  at  both  subsonic  and 
supersonic  external  flow  speeds.  Karamacheti  noted  that  acoustic 
intensities  were  higher  when  the--  boundary  layer  unstream  of  the 
cavity  was  laminar. 

Later  studies  of  the  unsteady  environment  in  and  around 
rectangular  cavities  were  conducted  by  Gibson  (1958),  Ingard 
and  Dear.  (1958),  Leupold  arid  Baker  (1959),  Plumblee  et  at. 

(1962),  Quinn  (1963),  Last  (1966),  Hossiter  (1966),  f nee  (1966), 
White  and  McGregor  (1970),  Covert  (1970),  Heller  et  al . (1970), 
Bilanin  and  Covert  (1973),  ana  Smith  et  al.  (197*0.  T’he  Mach 
number  rang,e  covered  in  these  studies  is  0.1  to  5-0,  and  the 
cavities  had  length-to-aepth  ratios  (L/D)  between  0.12  and  10. 

In  general,  these  studies  showed  that  all  but  the  shallowest 
cavities  resonated  at  discrete  frequencies . 

2.2  Excitation  Mechanisms 

Cavities  are  often  classified  as  being  shallow  (L/D  > 1) 
or  deep  (L/D  < 1).  Covert  (1970)  explained  this  distinction 
in  terns  of  the  effect  of  the  image  vorticity  reouired  to 
satisfy  the  cavity’s  rigid  wall  boundary  conditions.  The  image 
vorticity  on  the  floor  of  a shallow  cavity  tends  to  excite 
longitudinal  modes,  whereas  the  image  vorticity  in  tne  fore  and 
aft  bulkheads  of  deep  cavities  tends  to  excite  depth  modes. 

Using  an  acoustic  monopole  distribution  to  model  the  pressure 
environment  in  an  oscillating  cavity,  Bilanin  (1973)  also 
found  that  the  cavity  tends  to  respond  in  the  direction  of 
greatest  dimension. 

Several  physically  motivated  models  have  been  studied  to 
Identify  the  mechanism  of  oscillation  for  shallow  cavities. 
Plumblee  et  al.  (1962)  hypothesised  that  the  environment  inside 
the  cavity  must  follow  the  characteristic  acoustic  response  of 
the  cavitv  (a  result  which,  however,  appears  to  hold  only  for 
deep  cavities).  The  mathenati cal  problem  consisted  of  solving 
for  the  radiation  impedance  of  the  cavity  opening.  However, 
the  shear  layer  spanning  the  cavity  is  capable  of  sustaining 
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wave  Miction,  and  it  is  likely  to  be  unstable.  Howe  (1970) 
thooret i cal ly  examined  the  transmission  of  an  acoustic  pulse 
through  a plane  vortex  sheer  and  showed  that  a grossly  in- 
accurate description  of  t.he  radiation  field  could  result  unless 
the  analysis  considers  the  vortex  shoot  dynamics. 

experimental  results  of  Kararrachoti  0955)  and  Heller  et  a l. 
(1970)  have  cast  doubt  on  the  assumption  that  the  forcin' 
mechanisms  which  drive  t.he  cavity  are  provided  by  fluctuations 
in  the  turbulent  boundary  layer.  Cavity  oscillations  are  most 
intense  when  the  boundary  layer  upstream  is  laminar.  "’hickeninm. 
the  boundary  layer  has  the  effect  of  reriuci nr  oscillation  In- 
tensity (Heller,  197*;)-  Rossiter  (1955)  also  concurs,  nolntinr 
out  that  when  a small  spoiler  is  placed  upstream  of  the  cavity  it 
reduces  oscillation  intensities. 

Pinca  excitation  frequencies  o'-'  shallow  cavities  do  not 
agree  with  the  acoustic  modal  frequencies  of  the  enclosure, 
investigators  have  sought  an  excitation  mechanism  that  is  not 
st.ron.~ly  dependent  or,  the  detailed  acoustics  in  the  shallow 
cavity.  Rossiter  (1966)  suggested  such  a mechanism  based  on 
feedback  similar  to  that  used  by  Powell  (idol)  to  describe  the 
production  of  edgetenes.  'using  a shadowrranh , Rossiter  observed 
that  periodic  density  fluctuations  travel  downstream  over  the 
cavity  mouth,  '"base-  fluctuations  'were  assumed  to  be  vortices 
shed  from  the  leadinr  edge  of  the  cavity.  The  shadowgraphs 
also  show  acoustic  waves  in  the  cavity  whose  primary  source  is 
near  the  trailing  edge.  Assuming  that  the  vertices  were  shed 
when  an  acoustic  disturbance  reached  the  leading  edge  of  the 
cavity,  and  that  acoustic  disturbances  were  generated  when  the 
vertices  reached  the  cavi-y's  trailing  edre,  he  was  able  to 
determine  a frequency  relation  at  which  oscillation  might  occur. 
Although  this  frequency  relation  requires  two  empirically 
determined  constants,  with  a suitable  choice  of  constants  a 
reasonable  agreement  with  experimental  data  can  be  obtained  over 
a moderate  Mach  number  range . Heller  et  al . (1970)  improved 
this  result  by  correctin',  the  sound  speed  in  the  cav..  ,v. 

2.3  Oscillation  Mechanism  for  Shallow  Cavities 

Pressure  oscillation  ir.  shallow  cavities  was  believed  to 
result  from  the  unstable  shear  layer,  which  spans  the  cavity, 
bein'  forcibly  displaced  by  disturbances  that  are  generated  at- 
or  near  the  cavity's  trailinr.  edre.  These  disturbances  were 
thought  to  arise  from  the  interaction  of  the  oscillating  shear 
layer  with  the  trailing  edge;  this  constitutes  a source  mechanism. 
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Alt  hough  these  dl  sturburjoes  should  force  1.  he  shear  layer  alone 
the  entire  cavity  length,  models  that,  assured  coup]  1 in"  to  occur 
only  at  the  leading  edr.e  of  the  cavity  have  riven  rood  predic- 
tions for  allowable  frequencies  of  oscillation.  Thus,  the 
shear  layer  Is  most  sensitive  to  pert urhat 3 ons  at.  the  leadlnr 
edr.e,  where  it  Is  thinnest.  An  explanation  Is  that  an  acoustic 
disturbance  reflectinr  from  the  forward  bulkhead  causes  a 
pressure  doublinp,  while  disturbances  in  the  external  flow,  if 
present.,  are  not  reflected,  and  a pressure  .lump  across  the  shear 
layer  must  be  balanced  by  deflecting  toe  shear  layer. 

Recently,  hilanin  and  Covert  (1973)  analysed  a model  of  the 
excitation  mechanism  suggested  by  Rossiter.  In  their  model,  the 
complicated  flow  process  at  the  trail!  nr,  edpe  of  the  cavity  was 
modeled  with  an  acoustic  monopole.  The  pressure  field  of  the 
monopole  is  used  to  drive  the  shear  layer*  spanning  the  cavity, 
which  yields  an  eigenvalue  relation  for  excitation  frequency. 

For  shallow  cavities,  the  acoustic  field  In  the  cavity  can  be 
adequately  represented  bv  a single  monopole  at  the  cavity’s 
trailing  edge.  The  effect  of  the  cavity  floor  was  to  influence 
the  Instability  of  the  shear  layer  by  changing  both  the  phase 
velocity  (this,  In  turn,  changes  the  excitation  frequency)  and 
the  amplification  rate  of  the  shear  layer  displacement.  Coupling, 
of  the  acoustic  cavity  with  the  shear*  layer  without  considering 
forcing  from  the  external  flow  limits  trie  validity  of  the  model 
to  y > 1.  Analytic  results  v/ere  In  excellent  agreement  with 
experimental  data. 

2.4  Onset  Mechanism 

Although  the  above  feedback  mechanism,  as  suggested  by  Rossiter 
and  modeled  by  Rilanin,  estimated  possible  excitation  frequencies, 
it  aid  not  predict  whether,  in  fact,  any  of  these  will  occur. 
Selection  of  the  excited  frequency  was  thought  to  be  based  on  a 
gain  criteria;  i.e.,  the  frequency  or  frequencies  at  which  a 
cavity  responds  must  correspond  to  the  mode  or  modes  receiving 
sufficient  gain  along  the  feedback  loop. 

To  illustrate  this  idea,  let  us  suppose  that  for  a given 
geometry  and  external  velocity  the  possible  frequencies  of 
oscillation  are  estimated  using  the  above  feedback  model.  ^he 
shear  layer  is  then  forced  at  the  leading  edge  with  a periodic 
pressure  fluctuation  with  frequency  corresponding  to  a possible 
excitation  frequency.  The  shear  layer  is  displaced,  and  vorticlty 
waves  propagate  aft  and  interact  with  the  trailing  edge  of  the 
cavity,  which  cause  a pressure  disturbance  to  be  generated.  Rpor. 
reaching  the  leading  edge  of  the  cavity,  this  pressure  disturbance 
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v.  '11  navi-  'li''  phase  a:-,  the  fore  1 nr  pressure  (the  frequency 

I'lii'.  i't:  1 :■  a possible  t al  ion  I'n’riu-ncv)  , but  I hr  amplitude 

i ■!'  ill'-  ii I url-unco  will  not,  in  renera!  , eauu]  the  fnre  Iti’ 
a:;  V.-.udv.  IT  I’m  ret  until  nr  pressure  disturbance  is  of  smaller 
amplitude  ill.!;:  that  o!’  tin'  fore  5 nr  pressure  (■’•ait:  < 1.0),  j t 
car.  I ■.?  at  it  icipaled  that  this  node  cannot,  r.u.'.taln  ore  1 1 1 at  5 on  . 

1 the  jvtumii.r  di  :*.t  urbat.ee  It  of  rreat.or  amplitude 
then;  ’.lie  herein-'  (i-.ain  > 1.0),  t!;r  osc  1 1 J at  1 or:  annli  t.ude  will  he 
3 ncreuseh  until  non]  3 no  if  orf'’cm  hav"  reduced  the  rain  tc 
nn!:y.  Oh1  ar.t'l  i t ude  at  which  a cavity  will  oscillate  3 s-  that 
ar.;  11‘ude  at  which  the  rai ::  j un'l.v:  1 . , wht-n  the?  merry 

addition  to  the  feed!  ad:  loop  I'ron  the  externa]  flew  hist 
I aiar.ee:'  the  ener.-.y  losses  1 hrourh  viscosity  and  radial,  i oh . If 
other  pea. title  frequencies  or  excitation  can  meet  this  pain 
ci’it>?:>ia,  i tier;  the  e xo  1 tat  1 on  of  si  null  ancons  discrete  frequencies, 
as  observed  in  practice,  1st  possible.* 

Covert  (1970)  examined  the  oscillation  phenomenon  as  a 
problem  In  hydrodynamic  st  ability.  lie  supped  cd  that  while 
the  external  flow  is  capable  of  supply  3 np  merry  to  the  en- 
closure, the  cavity  is  capable  of  dissipatinr  onerrv  t hrourh 
viscosity  and  radiation.  If  the  rate  of  onerrv  supply  is  less 
than  the  rate  at  which  onerrv  can  be  dissipated  fror  tdv>  cavity, 
oscillation  cannot  occur-  (rain  < 1).  As  the  external  flow 
velocity  is  Increased,  a point  is  reached  where  enerry  supply 
equals  cavity  enerry  dissipation  ( •••a  in  - 1).  Thi  s velocity  is 
termed  the  critical  or  onset  velocity.  further  increases  in 
external  flow  velocity  result  in  sustained  cavity  oscillation. 

A stability  analysis  of  the  vortex  sheet  across  the  cavity 
open  i nr  supports  this  description  of  onset..  However , '"overt  ’ s 
model  did  net  predict  excitation  frequencies  for  shallow  cavities 
since  it.  did  not  adequately  rodel  the  interaction  of  the  shear 
layer  with  the  trailing  cdpe. 


*In  Secs.  3 and  if , a physical  model  of  cavity  pressure  oscilla- 
tion will  bo  discussed  that  deviates  in  several  important 
aspects  from  the  "Rossi  ter  Model." 


SECTION  3 

PHYSICAL  MECHANISMS 


3.1  Introduction 


vhis  charter,  a dotal  led , but  csstnt  1 


• ?.]  '.’  qualitative, 
■:  rhenor-^na  1 s 


d.-sori  ption  of  steady  and  unsteady  cavity  "low  phenomena  Is 
presented.  t'lils  description  provides  physical  insirht  that 
"will  be  helpful  In  underst and!  nr  and  j nt  erpret  in.’  the  analytical 
offer;  s and  experimental  results,  f-uch  an.  under  standi  nr.  is  also 
required  to  appreciate  the  operation  of  oscillation  suppression 
devices  . 


3.2  Steady-Flow  Considerations 

Although  v:e  are  primarily  concerned  with  unsteady  effects, 
it  is  worthwhile  to  consider  briefly  the  somewhat  hypothetical 
case  of  steady  flew  over  a cavity.  This  flow  is  steady  only 
in  an  overall  sense,  since  shear  layer  turbulence  is  included 
in  the  discussion.  Pi  cure  1 is  a sketch  cf  a steady  cavity 
flow.  If  tite  flow  were  supersonic,  then  a steady  wave  structure 
at  the  leading  ar.d  trail! nr.  edf.es  would  also  be  present. 

The  flow  separates  at  the  sharp  leading  edge  of  the  cavity, 
arid  a shear  layer  prows  downstream . At  the  rear  bulkhead,  the 
flew  stagnates  and  splits,  fart  of  the  shear  Inver  enters  the 
cavity,  and  the  rest  passes  over  the  cavi'v  trailing  odro  and 
continues  downstream  as  part  of  the  downstream  boundary  lavor. 

In  the  laminar  flow  case,  the  dashed  line  across  the  cavity 
mouth  (see  Fig.  1)  denotes  the  streamline  that  divides  the 
flow  which  enters  the  cavity  from  that  which  does  not.  The  same 
fluid  always  remains  in  the  cavity.  The  fluid  that  is  entrained 
from  the  cavity  by  the  shear  layer  over  the  cavity  mouth  is 
returned  to  the  cavity  by  the  trailing-edre  stagnation  process. 
T!ie  same  is  essentially  true  when  the  shear  layer  is  turbulent, 
although  there  will  actually  be  some  change  in  which  fluid 
elements  occupy  the  cavity.  This  exchange  of  fluid  occurs 
because  of  the  turbulent  diffusion  across  the  dividing  stream- 
line. On  the  average,  however , the  process  remains  as 
illustrated  in  the  figure. 

It  can  be  concluded  that  this  must  be  the  flow  configura- 
tion for  the  steady  case.  The  shear  layer  always  entrains 
fluid;  therefore,  there  must  always  be  a stagnation  on  the 
rear  bulkhead  to  return  this  fluid  to  the  cavity.  A similar 
argument  establishes  that  the  entire  shear  layer  does  not  enter 
the  cavity.  Figure  2 illustrates  this  point.  It  is  important 
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ret  urn  i e t he  cavity  at  the  t ratline  edrc  . The  murnitudo  of 
tin:  entrainment  tv  t hn  shear  layer  fixes  the  level  of  t.  M s 
rec  1 rcu  1 at  i on . fills  out  ralnr^ut  also  controls  the  location  of 
the  s : riitnal  1 on  s '.roar  1 i tie  within  1 ho  shear  layer. 

Tin.'  presence  of  this  ( sometimes  substantial)  recirculation 
within  t he  cavity  affects,  the  unsteady  phenomena.  Flow  re- 
circulation acts  to  chance  the  notion  of  waves  within  the 
cavity  and  alters  the  character  and  the  dynamics  of  the  shear 
layer.  This  recirculation  flow  is  essentially  a vorcox  flov; 
contained  within  the  cavity  and  bounded  by  its  side  walls. 

The  sidewall  boundary  layers  provide  considerable  damn Inn  to 
the  recirculation  flow;  they  cause  secondary  flows  to  he  set  un 
within  the  recirculation  vortex,  which  leads  to  fluid  exchange 
between  regions  with  hir.h  and  low  anrular  :ncrer.tun.  This  is  a 
fa:;:i  liar  process  in  viscous  rofatinr  flov;  boundary  layers.  It 
is  possible,  therefore,  that  the  presence  of  side  walls  contri- 
butes as  t.iuch  resistance  to  flow  recirculation  as  the  cavity 
floor  and  front  and  rear  bulkheads.  finer  the  level  of  rlow 
recirculation  affects  the  shear  stress  at  the  cavitv  mouth, 

It  can  be  concluded  that  cavity  width  will  affect,  the  cavity 
drap . Since  the  presence  of  flow  recirculation  affects  unsteady 
cavity  behavior,  this  behavior  must  be  dependent  to  some  extent 
on  cavity  width. 

Of  course,  the  cavity  width  affects  the  shear  l.aver 
characteristics  directly  in  addition  to  influoncinr,  flow 
recirculation.  Other  factors  that  are  important  to  the  shear 
layer  properties  are  the  upstream  boundary  layer  thickness, 
the  flow  Mach  number,  and  the  Reynolds  number. 

Two  requirements  must  be  met  io  achieve  a steady  cavity 
flow:  a steady-flow  solution  must  be  possible;  and  it  must 

be  stable  to  disturbances.  To  find  a steady  solution  means 
that,  the  mass  flow  balance  requirement  must  be  satisfied,  and 
that,  simultaneously,  the  shear  layer  must  assume  a shape  for 
which  it  is  in  force  equilibrium  everywhere:  i.e.,  the  internal 
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.3.1  Wave  Motion  in  a Cavity 


The  wa‘.  or  tai  ls.-  si : .ula l i on  study  reveals  that  unsteady 
c-tier.  of  the  shear  layer  loads  tc  periodic  trass  nod’ t Jon  and 
e-ova-  a*  th.o  cavity  trailin-*  ociro . Tr.  a nh.nl  low  cavity 

3 practice,  a loneth-to-coplt  ratio  os'  \ wo  or  'renter), 
his  mass  addition  and  rer  oval  producer  an  effect  i ha*  3 c-  similar 
o replacing  the  cavity  rear  bulkhead  with  an.  c?s:i  1 1 a*  3 i.r  piston, 
t fact,  this  olt  . i 1 nrl  t y has  h-'>on  d-r  tons1  rut  vd  o oi:i  mo » . t al  1 v 
'in,'  t he  water  table.  Therefore,  the  effect  cf  r .nss  addition 
id  removal  at  t he  rear  i ulk.head  In  sor-e*  ir.es  referred  to  as  t he 
'sender  i sV-on"  effect.  Thl  s effect  ronera’ os  forward  i ravel  in,- 
avec  1st  t he  cavity  that  reflect  frets  the  front  !",i3i:!)-:nc  a nr. 
eoore  rearward  travel  In';  w a v < ' s . The  rorc.lt  inr  wave  structure 
r.  \ ho  cav3.iv  fore  c-rs  tire  shear  layer  In  ar,  unsteady  runner. 
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FIG.  3.  TYPICAL  PRESSURE  OSCILLATION  CYCLE.  (3)  Mass  Removal 
Stops;  (4)  Mass  Addition  Starts;  (12)  Mass  Addition 
Ends;  (13)  Mass  Removal  Starts. 
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however , that  the  process  ir.  essentially 
low,  particularly  as  regards  the  Internal 
subsonic  case,  the  forward  travelinr  wave 
c relative  to  the  external  flow.  The  other 


(see  Sec.  6). 


Some  of  the  irregularities  ci 


where  they  cross . An  Jr.no rtant  feature  of  the  diap-.ram  is  t'ne 
staggered  arranr.er.ent  of  the  nattern,  which  is  caused  by  the 
difference  in  the  phase  speeds  cf  the  upstream  and.  downstream 
waves.  The  upstream  wave  (moving  toward  the  leading  edge) 
travels  faster  than  the  downstream  wave  (moving  toward  the 
trailing  edge ) . This  difference  between  wave  speeds  results 
from  the  complicated  interaction  of  the  shear  layer  motion  and 
the  internal  wave  system.  Essentially,  the  shear  laver  appears 
"stiffer"  to  the  upstream  wave  (supersonic  relative  to  the 
external  flow)  than  to  the  downstream  wave  (subsonic  relative 
tc  the  external  flow).  This  interaction  between  the  cavity 
wave  structure  end  the  shear  layer  will  be  discussed  further 
in  subsequent  sections. 


Because  the  upstream  wave  is  supersonic  relative  to  the 
external  flow,  it  radiates  an  external  wave,  "hus , energy  is 
radiated  away,  and  the  upstream  wave  can  be  expected  to  decay 
as  it  travels  forward.  .Similarly,  since  the  downstream,  wave 
produces  a shear  layer  disturbance  that  is  subsonic  relative 
tc  the  external  flow,  it  can  be  expected  to  draw  energy  from 
the  mean  flow  and  amplify  as  it  travels  rearward. 

Therefore,  the  ’wave  pattern  in  a cavity  is  a medal  structure 
formed  essentially  from  the  two  waves  traveling  in  opposite 
directions  with  slightly  different  wavelengths  and  snatial 
amplification  rates.  The  spatial  envelope  of  both  waves  will 
expand  towards  the  cavity  trailing  edge.  The  results  of  the 
water  table  simulation  correspond,  es3ential.lv  to  a first  mode 
o('  osc  i 1 iauior. , i.e.,  only  one  wavefront  o f'  each  tvoe  in  the 
cavity  and  me  pressure  r;ode . Real  cavities  in  air  crten 
exhibit  hir.'.er  modes;  ir.  "act , the  second  rode  usuallv  dominates. 
These  modes  are  seldom  seer,  natural  lv  ir  the  ’water  table  experi- 
ment, probably  because  of  hirri  damning. 
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FIG.  4.  WAVE  DIAGRAM  BASED  ON  WATER  TABLE  EXPERIMENT. 


3.3.2  External  Radiation  Pattern 


From  water  table  simulation  results  and  the  Sc'nlieren 
photographs  available  in  the  literature,  it  is  possible  to 
construct,  a picture  of  the  external  radiation  net  tern  of  an 
oscillating  cavity.  This  nicture  is  consistent  with  the 
internal  wave  structure  described  above. 

Consider  a supersonic  flow  over  a cavity.  Figure  5 shows 
the  expected  wave  patterns  at  the  leading  and  trailing  edges 
of  the  cavity  for  different  shear  layer  positions.  When  the 
shear  layer  is  deflected  downward  at  the  trailing  edge,  thrc 
basic  waves  can  be  expected  (Fig.  5a).  As  described  previously, 
one  of  these  waves  is  the  unstrean  traveling  wave  and  its 
corresponding  shear  layer  disturbance,  which  trails  a traveling 
wave  into  the  external  flow.  The  second  of  these  waves  is  a 
quasi-steady  "bow-wave,"  which  is  caused  by  the  supersonic  ex- 
ternal flow  interacting  with  the  relatively  blunt  trailing-edge 
region.  The  third  wave  is  an  oblique  shock  that  occurs  when 
the  flow  readjusts  just  behind  the  cavity  trailing  edge.  V'hen 
the  shear  layer  is  above  the  trailing  edge,  as  when  a downstream 
wave  is  arriving,  there  seems  to  be  no  significant  external 
wave  system  (see  Fig.  5b).  It  must  be  recalled  that  the  shear 
layer  disturbance  associated  with  this  wave  moves  at  a subsonic 
speed  relative  to  the  external  flow. 

At  the  leading  edge,  the  picture  is  simpler  (Fig..  5o  and  d). 
As  an  upstream  wave  approaches  the  leading  edge,  the  shear  layer 
is  bent  downward  and  there  is  an  expansion  wave  at  the  edge  as 
well  as  the  upstream  traveling  compression  wave.  After  the 
upstream  wave  is  reflected  to  become  a downstream  traveling  wave, 
the  external  portion  trails  away  since  the  downstream  wave  is 
subsonic  relative  to  the  external  flow.  At  this  time,  the  shear 
layer  is  deflected  upward  and  a leading-edge  oblioue  shock  occurs. 

This  sequence  of  occurrences  can  be  compared  with  observed 
radiation  patterns  (see  McGregor,  1969,  and  Rossiter,  1966)  to 
arrive  at  Fig.  6a,  for  supersonic  flow.  Figure  6b  is  a 
similarly  derived  sketch  at  high-speed  subsonic  and  transonic 
external  flow.  In  this  case,  the  radiation  caused  bv  the  leading- 
and  trailing-edge  oblique  waves  is  absent.  Although  the  pr'marv 
contributor  to  the  subsonic  flow  radiation  pattern  is  the  external 
portion  of  the  upstream  traveling  wave,  another  wave  is 
also  seen.  This  wave  is  caused  by  direct  external  radiation  from 
the  unsteady  shear  layer  impingement  at  the  trailing  edge.  Be- 
cause the  flow  is  subsonic,  this  wave  can  also  move  upstream. 


15 


-»..U 


TRAILI NG-EDGE 


FIG.  5.  STAGES  OF  TRAI L I NG- EDGE  FLOW  IMPINGEMENT  (a  and  b), 
AND  LEADING-EDGE  FLOW  SEPARATION  (c  and  d)  WITHIN 
OSCILLATION  CYCLE  . 
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3.4  Physical  Model 


i !: 


1 1 i s l po.' :;  i !■  i •■?  V;;  .i .; r.i.-./i i ■ T " < * Mi-.  ■ 
o'cMiatioii  { ,-r*oo »'.*::•  and  discuss  core 
phene:-.' -not; . As  indicated  in  the  nn1 
addition  and  removal  process  at  t ho 


yd-::  t ion--:]  a s pec  t « of 
1 o;;s  ?pc<  i on:;  the 
c n v ' tv  t ra*  ] i n<-  od  re- 


in caused  bv  unstcad 
: rod  ..cos  a pint  on-li 
up  t!i-'  internal  wave 


t ot  j on  o:'  the  shear  layer.  '"hi  a process 
o effect,  nt  the  roar  bulkhead , which  sets 
structure  that,  forces  the  shear  layer. 


This  shear  layer  notion  is  responsible  for*  the  tral  1 ' nr-edre 
mass  addition,  ami  removal  . Tho  wave  motion  of  the  shear  layer 
and  the  wave  strscvurr-  wit.hin  the  cavitv  are  strenrlv  counled 
and  cannot  he  considered  separately  as  lorn'  as  the  wave  ler.rth 
is  comparable  to  or  exceeds  the  cavity  rier.th.  The  cavity  in- 
ternal wave  structure  is  composed  primarily  cf  upstream  and 
downstream  traveling  wave trains . Their  combination  produces 
a modal  unsteady  pressure  field  in  the  cavity.  nf  course, 
other  waves  are  present,  to  satisfy  the  detailed  boundary  condi- 
tions imposed  by  the  details  of  cavity  reor.etrv , shear  laver 
shape,  and  the  mass  addition  nrocess.  The  upstream  waves  radiate 
enerry  away  from  the  cavity  and,  therefore,  decay  as  thev  move 
upstream.  These  waves  reflect  from  the  front  bulkhead  to  become 
the  downstream  travelinr  wave  system.  The  downstream  waves  draw 
‘?!if-rrv  from  i hn  near,  flow  and  amplify  as  they  move  downstream . 
This  amplification  must  provide  erourh  shear  layer  displacement 
at  the  trailinr  edre  tc  sustain,  the  process  , overconinm  the  many 
sources  of  losses  in  the  system.  At  least  for  the  simplified 
analytical  model  cf  the  shear  layer  dynamics,  presented  in 
to.  the  spatial  amplification  of  the  downstream  wave  is 
considerable  ana  easily  exceeds  the  decay  rate  of  the  upstream 
wave  in  almost  all  cases.  This  fact.  orobabiv  explains  the 
si  arm  tendency  of  cavities  to  oscillate  in  one  or  more  modes 
ov  r a broad  and  continuous  ranre  of  cavity  reomotrlos  and 


f 1 ow  speeds. 


The  modal  frequencies  at  which  a riven  cavity 
for  a riven  set  of  flow  conditions  are  determined 
phase  and  amplitude  conditions  of  the  shear  laver 
trailinr  edre  as  they  relate  tc  the  mass  addition 
the  internal  wave  motion  of  the  rear  bulkhead.  '"r. 
edre  mass  add!  tier:  and  removal  pr-  ;ess  is  essentia 
dynamic  pnenomonon ; as  such,  it  is  quite  complex, 
where  vortex  roll  up  of  the  shear  layer  has.  o-’curr 
complexity  would  be  further  increased. 
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water  table  results, 
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d 1 st  r 


lenrlh.  The  apparent  traveller  wave 
layer  is  a result  of  t he  dotal  ran4 


a:;  p;  i f'i  cat  i on  of  the  downs;.  roan  traveiir.r.  wave  in  the  cavity. 

From  the  analytical  results  which  are  documented  later,  it  Is 
apparent,  that  the  spatial  a no!  i *'i  cat,  ior.  of  the  down  stream  wave 
is  preatest  relative  to  that  of  the  upstream  wave  when  the 
external  flow  is  subsonic.  1 n this  case,  it  is  reasonable  that 
the  shear  layer,  subjected  to  forcinr  by  the  cavity  Internal 
wave  structure,  would  roll  up  into  downstream  t ravel i nr  vortices. 
Thus,  the  appearance  of  discrete  vortices  on  the  cavity  shear 
layer  is  completely  consistent  w'th  me  oscillation  mechanism 
previously  described.  In  fac;  , at  subsonic  flow  speeds,  It  is 
likely  that  the  nonlinear  process  of  vortex  roll  up  provides 
an  important  am.plitudo-llmilinr  r.r  chan  i nr:.  Vie  can  conclude  that 
"vortex  sneddinr"  is  a manifestation  of*  the  oscillation  process, 
but  it  is  not  essential  te  the  under lyi nr.  mechanism. 


Finally , fhe  trailinr-edre  mass  addition  and  removal  pro- 
cess deserves  some  detailed  examination.  •"•ur  discussion  of 
steady- flew  considerations  revealed  that,  a mass  flow  balance 
must  occur,  whereby  fluid  entrained  by  the  cavity  shear  layer  must 
be  returned  to  the  cavity  thro  urn  a flow  sr.amnat  ion  process 
at  the  trail! nr  edye.  Clearly,  when  fhe  shear  laver  is  deflected 
downward  relative  to  Its  steady-state  position,  additional  mass 
enters  the  cavity  at  the  trailinm  edro.  '"'his  mass  add  it  lor.  is 
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fluid  r.orn.ailv  in  t.’nat  melon.  /•  s rioution.ed,  Tills  di  snlncerent 
has  an  effect  similar  to  reniacinm  the  t.  mi  iter-  ed>~e  bulkhead 


rear  ot 

ft  « 


wi  t 


a piston. 


w h 1 c hi  is  then  d i s i 


laced  by  a correspond i nr  amount, 
in  oscillatir.r  cavities 
show  that  the  trai  1 i nr;— core  rericn  does  not  behave  as  n hard 
v/all  ; 1 .e.  , pressure  inodes  do  no  .-enerallv  oxhlbii  a maxi  mum 
at  the  rear  bulkhead.  Tnc idental ly , the  mass  entrainment  process 


bv  ‘ he  shear  layer  is  not  necessarily  affect eo 
bv  the  occurrence  of  shear  layer  one  1 !_at ’ on . 
illustrates  the  mass  addition  process. 
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FIG.  7.  MASS  ADDITION  AND  REMOVAL  PROCESS  BY  SHEAR  LAYER 
* DEFLECTION  AT  THE  CAVITY  TRAILING  EDGE. 
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;•:<>?! ia a!  :-i:.  for  v it  1 s l-asr.  removal  3 n not  as  en « 1 1 y understood 
add  1 t \ on  . 1 J nonet  J in'  :*  said  that  when  the 

h.'.'ar  luv.'f  Is  deflected  urwar-d , the  fluid  non*1  the  t.ru  3 1J  nr. 

,1,-,.  * ;;  or.t  ra'ned  'Uf  of  'he  cavity.  "h!  • si  ator-oni  Is,  at 

1 c ! v : ! ‘ i d'":’;cu1.t  t.c  J ra--' tie  <:>  .tra  1 n-eth  at  t he 
eal  I tit.--  edre  K-inr  capable  of  reinovin.*  us  rue!:  mass  1 n a tiro 
o:t::  via  a ndd-*d  by  a downward  deflect  ion  in  n comparable  tine 

i a he  enure  t he  nans  added  durJnr.  downward  deflection 
rtera  a'  a print  1 vvl  y hirh  velocity  due  to  the  velocity  rradier.t 
;i  t > ..  laver*;  this  prosunab  1 *•■  has  a nuch  r.fronrer  effect 

*y. : . rorova’.  hv  entra* n-env . It  la  possible  that  waver  arrivinr 
.•  •-,,<>  bullhead  from  the  front  portion  o'*  1 he  cavitv  heln 

m;o  .-.a  s s oy*.  at  t!‘.e  trail  inr  edre.  Vhile  t Ml  s effect  is 
■.Ir-oro  carlalr.lv  present  , it  does  not  veer:  sufficiently  strong. 
;ert  urbut  i or  velocities  frotr.  these  waves  are  much  smaller  than 
he  freest ream  velocity,  furthermore,  the  types  of  overpressures 
'reduced  by  waves  nrrivinr  at  the  rear  bulkhead  seen  inadequate 
,0  achieve  the  required  level  of  trass  removal. 


The  correct  explanation  for  the  mass  removal  mechanist", 
lives  the  way  an  upward  shear  layer  deflection  affects  the 
:,al  steady-state  mass  addition.  Certainly,  when  the  shear 
*r  is  moved  unward,  the  normal  mass  addition  process  Is 
3t.inti.ally  reduced,  "'he  re  fore , the  "mss  removal"  process 
reailv  be  a lack  of  the  steady-state  mass  addition,  fince 
average  situation  in  the  cavity  involves  continuous  mass 
it  ion  at  the  trail  inr  ohpe  to  balance  the  shear  layer  entrain- 
a reduction  of  this  addition  acts  as  a mass  removal . 

: her  mo  re , this  "mass  removal"  effect  ear.  be  the  same  mar.nitudc 
i ’ne  extra  mass  addition  caused  by  downward  shear  laver  de- 
ction.  The  mass  removal  process  is  illustrated  in  i ~ . 7b. 
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ANALYSIS 


4.1  Introduction 

The  physical  mechanl  sirs  at  work  in  oscillatinr  cavities 
have  been  sufficiently  described  so  that  an  analytical  model 
for  the  cavity  wave  structure  can  now  bo  considered.  The  first 
two  sections  will  discuss  the  formulation  of  the  problem  and 
introduce  the  ideas  that  are  the  basis  for  the  simplified  model. 

The  following  sections  will  then  develop  the  analysis  in  detail 
and  discuss  the  results.  The  final  section  will  analyze  a 
semiempi  ri  cal  frequency  prediction  scheme  showinp  lenpth-to- 
depth  ratio  dependence. 

4.2  Analytical  Yodel  for  Cavity  Oscillation 
4.2.1  General  Considerations 

From  Sec.  3,  it  is  clear  that  the  cavity  wave  structure 
can  be  modeled  as  follows.  The  fluid  motion  within  the  cavity 
is  poverned  by  the  solution  of  the  wave  equation  in  a rep, ion 
bounded  by  a rlrid  front  bulkhead  and  floor,  with  a free  shear 
layer  on  top,  and  a rear  bulkhead,  which  oscillates  like  a 
piston.  This  oseillatinp  bulkhead  approximates  the  pseudopiston 
effect  of  mass  addition  and  removal  at  the  trailinp,  edp.e . The 
external  f low , which,  to  a pood  approximation,  is  poverned  by 
the  convective  wave  equation,  is  on  the  other  side  of  the 
shear  layer,  which  bepins  at  the  leadinp  edpe.  The  Internal 
and  external  solutions  should  be  .icined  by  annlvinp  the  proper 
dynamic  and  kinematic  conditions  across  the  shear  laver.  Tr 
the  cavity  leadinp  edpe  is  a sharp  corner-,  then  the  Kutta  condition 
should  be  applied  to  the  shear  layer  at  this  edp.e,  which  makes  its 
displacement  equal  to  zero  at  the  leadinp.  edpe.  Although  the 
shear  layer  motion  should  not  be  constrained  at  the  trailinp  edp.e, 
the  internal  and  external  repions  should  remain  separated.  The 
effect  of  the  trai linp-edpe  mass  addition  and  removal  process 
is  accounted  for  entirely  by  the  equivalent  oscillation  of  tire 
rear  bulkhead.  The  relationship  between  the  shear  layer  displace- 
ment at  the  trailinp  edpe  and  the  oscillation  of  the  rear  bulk- 
head is  considered  in  the  pseudopiston  analysis  in  Appendix  A. 


The  exact  shape  of  the  oseillatinp  rear  bulkhead  Is  unknown. 
We  need  not  assume  that  it.  neves  as  a ripid  wall,  since  it  is 
only  an  approximation  to  the  actual  trai linp-edpe  mass  addition 
ana  removal  process.  The  most  important  feature  of  its  motion 
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4.2.2  Real  Flow  Effects 
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As  mentioned  earlier,  t fie  rr-circulat  icr:  level  depends  not 
only  on  the  external  r 1 0..,-  Mach  number  and  hevnolds  number  and 
the  cavity  lenpt  h-t  c-dent  b.  ratio,  but  also  on  the  cavity  width 
and  the  thickness  and  state  (laminar  or  turbulent)  of  the  shear 
layer.  In  turn,  these  shear  layer  characteristics  depend  net 
only  on  external  flew  conditions,  but  also  on  the  thickness  and 
state  of  the  upstream  boundary  layer.  Although  the  recirculation 
velocities  car.  be  fairly  larre  , it  1 s believed  that  their  e ff  col- 
on the  internal  wave  motion  is.  secondary.  Mean  the  cavity  floor, 
the  recirculat ion  flow  will  move  forward;  while  near  the  shear 
layer,  where  entrainment  occurs,  the  recirculation  motion  will 
be  rearward.  The  net  effect  is  probably  a tiltinr  of  the  wave- 
fronts.  If  this  effect  acts  or.  an  upstream  wave  that  then  re- 
flects to  become  a downstream  wave,  the  lilt  direction  will  be 
reversed  in  the  (mirror  im.are)  reflection  process.  The  effect 
of  recirculation  flow  on  the  downstream  wave  will  then  be  to 


rotate  the  wave front  to  a more  uorimhi 
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.on.  !n  other  words, 


there  is  a certain  rierree  of  cancel  lat  ion  associated  with,  the 
effect  of  flow  recirculation.  It  seems  likely  that,  the  effect 
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pseudopiston  due  to  the  mass  addition  process). 


Another  concern  introduced  by  real  flow  effects  is  the  prone: 
treatment  of  the  shear  layer  dynamics.  The  simples?  approach  is 
to  treat  the  shear  layer  as  an  infinitely  thin,  vort  ex  sheet. , 


i.e.,  an  idealised  veloc: 
the  freestream  velocity  : 
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■i'.'.'.  i i fi  cat  ior. , however,  the  shear  layer  thickness  and 
r.ro  r.oi  easily  predicted.  Tt  is  believed  that  the 
shear  layer  t sickness  produces  a r.lrnJfirar.t  offer*  on  the 
cavity  dvnnt: ! c l\v  alterinr  t ho  linonal  ic  and  dynntrJe  conditions 
t hat  match  the  internal  and  external  sol  unions . This  effect 
would  bo  rou.’hly  si:-. liar  to  the  inclusion  of  an  effective  inertia 
due  to  i he  shear  Inver  thickness.  If  the  spatial  variation  of 
shear  Inver  thickness  is  not  too  j treat  ever  the  wavolonrth,  then 
the  probier.  oar.  be  nanayed  analytically  1 •’  the  shear  layer 
structure  is  ronoled  i u a simple  : .ar.r.er . 'the  analytical  re- 
st; 1 1. s indicate  that  the  inclusion  of  shear  layer  thickness 
effects  (which  will  reduce  shear  layer  anpl if i cation  rates) 
is  an  important  topic,  which  is  worthy  of  future  theoretical 
inv on  t ir.at  ion  . 

4.2.3  A Simple  Analytical  Model 

-nr.ee  ’ho  important  features  of  an  analytical  nodel  are 
reviewed,  it  is  appropriate  t o discuss  the  actual  solution  of 
tin-  p rob  lev.  or  an  accept  able  appro  >:1r:at  e solution.  The  ext  err. al 
and  internal  re, "ions  can  both  be  solved  ty  transform  methods 
with,  the  application  of  appropriate  boundary  conditions.  The 


two  solutions  must  be  matched  kinematically  and  dynamically  across 
the  shear  layer.  If  the  shear  layer  is  treated  as  an  infinitely 
thin  discontinuity,  then  this  matching  corresponds  to  applying 
the  material  surface  condition  and  requiring,  equal  static  pres- 
sures. If  shear  layer  thickness  effects  are  included,  then  a 
more  complicated  matching,  condition  must  be  applied.  The 
physical  solution  is  then  obtained  by  inverse  transforms. 

rphis  formal  approach  is  subject  to  severe  mathematical 
difficulties.  Indeed,  the  effort  required  may  not  be  justifiable 
on  physical  grounds.  Because  the  details  of  the.  trail ing-edge 
mass  addition  process  cannot  be  determined,  any  effort  SDent 
assuring  that  other  cavity  boundary  conditions  are  satisfied 
in  detail  Is  wasted.  Specifically,  a different  set  of  mass 
flow  moments  at  the  rear  bulkhead  will  change  the  detailed 
wave  structure  elsewhere;  hence,  the  strength  of  waves  required 
to  satisfy  the  boundary  conditions  will  be  altered.  Therefore, 

It  is  proposed  to  use  a simple  physically  motivated  solution 
that  satisfies  only  the  most  basic  conditions  and  that  does  not 
include  the  finer  details.  In  fact,  this  is  the  onlv  approach 
that  is  really  justifiable.  When  the  cavity  problem  is  reduced 
to  its  essence,  only  the  basic  upstream  and  basic  downstream 
waves  are  considered.  The  reflection  condition  at  the  front 
bulkhead  will  be  satisfied  only  to  the  extent  that  the  instantaneous 
net  mass  flow  induced  by  the  waves  be  zero;  moments  of  mass  flow 
are  not  required  to  vanish.  This  approach  is  consistent  with  the 
rear  bulkhead  condition,  where  there  is  an  oscillating  net  mass 
flow.  The  Kutta  condition  will  be  applied  at  the  leading  edge 
only  in  the  sense  that  the  spatial  envelope  of  the  waves  decays 
exponentially  upstream.  This  approach  allows  the  problem  to 
be  reduced  to  the  motion  of  waves  between  a shear  layer  and  a 
boundary.  The  frequency  of  these  waves  is  real  corresponding 
to  an  oscillating  phenomenon;  therefore,  their  wave  number  is 
complex,  corresponding  to  propagation  and  spatial  amplification 
or  decay.  Since  only  basic  waves  are  considered,  and  only  one 
condition  (front  bulkhead  reflection)  is  required  to  relate  the 
wave  strengths,  transform  techniques  are  no  longer  required. 

At  this  point,  the  oscillation  frequency  of  the  rear  bulk- 
head or  pseudopistcn  is  still  undetermined.  Determination  or 
this  frequency  may  come  through  a set  or  phase  and  amplitude- 
matching  conditions,  which  are  arrived  at  by  separate  considera- 
tion of  the  trailing-Cidge  mass  addition  process  and  then  the  syn- 
thesis of  the  two  analyses.  For  instance,  a simple  condition  would 
require  the  shear  layer  motion  to  be  in  phase  with  the  oscillating 
rear  bulkhead.  However,  because  of  the  complexity  of  the  processes 
at  the  trailing  edge,  especially  in  view  of  real  flow  effects,  it 
Is  easy  to  believe  that  significant  phase  shifts  between  the  shear 
layer  and  pseudopiston  may  occur. 
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Given  this  latitude,  it  is  possible  to  suppose,  alternatively, 
that  within  some  limits  the  phase  and/or  amplitude  of  the  mass 
addition  process  adjusts  to  a condition  most  favored  by  the  in- 
ternal wave  dynamics  of  the  cavity.  This  adjustment  might  serve, 
for  instance,  to  accommodate  a resonance  of  the  system.  Or,  as 
another  possibility,  the  condition  might  be  the  frequency  that 
gives  the  maximum  shear  layer  displacement  for  a given  pseudo- 
piston displacement,  thereby  maximizing  the  energy  removal  from 
the  mean  flow  for  a given  level  of  cavity  pressure  oscillations. 
This  last  condition  seems  to  be  the  correct  one  in  that,  in  later 
sections,  it  will  be  shown  to  give  the  best  agreement  between  the 
present  theory  and  the  experimental  measurements.  Such  conditions 
are  actually  satisfied  for  a set  of  frequencies,  given  the  cavity 
geometry.  Each  frequency  corresponds  to  a different  mode  of 
oscillation.  Experiments  show  that  these  modes  often  occur 
simultaneously.  The  cavity  geometry  and  the  flow  Mach  number 
determine  which  mode  level  dominates;  however,  it  is  frequently 
the  second  mode.  Presumably,  this  dominance  can  be  determined 
theoretically.  It  may  correspond  to  which  mode  can  extract  the 
most  mean  flow  energy  or  to  which  mode  produces  the  lowest  im- 
pedance of  the  pseudopiston. 

4.3  Wave  Motion  of  a Shear  Layer  Above  a Boundary 

In  this  section,  we  discuss  the  wave  motion  between  a shear 
layer  and  a boundary.  The  ultimate  goal  Is  to  understand  the 
speed  and  amplification  characteristics  of  the  type  of  waves 
that  are  responsible  for  the  cavity  oscillation  phenomenon. 

The  basic  solutions  developed. here 'will  eventually  be  super- 
imposed to  represent  the  boundaries  of  a cavity  being  forced 
by  a piston  at  the  rear  bulkhead.  In  this  manner,  the  modal 
structure  within  the  cavity  can  be  established. 

Consider  a shear  layer  and  a parallel  solid  boundary 
separated  by  a distance,  D.  Between  the  shear  layer  and  the 
boundary  is  an  inviscld,  compressible  fluid  with  a sound  speed, 
ac.  On  the  other  side  of  the  shear  layer  is  an  Inviscld, 
compressible  fluid  moving  with  uniform  velocity,  U<»,  and 
having  a sound  speed,  a®.  It  is  of  interest  to  find  the  re- 
lation that  governs  the  motion  of  small  amplitude  waves  under 
the  shear  layer  (see  Fig.  8). 

The  wave  motion  within  the  cavity  is  analyzed  using  the 
coordinate  system  (I)  in  Fig.  8.  The  perturbation  velocity 
potential  is  governed  by  the  simple  wave  equation  in  two 
dimensions : 
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FIG.  8.  WAVES  ON  A SHEAR  LAYER  ABOVE  A BOUNDARY. 
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Assume  solutions  of  the  form. 
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where 
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(4.3) 


A more  general  solution  for  fixed  frequency  can  be  constructed 
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where  Cn  represent  arbitrary  constants.  .Similarly,  solutions 
with  different  frequencies  could  be  added  if  desired.  Notice 
that  the  condition  of  zero  normal  velocity  at  the  boundary  is 
satisfied  by  using  the  hyperbolic  function  to  re-express  the 
y-dependence . 


The  perturbation  pressure  at  the  shear  layer  is 
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The  displacement  of  the  shear  layer  is 
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The  flow  phenomena  above  the  shear  layer  Is  governed  by  the 
convective  wave  equation: 
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The  coordinate  system  (II)  in  Pig.  8 is  now  used.  Assume 
solutions  of  the  form 
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As  before,  these  solutions  can  be  added  to  achieve  a more  general 
form, 
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where  the  An  are  arbitrary  constants.  The  perturbation  pressure 
at  the  shear  layer  is  now  given  by 


(- 


, - U $ 
t “Tx 


^ y=0+ 


. . m Ik  x 

= eltot  2 _ i (cu+U  k )A  e n 
n=l  00  n n 


(4.11) 


The  normal  velocity  at  the  shear  layer  is 
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We  can  apply  the  boundary  conditions  that  static  pressures 
match  across  the  shear  layer  and  that  the  shear  layer  be  a 
material  surface.  The  pressure-matching  condition  requires 
that  the  two  pressures  found  by  considering  the  regions  above 
and  below  the  shear  layer  should  be  equated.  The  material 
surface  condition  is  determined  by  requiring  that  the  vertical 
velocity  in  the  upper  region  calculated  from  shear  layer  displace- 
ment equals  that  calculated  directly: 


V ( —t>  \ = rn  + U n 

[ y y=o+ / " x 


(4.13) 


The  first  term  on  the  right-hand  side  arises  from  unsteady  shear 
layer  motion  and  the  second  from  tilting  the  free  stream  so  it 
is  tangent  to  the  displaced  shear  layer.  The  completion  of  the 
details  leads  to  the  following  dispersion  relation  between 
frequency  and  wavenumber: 


it u2B  coshaD  = — (to  + U k)2  a sinhaD  , 
pc 


where 
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For  disturbances  that  originate  in  the  vicinity  of  the 
shear  layer  and  the  boundary,  only  some  of  the  solutions  of 
the  dispersion  relation  will  be  admissible.  A radiation 
condition  will  be  used  to  exclude  solutions  that  represent 
waves  traveling  inward  from  infinity  and  solutions  that  become 
singular  at  infinity.  Furthermore,  the  case  of  interest  will 
be  tne  determination  of  wavenumber  k for  a fixed  real  frequency, 
to.  This  corresponds  to  the  case  of  a forced  system  that  is  of 
interest  for  the  future  application  of  these  solutions.  When 
the  frequency  is  real,  the  wavenumber  will  generally  be  complex. 
This  corresponds  to  waves  that  travel  and  amplify  or  decay 
spatially . 


f. 


Because  of  the  transcendental  nature  of  the  dispersion 
relation>  there  are  infinitely  many  solutions  for  a fixed 
frequency.  Those  that  are  admissible  by  the  radiation  condition 
can  be  divided  according  to  whether  they  travel  upstream  or 
downstream.  The  main  interest  will  be  in  the  speed  and  amplifica- 
tion of  waves  that  correspond  most  nearly  to  plane  waves 
traveling  between  the  shear  layer  and  the  boundary.  It  is  these 
waves  that  seem  to  dominate  in  the  oscillation  nhenomenon  of 
shallow  cavities. 

The  dispersion  relation  will  now  be  re-expressed  in  a more 
appropriate  form.  If  we  use  the  equation  of  state  and  the 
definition  of  the  speed  of  sound,  and  assume  equal  static  pres- 
sures prior  to  perturbation,  we  obtain 


If  the  ratio  of  specific  heats  is  the  same  on  both  sides  of 
the  shear  layer,  then 


(4.16) 


If  we  assume  that  the  static  temperature  below  the  shear  layer 
equals  the  stagnation  temperature  above  (the  appropriate 
assumption  for  the  application  of  this  analysis  to  cavity 
oscillations),  then 
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where  y - 1 . 4 for  air. 

It  is  convenient  to  define  the  following  nondimensional 
quantities : 
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K = kD 

(4.18) 

T = aD  and  T2  = K2  - S2 
Q = (MK  + aS)2 
B = 6D  and  B2  = Q-K2 

Using  the  above,  the  dispersion  relation  can  be  manipulated  into 
the  following  form. 


Q2T2sinh2T  + S4 (Q-K2 )cosh2T  = 0 


(4.19) 


In  the  process  of  obtaining  this  form,  the  original  relation  has 
been  squared.  Although  the  order  is  now  doubled,  the  ambiguity 
of  signs  for  the  square  roots  is  no  longer  present  and  the  total 
number  of  roots  is  unchanged. 

Solutions  of  the  above  form  of  the  dispersion  relation  were 
obtained  using  the  computer.  Analytical  techniques  that  were 
also  developed  were  either  unsuccessful  or  not  sufficiently 
accurate  over  a wide  enough  range  of  the  parameters  to  justify 
their  further  use.  The  general  spatial  configuration  of  roots 
in  the  complex  K plane  is  Illustrated  In  Pig.  9.  A detailed 
map  of  ro-  ■> , A,  B,  and  C,  is  presented  In  Pig.  10.  The  roots 
D and  E faal  on  the  negative  real  K axis;  their  presentation 
in  Pig.  10  would  not  be  Informative  since  everything  falls  on 
a single  line.  The  additional  roots  indicated  in  Pig.  9 were 
computed.  However,  such  a computation  follows  easily,  either  by 
numerical  methods  or  by  assumption  of  large  |K|  in  the  dispersion 
relation  in  order  to  achieve  an  analytical  approximation.  These 
roots  would  be  useful  only  if  the  Internal  wave  structure  of 
the  cavity  were  to  be  constructed  in  greater  detail  than  is  the 
present  intention.  Appendix  B presents  all  of  the  roots  that 
were  computed. 

The  next  step  Is  to  choose  the  waves  that  satisfy  the 
boundary  conditions  of  the  problem,  so  that  they  can  be  super- 
imposed to  (approximately)  form  the  cavity  peometrv.  ^hen. 
the  detailed  dynamics  of  the  system  can  be  investigated.  Upstream 
traveling  waves  h^ve  ReK  > 0,  and  downstream  waves  have  ReK  < 0. 
The  condition  th-  aisturbances  do  not  amplify  away  from  the 
boundary,  namely  chat  amplitudes  are  finite  or  zero  as  y + 
requires  ImB  >_  0 . It  Is  also  necessary  that  the  spatial  envelooe 
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of  the  waves  decays  upstream  since  superfluous  upstream  influence 
in  this  simple  cavity  model  would  then  be  reduced,  furthermore } 
this  requirement;  is  consistent  with  the  idea  that  the  upstream 
wave  radiates  energy  anc.  decays,  whereas  the  downstream  wave 
draws  energy  from  mean  flow  and  amplifies.  To  satisfy  these  condi- 
tions, the  roots  K*  (the  complex  conjugates)  should  be  used. 

The  roots  A*,  B# , and  C*  are  the  prime  candidates  for  use 
in  a simple  cavity  oscillation  model.  3oot  C*  is  the  only 
upstream  traveling  wave,  so  it  must  be  included  in  the  analyt- 
ical model.  Both  A*  and  B#  are  downstream  waves,  and  It 
must  be  determined  which  is  the  important  wave  in  the  cavity 
oscillation  phenomena.  If  we  look  at  the  nondimensional  phase 
speed  of  these  waves  at  the  highest  Mach  number  (see  M=3  in 
the  table),  it  is  apparent  that  root  A*  behaves  as  a pressure 
wave  traveling  below  the  shear  layer  at  nearly  the  speed  of 
sound,  while  root  B*  behaves  as  a convection  wave  on  the  shear 
layer  moving  at  nearly  half  the  mean  flow  speed  (supersonically 
relative  to  the  cavity  interior  when  M=3).  Note  that  as  M-+® 
the  shear  layer  becomes  infinitely  stiff,  and  the  pressure  waves 
below  the  shear  layer  behave  as  if  in  a hard-walled  enclosure 
and  travel  at  the  speed  of  sound.  Therefore,  it  is  concluded 
that  C#  is  the  relevant  upstream  wave  and  A*  is  the  relevant 
downstream  wave. 

4.4  Wave  Solution  Synthesi  zatiort 

Solutions  for  wave  motion  between  a shear  layer  and  a 
boundary  have  been  found  for  the  case  of  real  frequency.  Two 
of  these  waves  will  be  combined  to  approximate  a cavity  with 
a piston  at  the  rear  wall,  as  shown  in  fig.  11.  ^he  wave 
solutions  must  satisfy  the  following  conditions; 

1.  They  must  be  essential  to  the  phenomenon;  i.e.,  they 
must  represent  the  basic  physical  effects  and  not 

be  needed  only  to  provide  higher  order  structure. 

2.  Both  upstream  and  downstream  waves  must  be  used. 

3.  They  must  satisfy  the  boundary  conditions  at  infinity: 
i.e.,  these  waves  must  originate  in  the  region  of 

the  shear  layer  and  boundary  and  not  at  infinity. 

4.  They  must  decay  exponentially  in  the  upstream  direction. 
This  reduces  the  upstream  influence  in  this  simple  model 
and,  in  effect,  provides  an  approximate  Kutta  condition. 

In  the  previous  section.  It  was  found  that  the  waves 
designated  as  A*  and  C*  satisfy  these  requirements. 
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The  upstream  wave  C#  and  the  downstream  wave  A#  can  be 
combined  to  form  an  approximation  to  wave  motion  in  an 
oscillating  cavity.  With  only  two  waves,  .lust,  one  boundary 
condition  at  the  cavity  leading  edge  can  be  aonlied.  mhis 
boundary  condition  will  be  that  the  instantaneous  net  mass 
flow  through  the  front  bulkhead  is  always  zero.  Because  these 
two  waves  do  not  have  the  same  phase  dependence  with  respect 
to  depth,  their  local  x-components  of  velocity  cannot  cancel 
over  the  entire  front  bulkhead.  Thus,  the  instantaneous 
moments  of  mass  flow  through  the  bulkhead  vrill  generally  not 
vanish.  It  is  felt  that  these  nonzero  moments  of  mass  flow 
only  change  the  details  of  the  cavity  wave  structure,  and  they 
do  not  produce  any  alteration  of  the  basic  process. 


Similarly,  the  motion  of  a piston  at  the  cavity  trailing 
edge  is  defined  in  terms  of  the  net  mass  flow  at  the  rear 
bulkhead  position.  Again,  moments  of  mass  flow  will  not 
vanish.  This  approach  is  Justified  because  the  actual  mass 
addition  process  does  not  correspond  to  a rigid  piston  shape. 


Within  the  cavity,  the  velocity  potential  of  this  two-wave 
system  is 
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where  A and  C are  constants.  The  velocity  components  are 
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The  shear  layer  deflection  is  derived  from  the  vertical  velocity 
component : 
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The  pressure  is : 
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The  front  bulkhead  reflection  condition  requires  that  the 
instantaneous  net  mass  flow  vanish,  namely 

D ...  D 
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Thus,  the  relationship  between  the  constants  is 
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Similarly,  the  depth  average  of  the  lengthwise  velocity  component 
is 
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The  depth  average  pressure  is 
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When  the  depth-averaged  velocity  and  pressure  are  evaluated  at 
x = L,  they  become  the  mean  velocity  and  mean  pressure  on  the 
pseudopiston . 


It  is  convenient  to  recast  these  results  in  nondimensional 
form.  As  in  Sec.  4.3,  let 


K = kD,  S = , and  T = ctD  . 

ac 


Recall  that  the  quantities  K and  T are  functions  of  and 


S.  Redefine  the  constants  as 
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The  depth-averaged  Mach  number  of  the  wave-induced  longitudinal 
velocity  perturbation  is 
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The  nondimenslonal  depth-averaged  pressure  is 
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and  the  nondimenslonal  shear  layer  displacement  is  defined  as 
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If  we  eliminate  one  of  the  constants  and  nondlmensionalise 
the  results,  we  obtain 
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There  are  two  other  quantities  of  interest.  The  complex 
acoustic  impedance  of  the  pseudopiston  is  given  by 
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The  complex  ratio  of  n/Mm  is  also  of  interest,  since  its  phase 
indicates  the  relation  between  shear  layer  position  at  the 
trailing  edge  and  the  pseudopiston  motion.  The  phase  angle  is 
zero  when  the  shear  layer  deflection  becomes  negative  Just  as 
the  pseudopiston  begins  to  move  into  the  cavity.  The  magnitude 
of  this  quantity  is  the  reciprocal  of  piston  motion  per  unit 
shear  layer  deflection. 


These  formulas  can  now  be  used  with  the  results  of  the 
previous  section  to  calculate  the  wave  structure  in  oscillating 
cavities . 

4.5  Discussion  of  Results 


Using  the  analyses  of  the  previous  sections,  calculations 
were  carried  out  to  predict  pressure-mode  shapes  and  other 
properties  of  interest  in  shallow  cavities.  For  mathematical 
convenience,  the  calculations  were  carried  out  with  nondimensional 
frequency,  S,  held  constant,  while  the  distance  from  the  leading 
edge  was  varied.  In  this  way  the  mode  shapes  are  immediately 
obtained.  A certain  discrete  set  of  distances  from  the  leading 
edge  corresponds  to  the  length-to-depth  ratios  of  cavities  that 
will  oscillate  freely  at  frequency  S.  When  the  calculation  is 
repeated  for  a different  nondimensional  frequency,  a different 
set  of  length-to-depth  ratios  will  be  found.  Thus,  the 
dependence  of  frequency  on  mode  number,  length-to-depth  ratio, 
and  flow  conditions  can  be  generated  through  this  calculation 
procedure.  Other  methods,  such  as  fixing  length-to-depth  ratio 
and  varying  frequency,  could  also  be  used.  It  should  be  noted 
here  that  the  nondimensional  frequency,  S = toD/ac,  is  not  the 
same  as  S*  = fL/Uo,,  the  Strouhal  frequency  which  is  often  used 
to  display  cavity  oscillation  results.  Of  course,  they  are 
related : 
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where  a 


A set  of  results  was  computed  for  the  case  of  = 3.0  and 
S = 1.0.  Although  this  Mach  number  is  really  above  the  range 
of  greatest  interest,  the  results  are  instructive  for  this  case; 
furthermore,  a difficulty  that  occurs  at  lower  Mach  numbers  is 
avoided.  Figure  12  shows  the  magnitude  of  unsteady  pressure 
fluctuations,  expressed  in  dB,  as  a function  of  distance  from 
the  leading  edge.  The  arrows  denote  the  possible  locations  of 
a rear  bulkhead.  Although  no  experimental  mode  share  data  are 
available  at  Mro  = 3.0,  it  is  encouraging  that  this  computed  shape 
bears  a considerable  qualitative  resemblance  to  mode  shaoe 
measurements  at  lower  Mach  numbers.  Figures  13  and  lb  show 
the  magnitude  of  longitudinal  velocity  and  shear  layer  displace- 
ment, respectively.  Since  both  magnitudes  are  only  determined 
to  within  a constant  multiple,  these  curves  should  be  of  interest 
primarily  for  their  shape  and  relative  size. 

~ Figures  15  and  16  show  the  magnitude  and  phase  of  the  ratio 
n/Mm,  which  is  the  complex  shear  layer  deflection  per  pseudo- 
piston velocity.  Here,  each  value  of  x/D  is  interpreted  as  L/D, 
i.e.,  as  a candidate  trailing-edge  location.  Similarly,  Figs.  17 
and  18  show  the  magnitude  and  phase  of  the  pseudopiston  impedance. 
Finally,  Figs.  19^and  20  show  the  magnitude  and  phase  of  the 
complex  ratio  p^/n,  which  is  a type  of  impedance  based  on  shear 
layer  displacement  rather  than  pseudopiston  velocity.  These 
quantities  are  presented  as  function  of  x/D  (or  L/D)  at  constant 
frequency  S.  (Note  that  this  manner  of  presentation  differs 
from  the  more  familiar  impedance  vs  frequency  often  seen  in 
acoustics . ) 

The  important  question  arises  as  to  how  the  trailing-edge 
locations  (L/D)  can  be  determined  from  the  information  presented. 
Experimental  results  show  that  a pressure  maximum  occurs  near, 
but  not  exactly  at,  the  trailing  edge;  very  often  the  pressure 
maximum  occurs  just  before  the  trailing  edge  (see  Sec.  6).  Of 
course,  the  pressure  maxima  correspond  closely  to  the  noints  of 
maximum  impedance.  The  vicinity  of  the  pressure  maxima  seems  to 
bear  no  special  relation  to  the  phase  of  the  shear  layer  relative 
to  the  pseudopiston,  arg  r/ffi_  (see  Fig.  16).  It  is  somewhat 
surprising  that  this  phase  difference  is  never  zero,  which  would 
correspond  to  the  intuitive  kinematic  relation  between  shear 
layer  position  and  pseudopiston  motion  for  an  oscillating  cavity. 
However,  as  mentioned  earlier,  there  is  reason  to  exreet  that 
significant  phase  shifts  are  introduced  by  the  actual  mass 
addition  process.  The  one  quantity  that  behaves  in  a distinct 
manner  near  the  qualitatively  expected  trailing-edge  locations 
is  | n/ftm| . This  quantity,  illustrated  in  Fig.  15,  shows  a 
series  of  distinct  and  rather  narrow  peaks  curring  slightly 
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behind  the  pressure  maxima.  The  rest  of  the  curve  is  relatively 
flat.  It  is  believed  that  the  maxima  of  |n/$m|  may  represent 
the  proper  trailing-edge  locations.  The  arrows  on  Fig.  12  were 
located  according  to  this  criteria.  The  maxima  of  |n/$m| 
correspond  to  the  most  mean  flow  addition  for  a given  pseudo- 
piston displacement.  Therefore,  this  condition  corresponds  to 
a maximum  removal  of  energy  from  the  mean  flow  by  the  trailing- 
edge  mass  addition  process  for  a given  motion  of  the  pseudo- 
piston. If  the  net  energy  losses  from  other  processes  in  the 
cavity  are  determined  essentially  by  its  level  of  activity  (as 
indicated  by  the  pseudopiston  motion) , then  these  criteria  provide 
the  most  addition  of  energy  to  overcome  these  losses. 

A somewhat  different  viewpoint  can  bring  us  to  the  same 
conclusion.  From  Fig.  13,  it  is  apparent  that  the  magnitude 
of  wave-induced  longitudinal  velocity  in  the  cavity  is  never 
zero,  except  at  the  leading  edge.  Therefore,  without  the  trailing- 
edge  mass  addition  and  removal  process,  the  proper  boundary 
condition  could  never  be  satisfied  at  the  rear  bulkhead.  However, 
it  is  nearly  satisfied  in  the  vicinity  of  the  velocity  minima, 
since  in  this  location  the  least  pseudopiston  notion  is  required 
to  provide  the  proper  boundary  condition.  Clearly,  the  trailing- 
edge  mass  addition  process  operates  most  efficiently  when  the 
maximum  shear  layer  deflection  is  available  for  a given  pseudo- 
piston motion.  Thus,  the  boundary  condition  at  the  trailing 
edge  is  most  easily  satisfied  whenever  |n/Mm|  Is  a maximum. 

Once  the  criterion  for  picking  the  trailing-edge  location 
is  established,  the  oscillation  frequencies  for  a fixed  cavity 
geometry  can  be  determined.  Then,  it  may  be  possible  to  determine 
the  relative  strength  of  these  modes  of  oscillation.  It  is 
believed  that  the  determination  of  the  dominant  mode  is  related 
to  a condition  such  as:  which  mode  gives  the  lowest  impedance 

of  the  pseudopiston;  or  which  mode  removes  the  most  energy  from 
the  mean  flow;  or  for  which  mode  does  the  mass  addition  process 
operate  most  efficiently.  This  important  area  of  research  should 
receive  more  attention  in  the  future,  particularly  since  the 
answer  has  an  Importance  which  extends  beyond  the  present  problem. 

The  present  analysis  is  less  successful  when  the  calculation 
Is  made  for  lower  Mach  number.  This  problem  Is  illustrated  in 
Fig.  21,  where  the  Mach  number  is  = 1.0  and  S = 1/2.  It 
can  be  seen  that  only  the  first  mode  appears,  and  that  the  magnitude 
of  unsteady  pressure  Increases  rapidly  with  increasing  x/D.  A 
similar  trend  is  found  when  the  calculation  is  repeated  for 
Mm  = 1.0  and  S = 1.0.  Although  the  first  mode  is  predicted 
successfully,  none  of  the  higher  modes  is  found  because  the 
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magnitude  of  pressure  increases  rapidly  without  additional  nodes. 
This  result  contradicts  the  experiments  where  modes  as  high  as 
the  fourth  are  detected.  In  fact,  the  second  mode  is  often  the 
one  with  dominant  pressure  levels.  Because  this  shortcoming  is 
found  everywhere  below  about  M = 2.0  for  all  S,  an  extensive 
investigation  was  undertaken  to  identify  the  problem.  During 
this  investigation,  several  alternate  boundary  conditions  at 
both  leading  and  trailing  edges  were  tried.  The  effect  of  in- 
cluding the  shear  layer  convection  wave  B*  was  also  considered 
for  several  different  boundary  conditions.  For  instance,  when 
this  additional  wave  Is  Included,  the  Kutta  condition  can  be 
applied  in  addition  to  the  reflection  condition  at  the  front 
bulkhead.  Finally,  the  possible  effects  of  adding  some  of  the 
other  waves,  and  of  doing  the  problem  more  formally,  were 
considered.  No  dramatic  improvements  over  the  basic  model  were 
made,  and  some  results  were  much  worse.  Nor  did  it  seem  worth- 
while to  make  the  problem  more  complicated. 

Careful  investigation  shows  that  the  source  of  the  problem 
is  that  the  spatial  amplification  rate  of  the  downstream  wave 
greatly  exceeds  tnat.  of  the  upstream  wave  over  much  of  the  Mach 
number  range.  This  fact  is  apparent  from  the  basic  wave  solu- 
tions which  show  that  the  imaginary  part  of  K*  exceeds  that  of  ; 
the  discrepancy  increases  with  decreasing  Mach  number.  Not  far 
behind  the  leading  edge,  the  exponential  amplitude  growth  of 
the  downstream  wave  dominates  so  completely  that  pressure  nodes 
can  no  longer  be  formed;  i.e.,  the  energy  extracted  by  the  down- 
stream wave  from  the  flow  greatly  exceeds  that  radiated  by  the 
upstream  wave,  ^his  fact  does  add  theoretical  support  to  the 
observation  that  cavity  oscillations  are  self-sustaining  under 
almost  all  flow  conditions.  However,  since  higher  modes  are 
commonly  observed  experimentally,  it  must  be  concluded  that 
either  the  amplification  rate  of  the  shear  layer  is  overestimated 
or  that  the  ultimate  amplification  of  the  downstream  wave  is 
limited  by  some  additional  process.  An  example  of  the  second 
possibility  would  be  the  nonlinearity  introduced  by  the  shear 
layer  rolling  up  into  discrete  vortices,  which  has  been  observed 
at  subsonic  flow  speeds.  The  present  analytical  difficulty 
extends  well  Into  the  supersonic  flow  ime  where  roll  up  is 
much  less  likely. 

It  is  believed,  therefore,  that  this  simple  model,  which 
assumes  an  infinitely  thin  shear  layer,  overestimates  the 
amplification  rate  of  the  downstream  wave.  This  difficulty 
can  be  resolved  by  accounting  for  the  finite  thickness  and 
different  flow  velocity  of  the  shear  layer.  This  effect  could 
be  included  analytically  in  a future  research  effort.  The 
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->ed  for  the  cavity  problem  does  not  involve 
i t leal  flow  solutions  in  the  shear  layer  region. 

layer  is  properly  accounted  for,  it  is  expected 
;r  results  could  be  obtained  at  lower-  Mach  numbers, 
/ the  same  analytical  approach  as  described.  If 
rh  Mach  number  results  are  any  indication,  this 
cal  approach  may  eventually  lead  to  a complete  and 
ctical  understanding  of  the  physics  of  the  oscilla- 


4.6  Frequency  Prediction 

In  this  section,  the  wave  solutions  found  in  Sec.  4.2  are 
used  as  the  basis  for  a semiempirical  frequency  prediction 
scheme.  The  approach  taken  is  similar  to  that  of  Rossiter 
(1966).  This  discussion  helps  to  illustrate  the  relationship 
between  Rossiter' s assumptions  and  the  present  physical  model. 

The  pressure  modes  in  the  cavity  are  made  from  the  super- 
position of  upstream  and  downstream  waves  whose  phase  speeds 
are  different.  Suppose  the  averap;e  distance  between  pressure 
maxima  is  X . Because  of  the  different  phase  speeds  and 
amplification  rates,  the  actual  distances  between  pressure 
maxima  need  not  all  be  the  same.  Let  L + = nXp>  where  n is 

an  integer  and  L is  the  cavity  length.  The  quantity  6 is  an 
undetermined  value  that  accounts  for  the  facts  that  the  cavity 
trailing  edge  is  not  actually  a reflection  plane  (pressure  maxima), 
that  there  are  phase  shifts  in  the  system,  and  that  Ap  is  only  an 
average  value. 

If  T is  the  oscillation  period. 


-L+-6- . + L±i_  = nT 

c c 


( 4 . 40) 


where  the  mode  number  n equals  the  number  of  pressure  minima  in 
the  cavity,  Cn  and  Cn  are  the  real  part  of  the  phase  speed 
1 Ar  HCr 

of  the  downstream  wave  and  the  upstream  wave,  respectively. 


(4.41) 
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Suppose  the  factor  6 is  a fraction  (positive  or  negative)  of  a 
wavelength,  namely  6 - aX  . Since  1 + 6 = nX  , and  L = nX  -6: 

P P P 


L+6 

L 


n 

n-a 


(4.45) 


Hence 


'kaP 
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Kc  * * \ _ 2ir(n-q) 


K 


Cr 
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(4.46) 


where  and  Kc  are  functions  of  Mm  and  S.  Tf  !*!„,  n,  and  L/D 
are  specified,  and  a value  of  the  empirical  constant  a js  assigned, 
then  the  above  equation  determines  S.  From  this,  the  cavity 
Strouhal  number,  S*  = fL/Uw,  can  be  determined. 


S* 


- a(L/D) 
2ttM 


(4.47) 


T'*"  yrr  -wy. 
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Figure  22  shows  these  results  plotted  for  a = 0.1  for 
n = 1.0,  and  a = 0.25  for  n * 2,  3,  and  4.  The  Rossiter  formula 
is  also  shown. 

The  results  show  the  effect  of  length-to-depth  ratio.  The 
predicted  Strouhal  numbers  begin  to  fall  below  the  measured 
results  when  M < 1.0.  This  problem  can  be  traced  once  again  to 
the  excessive  amplification  rate  of  the  downstream  wave,  which 
affects  the  phase  speed.  If  the  calculation  is  repeated  by 
setting  the  imaginary  part  of  K equal  to  zero  (no  amplification), 
then  the  computed  frequencies  fall  above  the  measured  results  at 
low  Mach  number.  This  fact  suggests  that  lowering  the  amplifica- 
tion rate  will  raise  the  present  results  at  low  Mach  number, 
which  then  gives  better  agreement  with  the  experiment.  It  is 
noted  that  there  is  no  guarantee  that  the  empirical  constant  a 
is  independent  of  Mach  number.  However,  substantial  variations 
in  a must  be  introduced  to  account  for  the  discrepancy  at  low 
Mach  number. 
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EXPERIMENTAL  CONSIDERATIONS 
5.1  Scope  of  Experimental  Studies 

The  analytical  considerations,  presented  in  the  previous 
sections,  provided  insight  into  the  basic  physical  phenomena 
that  ar  responsible  for  the  occurrence  of  flow-induced  pres- 
sure oscillations  in  shallow  cavities.  In  particular,  it  was 
determined  that  the  onset  of  pressure  fluctuations  depends  on 
the  balance  between  energy  supplied  by  the  external  flow  and 
energy  dissipated  by  viscous  losses  and  acoustic  radiation, 
so  that  an  excess  of  energy  drawn  from  the  external  flow 
results  in  the  onset  and  sustenance  of  oscillations  and  vice 
versa. 

While  the  oscillation  onset  is  linked  to  both  phase 
and  gain  criteria,  no  particular  onset  velocity  was  found 
to  exist  for  shallow  cavities.  However,  a set  of  general 
aerodynamic  criteria  was  established  that  define  a region 
(in  terms  of  geometry  and  flow  parameters)  where  the  least 
damping  of  shear  layer  wave  motion  occurs;  thus,  the  oscilla- 
tion process  sustains  itself.  The  analysis  had  a twofold 
purpose:  (1)  to  determine  the  oscillation  frequencies,  the 

mode  shapes,  and  mode  levels  in  a resonant  cavitv;  and  (2)  to 
establish  criteria  for  stabilizing  the  flow  in  order  to  suppress 
the  oscillation.  The  analysis  was  reasonably  successful  since 
it  provided  resonant  frequencies  as  a function  of  geometric  and 
aerodynamic  parameters  and,  to  some  extent,  the  mode  shapes; 
however,  because  of  the  extreme  complexity  of  the  flow/cavity 
interaction,  no  conclusive  model  could  be  developed  to  predict 
the  mode  levels  accurately.  The  analysis  was  successful  in 
establishing  criteria  for  shear  layer  stabilization  (a  reouire- 
ment  for  the  suppression  of  pressure  oscillations). 

To  evaluate  the  analytical  predictions  and  to  gain  practical 
information,  an  experimental  program  was  initiated.  One  important 
purpose  of  this  program  was  to  verify  the  resonant  frequency 
predictions  and  the  theoretically  derived  mode  shapes,  as  well 
as  to  provide  reliable  information  on  maximum  levels  and  mode- 
amplitude  distribution  in  the  cavity  for  a variety  of  cavity 
length-to-depth  ratios  and  freestream  Mach  numbers.  Another 
purpose  of  the  experiments  was  to  verify  the  criteria  for  shear 
layer  stabilization  and  to  find  and  optimize  geometric  changes 
of  the  basic  rectangular  cavity  configuration  for  maximum  oscilla- 
tion suppression. 
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In  order  to  attain  these  goals,  it  was  necessary  to 
conduct  the  experiments  In  a large  flow  facility  primarily 
to  eliminate  the  uncertainty  usually  associated  with  scaling. 

It  was  also  desirable  to  conduct  the  experiments  over  a large 
Mach  number  range,  to  duplicate  typical  flight  envelopes  of 
modern  aircraft.  Finally,  a test  setup  was  necessary  that 
allowed  quick  changes  in  cavity  geometry. 

The  experiments  were  conducted  in  the  NASA  Lewis  8 x 6 ft 
Transonic  Wind  Tunnel,  which  met  all  the  above  conditions. 

In  preparation  for  these  large-scale  experiments,  however, 
tests  were  conducted  in  a small-scale  wall-jet  flow  facility 
at  BBN . This  facility  allowed  a quick  survey  of  a large  number 
of  configurations,  so  that  some  optimization  of  geometric  changes 
could  be  made  prior  to  large-scale  testing.  These  small-scale 
tests,  however,  comprised  only  the  subsonic  Mach  number  ranee . 

In  a supplementary  effort,  therefore,  a water  table  flow  visualiza- 
tion apparatus  was  constructed  that  proved  to  be  helpful  In 
visualizing  the  unsteady  phenomena  inside  and  outside  the  cavity 
for  simulated  supersonic  external  flow.  The  water  table  tests 
also  served  to  check  stabilization  configurations. 

5.2  Experimental  Techniques:  Small-Scale  Studies 

5.2.1  Exploratory  Studies:  Wall-Jet  Flow  Facility 

Subsonic  flow  experiments  were  conducted  using  the  test 
setup  shown  in  Fig.  23.  High-pressure  air  is  fed  through  an 

absorptive  muffler  into  a converging  nozzle  with  a square 
opening  of  3 x 3 in2.  The  lower  nozzle  rim  is  flush  with  a 
plate,  which  contains  a rectangular  cavity  1 In.  wide  by  2.5  in. 
deep  by  8 in.  long.  When  Insert  blocks  are  used,  the  geometry 
can  be  readily  changed  to  form  cavities  of  length-to-depth 
ratios  from  0.5  to  20.  In  all  cases,  however,  the  width  Is 
constant  at  1 in.  The  flov;  emanating  from  the  nozzle  forms  a 
wall  jet,  and  the  flow  on  the  plate  surface  has  the  aerodynamic 
properties  of  any  surface  flow  with  a slow  boundary  layer 
growth  (see  Fig.  24).  Along  the  unconstrained  jet  boundaries, 
the  shear  layer  grows  in  accordance  with  the  free  mixing  process, 
surrounding  a potential  core  of  low  turbulence,  which  is 
reminiscent  of  quiescent  air  above  a surface  boundary  layer.  In 
all  experiments,  care  was  taken  that  the  cavity  mouth  area  was 
well  within  this  potential  core  regime,  so  that  the  ^ree  shear 
layer  of  the  wall  jet  would  not  affect  the  cavity  pressure  oscilla- 
tion process.  Since  only  cavity  internal  fluctuating  pressures 
were  measured,  rather  than  farfield  radiated  acoustic  pressures, 
there  was  no  danger  of  shear  layer  acoustic  refraction  to  affect 
the  measurement  data. 
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Prior  to  exiting  from  the  nozzle,  the  flow  passes  through 
a honeycomb  structure  to  extract  any  swirl  and  through  several 
mesh  wire  grids  to  minimize  the  turbulence  intensity  in  the 
potential  core  flow.  Sandpaper  roughness  across  the  exit  lip 
of  the  nozzle  provides  turbulent  boundary  layer  flow  to  interact 
with  the  cavity  mouth.  Since  it  had  been  established  in  flight 
tests  (Smith  et  al . 3 197*0  that  inlet  flow  in  aircraft  bays 
under  most  flight  conditions  is  turbulent,  most  of  the  experi- 
ments were  conducted  with  the  surface  roughness  at  the  nozzle 
exit.  Exit  flow  speed  could  be  varied  from  0 to  over  700  ft/sec. 
Plow  speed  is  monitored  through  a pitot-static  tube  in  the 
nozzle  exit  cross  section. 

Cavity  internal  fluctuating  pressures  were  measured  with 
BBN  type  376  1/4  in.-diam.  piezoelectric  pressure  sensors 
located  at  the  leading-edge  and  the  trailing-edge  bulkheads. 
Sensor  signals  were  usually  fed  into  a GR  type  1564A  1/10-octave 
band  analyzer,  and  spectra  plotted  out  on  a GR  type  1521  Oraphic 
Level  Recorder  in  real  time. 

One  important  issue  of  the  study  was  the  evaluation  of  the 
steady-state  (dc)  drag  induced  by  the  basic  cavity  configuration 
and  the  change  in  drag  due  to  geometric  cavity  modifications 
through  oscillation  suppression  devices.  To  measure  dc  drag, 
the  plate  containing  the  cavity  was  supported  by  low  friction 
ball  bearings  that,  in  turn,  ran  on  two  horizontal  steel  rods 
(see  Pig.  25).  The  drag  forces  (i.e.,  forces  in  the  mean 
flow  direction)  were  transferred  via  a 1/4-in.  diameter  steel 
ball  to  a (BBN-developed ) force  link,  which  was  attached  to  a 
heavy  steel  frame.  A bias  force  was  provided  through  a spring 
pressing  against  the  cavity  block  in  the  downstream  direction 
with  several  times  the  force  expected  from  the  cavity  drag. 

Force  transfer  through  a steel  ball  eliminated  the  possibility, 
of  any  force  or  moment  transfer,  other  than  in  the  direction 
of  the  sensitivity  axis  of  the  force  link,  which,  of  course, 
coincides  with  the  direction  of  flow-induced  drag.  The  system 
was  calibrated  by  applying  a known  force  in  the  flow  direction 
and  reading  the  output  voltage  of  the  force  gauge,  after 
cancelling  out  the  bias  force.  The  force  measurement  system 
was  perfectly  linear  in  the  range  of  measurements  taken. 
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Three-dimensional  flow  visualization  was  obtained  by  feeding 
smoke  into  the  cavity  near  the  floor  and  the  leading-edge  bulk- 
head. Short -duration  flash  photographs  allowed  visualization  of 
the  oscillating  shear  layer  in  a qualitative  manner. 

5.2.2  Exploratory  Studies:  Water  Table  Facility 

Since  the  small-scale  flow  facility  permitted  only  subsonic 
flow  experiments,  some  exploratory  studies  were  conducted  in  a 
water  table  facility  which  allows  supersonic  flow  simulation. 

A simple  water-  table  setup  permits  two-dimensional  simula- 
tion of  the  pressure  oscillation  phenomenon.  Since  pressure 
oscillations  in  shalluw  and  long  cavities  are  basically  two- 
dimensional  in  character,  a high  degree  of  similarity,  especially 
on  the  unsteady,  01  periodic,  part  of  the  phenomenon  can  be 
achieved . 

Figure  26  shows  a schematic  of  the  setup,  which  essentially 
involves  a slightly  inclined  transparent  glass  plate  supporting 
the  water  flow.  Using  the  optical  system  of  an  overhead  pro- 
jector, a direct  projection  on  a screen  is  possible;  the  surface 
waves  and  their  motion  are  shown  at  a high  contrast.  Optimum 
simulation  of  the  phenomenon  is  achieved  by  using  a water  depth 
of  about  1/4  in.;  blocks  defined  a cavity  of  2-in.  depth  (in 
the  plane  of  the  water),  and  variable  lengths  (2  in.  to  8 in.). 
Direct  observation  of  the  unsteady  phenomena  is  possible; 
however,  a motion  film  camera  with  a moderately  high  speed 
(e.g.,  54  frames  per  sec)  sufficiently  slows  the  phenomenon  so 
that  the  motion  of  the  shear  layer,  the  external  leading-edge 
and  trailing-edge  shock  wave  formation,  and  the  cavity  internal 
pressure  wave  (Heller  et  al . , 1973)  can  be  followed.  The 
analysis  developed  in  Sec.  4 of  this  report  is  essentially  based 
on  the  water  table  visualization  results. 

5.3  Experimental  Techniques:  Large-Scale  Studies 

5.3.1  Wind  Tunnel  Facility 

The  NASA  Lewis  8 by  6 ft  Wind  Tunnel  is  a continuous  closed- 
circuit  flow  facility.  A remotely  variable  nozzle  contour  allows 
operation  in  the  Mach  number  range  from  about  0.8  to  2.1.  Since 
the  facility  is  an  atmospheric  tunnel,  stagnation  and/or  dynamic 
pressure  of  the  free  stream  is  a function  of  the  operational  Mach 
number  and  cannot  be  changed  Independently  of  Mach  number.  A 
plot  of  stagnation  pressure,  dynamic  pressure,  and  unit  Reynolds 
number  as  function  of  tunnel  Mach  number  appears  in  Fig.  27. 

The  drop  in  static  pressure  for  increasing  Mach  number  defines  a 
simulated  altitude  as  shown  in  Fig.  28. 
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5.3.2  Cavity  Model 

The  cavity  model  is  shown  in  Figs.  29  to  31.  Cavity  width 
is  9 in. s the  length  is  36  in.,  and  the  floor  position  can  be 
remotely  controlled  between  a depth  of  16  in.  and  6 in.,  resulting 
in  length-to-depth  ratios  of  2.25  < L/D  < 6.0  (Fig.  29). 

Figure  30  shows  the  partially  disassembled  cavity  with  the  driving 
mechanism  for  the  floor  movement.  The  cavity  was  made  entirely  of 
1-in.  aluminum,  with  Teflon  strips  to  seal  the  floor  against  the 
cavity  walls.  Figure  30  also  shows  the  traverse  rod,  which 
contains  a BBN  type  376  sensor,  5 in.  below  the  cavity  surface 
plane.  The  rod  could  also  be  driven  remotely  to  survey  the 
cavity  internal  pressure  field.  Figure  31  shows  the  assembled 
cavity  with  two  extensions,  one  at  each  of  the  leading-  and  the 
trailing-edge  bulkheads.  These  extensions  measure  6 in.  long  by 
6 in.  deep  and  8,5  in.  wide,  which  increases  the  total  cavitv 
length  to  48  in.  However,  they  were  not  used  to  increase  the 
cavity  length  but  to  allow  changes  in  the  cavity  geometry  at 
leading-  and  trailing-edge  areas,  such  as  slants  or  roundings, 
or  to  accept  the  drive  mechanism  for  changing  the  angle  of  up- 
stream canard  spoilers  (see  Sec.  7-3).  Figure  32  also  shows 
two  remotely  drivable  rods  near  the  rear  end  of  the  cavitv; 
these  rods  support  leading-  or  trailing-edge  cowls  or  guidevanes 
and  vary  the  position  of  such  suppression  devices  in  the  vertical 
di recti  on . 

Hence,  in  any  given  test,  there  was  the  possibility  to  vary 
independently  the  position  of  the  floor,  the  traverse  rod,  a 
leading-  or  trailing-edge  oscillation  suppression  device,  and 
the  canard  spoilers.  Figure  33  shows  the  cavity  model  installed 
in  the  tunnel  test  section.  Since  the  test  section  walls  are 
perforated,  the  cavity  mouth  area  is  surrounded  by  a flat  smooth 
plate,  providing  a surface  for  natural  boundary  layer  growth.  The 
distance  from  the  plate  leading  edge  to  the  cavity  leading  edge  is 
76  in.  A closeup  of  the  installed  cavity  aopears  In  Fig.  34, 
which  shows  the  traverse  rod  and  two  laterally  displaced  oitot- 
static  tubes  in  the  floor. 

5 3.3  Instrumentation 

Both  mean-flow  and  unsteady-flow  quantities  were  measured. 
Mean-Flo u)  Quantities 

Boundary  Layer  Profile . A boundary  layer  rake  extending 
1.5  in.  above  the  surface  at  the  cavity  leading  edge,  but 
laterally  displaced,  provides  information  on  steady-state  boundary 
layer  characteristics,  specifically,  velocity-profile,  boundary 
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layer  thickness,  and  displacement  thickness.  Similar  informa- 
tion is  obtained  from  a rake  downstream  of  the  cavity  located 
on  the  cavity  centerline.  This  rake  extends  13  in.  above  the 
surface  plane.  Figure  35  provides  details  or.  rake  location. 


Recirculation  Velocity . Qualitative  information  on  recircula- 
tion velocities  is  obtained  through  two  pitot-static  tubes 
attached  to  the  cavity  floor  facing  in  the  downstream  direction. 


Static  Pressures . Several  static  pressure  taps  are  provided 
along  the  cavity  floor  and  the  forward  and  aft  bulkhead  (see 
Fig.  35). 


Cavity  Internal  Temperatures . Two  thermocouples  at  the 
locations  indicated  in  Fig.  35  provide  information  on  cavity 
internal  temperature  and,  hence,  on  cavity  internal  speed  of  sound, 


Unsteady  Quantities 


Fluctuating  Pressures . Fluctuating-pressure  levels  (sound 
pressure  levels)  are  measured  through  9 BBN  type  376  1/4-in. 
diameter  piezoelectric  pressure  sensors.  Sensor  locations  are 
shown  to  scale  in  Fig.  36.  There  is  one  sensor  each  in  the 
surface  plane  upstream  and  downstream  of  the  cavity  leading 
and  trailing  edge  on  the  centerline  (Sensors  1 and  8);  one 
each  in  the  forward  and  aft  bulkhead  (Sensors  2 and  7);  four 
in  the  floor,  two  of  these  in  a corner  position  (Sensors  3 
and  6);  and  two  near  the  center  of  the  floor  but  near  opposite 
side  walls  (Sensors  4 and  5).  One  sensor  is  imbedded  in  the 
traverse  rod  (Sensor  9). 


Vibration.  One  accelerometer  is  located  at  a middle  loca- 
tion on  the  floor,  outside  of  the  cavity,  to  monitor  vibration 
for  possible  interference  with  fluctuating-pressure  measure- 
ments . 


5.3.4  Data  Acquisition  and  Reduction 


Figure  37  shows  the  data  acquisition  system,  consisting  of 

9 BBN  type  376  pressure  transducers,  1 BBN  type  501  accelerometer, 

10  Preston  8300  XWB  Amplifiers  with  stepwise  adjustable  1 to 
1000  gain,  10  voltage  dividers,  and  1 Sangamo  4700  14-char.nel 
tape  recorder. 
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FIG.  35.  INSTRUMENTATION  FOR  STEADY-FLOW  QUANTITY  MEASUREMENTS 
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The  voltage  dividers  compensate  for  the  different  sensor 
sensitivities,  so  that  the  input  voltage  to  the  tape  recorder 
for  a given  pressure  level  was  the  same  for  each  pressure 
sensor  (provided  the  sensor  signal  is  put  through  the  same 
amplifier  gain).  All  recordings  are  made  at  30-ips  tape  speed; 
the  system  frequency  response  ranges  from  20  Hz  to  beyond 
MO, 000  Hz.  The  data  analysis  system,  which  is  shown  in  Fig.  38, 
consists  of  one  Sangamo  4700  14-channel  tape  recorder,  one  switch 
for  channel  selection,  and  one  Preston  Amplifier  Type  12048.  The 
signal  can  then  be  analyzed  either  in:  1/3-octave  bands  using 

a General  Radio  type  1925/1926  Real  Time  Analyzer  and  type  1522 
Grap!  ic  Level  Recorder;  or  in  10-Hz  constant  bandwidth  using  a 
General  Radio  type  1900-A  Wave  Analyzer  with  a type  1521-A 
Graphic  Level  Recorder.  A Hewlett  Packard  type  521C  frequency 
counter  is  used  to  determine  resonant  frequencies  accurately. 
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SECTION  6 

TEST  RESULTS:  BASIC  CAVITv 

6.1  Flow  Visualization 

In  this  section,  both  subsonic  and  supersonic  flow  visual- 
isation results  will  be  discussed.  Cavity  pressure  oscillation 
mechanisms  for  external  subsonic  and  supersonic  flow  are  quite 
similar;  hence,  the  subsonic  visualization  results  can  be 
interpreted  on  the  basis  of  the  sunersonic  visualization  results, 
for  which  much  more  detailed  information  was  obtained. 

Figure  39  represents  a complete  sequence  of  the  water  table 
simulation*  for  a typical  oscillation  cycle  for  external  super- 
sonic flow  of  about  M = 1.5.  The  motion  pictures  were  edited  to 
show  the  mass  injection  process,  which  starts  with  the  shedding 
of  a shear  layer  from  the  leading  edge  (frames  1 to  7).  The 
separated  shear  layer  approaches  the  trailing  edge  (frames  8 to 
12),  and  flows  over  the  trailing  edge  (frames  13  to  17).  There- 
after, the  shear  layer  "whips"  down  (frames  l8  to  24),  thus 
exposing  the  trailing  edge  to  the  free  stream,  which  causes  a 
"bow  wave,"  i.e.,  a shock  front,  to  occur  at  the  trailing  edge  in 
(frames  24  to  30).  Simultaneously,  the  downward  motion  of  the 
shear  layer  causes  a pressure  wave  (frame  24)  to  appear  in  the 
joint  region  of  the  trailing-edge  bulkhead  and  the  floor;  this 
pressure  wave  subsequently  moves  towards  the  leading  edge 
(frames  24  to  35),  straightening  out  along  its  path.  The  pres- 
sure wave  also  causes  the  shear  layer  above  to  bend  outward  into 
the  free  stream.  This  outward  bend  causes  a shock  front  to 
trail  along  with  the  pressure  wave  in  the  upstream  direction 
(frames  26  to  40).  Since  this  shock  front  travels  against  the 
free  stream  with  the  speed  of  the  cavity  internal  pressure  wave 
(roughly  with  Mach  1,  the  speed  of  sound  in  the  cavity),  the 
effective  speed  against  the  free  stream  is  higher  than  the  free- 
stream  speed.  Consequently,  the  trailing  shock  front  is  inclined 
more  than  would  correspond  to  the  freestream  Mach  angle.  Upon 
arrival  at  the  leading  edge  (frame  39),  the  cavity  internal  pres- 
sure wave  is  reflected  and  travels  downstream  in  the  cavity 
(frames  40  to  50).  Once  the  external  shock  wave  has  arrived 
at  the  leading  edge,  it  separates  from  its  originator,  i.e.,  the 
outward  bump  in  the  shear  layer,  and  trails  out  In  the  freestream 
medium.  Since  the  downstream-traveling  cavity  internal  pressure 


*In  this  sequence,  flow  is  from  right  to  left. 


79 


wave  propagates  with  the  cavity  internal  speed  of  sound  (roughly 
Mach  1),  the  pressure  wave  speed  is  now  subtracted  from  the 
external  flow  speed  and  the  "bump"  in  the  shear  layer  travels 
with  an  effectively  subsonic  speed,  which  causes  no  external 
shockwave.*  The  pressure  wave  travels  further  downstream  and 
dissipates  its  energy  along  its  cavity  internal  path.  Meanwhile, 
probably  around  frame  50  ± 5,  a new  downlash  of  the  shear  layer 
at  the  trailing  edge  occurs,  which  causes  a new  pressure  wave  to 
originate,  whereupon  the  process  repeats  itself. 

The  above  sequence  shows  clearly  the  generating  process 
for  a fundamental  mode.  However,  higher  order  modes  can  be 
simulated  as  shown  in  Fig.  40  (which,  again,  is  an  excerpt 
from  the  motion  picture).  Here,  the  cavity  trailing-edge  bulk- 
head was  moved  forward  and  backward  in  a piston-like  manner 
forcing  the  cavity  at  a higher  frequency  than  the  naturally!" 
occurring  fundamental  frequency. 

The  generation  mechanisms  for  supersonic  and  subsonic 
external  flow  are  believed  to  be  quite  similar.  This  belief 
is  substantiated  by  the  experimental  data  that  relate  resonant 
Strouhal  numbers  to  freestream  Mach  numbers  on  one  hand,  and 
resonant  mode  levels  and  freestream  Mach  numbers  on  the  other 
hand  (see  Secs.  6,6  and  6.8).  In  both  sets  of  data,  there  is 
a smooth  transition  from  the  subsonic  over  the  transonic  to 
the  supersonic  flow  regime.  The  main  difference  lies  in  the 
appearance  of  the  external  shear  layer.  In  the  supersonic  case, 
the  external  shear  layer  exhibits  a wavy,  though  continuous, 
shape,  resulting  in  an  up-and-down  motion  at  the  trailing  edge; 
as  described  above,  this  results  from  the  forcing  of  the  internal 
pressure  wave  motion.  In  the  subsonic  case,  the  external  shear 
layer  tends  to  roll  up  to  form  Individual  vortices.  In  the 
three-dimensional  visualization  photograph,  shown  in  Fig.  41, 
two  individual  vortices  convect  and  grow  In  a downstream  direction, 
This  instantaneous  picture  was  taken  in  the  small-scale  experi- 
mental flow  facility  (see  Sec.  5.2.1)  for  a Mach  number  of  about 
0.25.  These  pictures  correspond  Qualitatively  with  those  obtained 
by  Rossiter  (1966),  which  are  reproduced  here  as  Fig.  42. 


*The  external  flow  speed  in  this  simulation  was  about  Mach  1.5; 
hence,  the  effective  speed  of  the  shear  layer  bump  is  1.5  - 1 
= 0.5,  i.e.,  subsonic. 

*j* 

In  the  water  table  visualization  setup,  only  the  fundamental 
mode  can  sustain  Itself. 


As  has  been  emphasized  In  Sec.  3,  these  vortices  are  not 
the  cause  of  the  oscillation  as  assumed  by  Rossiter,  but  the 
external  manifestation  of  the  oscillation  process,  which  is  a 
result  of  cavity  internal  pressure  wave  motions. 

6.2  Boundary  Layer  Characteristics 

As  shown  in  Pig.  27,  the  unit  Reynolds  number  in  the  NASA 
Lewis  Research  Center  8 x 6 ft  tunnel  at  Mach  numbers  above 
0.8  is  between  4,2  and  4.8  x 106 , which  results  in  Reynolds 
numbers  at  the  cavity  mouth  of  26  to  30  million.  These  Reynolds 
numbers  are  high  enough  to  result  in  fully  developed  turbulent 
flow  at  the  cavity  mouth. 

Boundary  layer  profiles  were  measured  at  the  sideline  ex- 
tension of  the  cavity  mouth,  and  the  results  are  shown  in 
Fig-  43.  The  profiles  in  the  Mach  number  range  or  0.8  < < 1.5 

coincide,  showing,  the  shape  of  a fully  developed  turbulent 
boundary  layer.  The  exponent  N in  the  profile  equation,  U/U 
= (y/6)l/N,  corresponds  to  7-3.  The  Mach  2 velocity  profile” 
deviates  in  shape  from  the  lower  Mach  number  profiles,  with  an 
exponent  N of  about  10. 5- 

The  boundary  layer  displacement  thickness  is  about  0.165 
+ 0.006  In.  for  0.8  < < 1.5;  thus,  it  is  a factor  of  35  to  95 

times  smaller  than  the  cavity  depth.  Displacement  thickness  at 
Ma  = 2 is  about  0.14  in.  (Fig.  44).  Correspondingly,  the 
boundary  layer  thickness  is  about  1.0  in.  ± 0.13  for  0.8  < M^ 

< 1.5  and  about  0.78  in.  for  M =2.  Hence,  for  the  shallow” 
cavity  experiments  (L/D  = 6),  ?he  ratio  of  cavitv  depth  to 
boundary  layer  thickness  (D/6)  is  6 to  7;  for  the  deep  cavitv 
experiments  (L/D  = 2.25),  D/6  is  16  to  20.  Thus,  the  boundary 
layer  dimensions  are  many  times  smaller  than  the  cavity  depth. 

6.3  -Mean  Pressure 

Static  pressures  in  the  cavity  were  measured  at  six  locations, 
four  on  the  floor,  and  one  each  on  the  forward  and  the  aft  bulk- 
head. Figure  45  presents  the  static  pressure,  on  the  floor 

at  0.1,  0.33,  0.66,  and  0.9  of  the  cavity  lengt h ,”normali zed 
with  the  freest-ream  total  pressure,  PD  . The  information  is 
riven  for  freestream  Mach  numbers  ranging  from  0.78  to  2.0, 
and  broken  down  for  three  L/D  ratios.  There  is  no  obvious 
dependence  on  the  L/D  ratio,  but  the  pressure  ratio  at  each 
Mach  number  rises  towards  the  trailing  edge  (x/L  = 1).  further- 
more, the  lower  the  Mach  number,  the  higher  is  the  pressure 
ratio . 
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FIG.  44.  MACH  NUMBER  DEPENDENCE  OF  BOUNDARY  LAYER 
DISPLACEMENT  THICKNESS  A1  CAVITY  LEADING 
EDGE  . 
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Figure  46  i' i1' ' *v! i i Un  i::'  r i-c a v 1 t y *•]  tx-r  si  n 1. 1 c pressure# 
.■■•'iju.'i 'd  by  1 In'  1. r i -an  f ;.i 1. 1 e nror-nure , and  referenced  to 

liic  mean -ea  v ] V y static  piViV-uri' , a:?  function  of  Mn-'h  number' 
per  three  l./D  rat  1 oi! . The  plot-  coni  ains  data  from  ov-hex'  invosfi- 
■ at  or:'-.  At  sub. font  c Mae);  numbor:; , there  1 s no  rone  I</D 
dependence . At  unperson  1 c Mach  numbers,  the  dal  a fan  out,  so 
vital  lil -“her  normal  1 ami  preonupe a are  observed  for  lower  L/D 
ratios,  and  vice  verso, 

6.4  Cavity  Temperature 

Temj  orature;;  wore  measured  at  two  locations.  on  the  floor'. 

IV it  h are  pr-us,  ntod  ar,  a recovery  factor  r 5 (,:<  -TO))/(",0-':-,«M)  v" 
roam  Mach  number;*  In  Fir.  47.  The  recovery  factor  seems 
to  reach  a relative  minimum  at  sonic  and  slirhtlv  sunersonic 
Mach  numborr  , i ncrea r.  1 up  both  toward  lower  and  hirher  Mach 
numbers . A systematic  dependence  on  lem-th-to-denth  ratio  if 
evident,  no  that  lowpr  L/D  rat ios  result  in  lower  recovery 
factors.  The  data,  which  were  obtained  in  the  NASA  Lewis 
Research  Center  S x 6 ft  funnel,  however,  am  systematically 
lowei'  1 han  those  obtained  In  a previous  study  (Heller,  1^70) 
in  ihe  MIT  Naval  Supersonic  Tunnel. 


Heller  at.  a 7 . (19  70)  have  arpued  that  the  speed  of  sound 
in  the  cavity  approaches  the  starvation  speed  of  sound  in  the 
Ills  was  a key  postulate  in  riodifyinr  Rossitor  ’ - 


1 1'iio  dti*earn . 


(1.966)  equation  that  relates  resonant  frequencies  In  the  cavity 
t.o  freest  roam  Mach  numbers.  The  measured  data  in  '-‘if.  4f  support 
this  earlier  conclusion;  recovery  factors  are  closer  to  unity 
than  to  sere,  as  had  been  assumed  by  Rossiter. 


6.5  Fluctuating-Pressure  Lata 


6.5.1  Spectra 


Forty  runs  were  conducted;  each  run  has  a different  Mach 
number  and  leaf th-to-der th  ratio  of  the  cavity.  Since  the  cavity 
contains  nine  ( f luetuatinp ) pressure  sensors  and  one  accelero- 
meter, close  to  SuO  individual  data  points  were  obtained. 

Tests  were  conducted  at  the  foliowinr  r,ach  numbers:  0.R, 

0.9,  0.95,  1.0,  1.05,  1.1,  l.P,  1.5,  and  ?.C,  ^or  most  of  these 
Mach  numbers,  data  were  recorded  at  three  L/D  ratios,  i.e.,  P.3, 


KAvei'aped  over  the  four  floor  positions. 
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TABLE  I.  CONVERSION  FACTORS 


M 

CO 

Ordinate 
Correcti on* 
(dB) 

Abscissa 
Correcti on+ 

0.8 

184 

0.00341 

0.9 

185 

0.00307 

1.0 

186 

0.00281 

1.1 

187 

0.00262 

1.2 

187.5 

0.00239 

1.5 

188.5 

0.00202 

2.0 

189.5 

0.00162 

Of  course,  the  frequency  scale  Is  only  valid  for  a cavity 
of  3-ft  length.  To  nondimensionalize  the  frequency  on  a Strouhal 
number  basis,  we  must  multiply  the  frequency  abscissas  by  the 
appropriate  numbers  in  Table  I. 

As  further  illustration  of  the  spectral  energy  distribution 
in  and  near  the  cavity,  Figs.  49  and  50  show  1/3-octave  band 
spectra  obtained  for  an  L/D  = 2.3  cavity  at  various  locations  at 
Mot  = 0.8  and  M„o  = 1.5,  respectively.  Since  sensor  signals  on  a 
cavity  leading-  or  trailing-edge  wall,  and  in  corners,  are  more 
easily  compared  than  signals  from  sensors  halfway  in  the  cavity, 
where  levels  are  affected  by  location, only  signals  from  sensors 
1,  2,  3 3 6,  7,  and  8 (see  Fig.  36)  are  shown. 


* To  convert  ordinate  levels  in  Fig.  49  to  read  20  log  (prin,,/qoo) , 
subtract  the  numbers  in  the  column  from  the  ordinate  levels. 

+ To  convert  abscissa  frequencies  in  Fig.  49  to  Strouhal 
frequencies,  S = fL/U^,  multiply  the  frequencies  by  the 
numbers  in  the  column. 
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Generally , we  can  conclude  that  spectral  energy  distribution 
and  levels  are  quite  uniform  inside  the  cavity  near  the  leading- 
edge  bulkhead.  At  subsonic  speeds,  acoustic  energy  is  radiated 
upstream,  so  that  a generically  similar  spectrum,  although  of 
lower  intensity,  is  observed  a short  distance  upstream  of  the 
leading,  edge.  Acoustic  characteristics  in  the  trai ling-edge 
region  are  also  quite  similar;  however,  the  discrete  tone  levels 
are  higher  than  at  the  leading-edge  bulkhead,  and  broadband  noise 
has  increased  significantly.  Acoustic  signatures,  a small 
distance  downstream  of  the  cavity  trailing  edge,  correspond  to 
those  inside  the  cavity,  but  the  levels  are  slightly  reduced. 

The  data  presented  in  Figs.  ^9  and  50  can  also  be  nondimen- 
sionalized  using  Table  I. 

A more  detailed  picture  about  the  energy  distribution  be- 
tween discrete  tones  and  broadband  energy  is  obtained  through 
narrowband  analysis.  In  Fig.  51 , 10-Hz  constant  bandwidth 
spectra  spanning  the  0 to  1000-Hz  frequency  range  are  presented 
for  three  L/D  ratios  at  the  leading-  and  trailing-edge  regions. 
The  Mach  number  in  this  test  was  unity.  These  spectra  show 
that  near  the  leading  edge  for  deep  cavities  (L/D  = 2.3)  almost 
all  energy  is  in  discrete  frequency  bands;  in  fact,  the  broad- 
band energy  (in  this  10-Hz  bandwidth  representation)  is  about 
hO  dB  down  from  the  dominant  mode-2.  Tone  levels  near  the 
trailing  edge  are  comparable  to,  or  exceed,  those  at  the  leading 
edge,  however,  the  broadband  energy  has  significantly  increased; 
it  is  only  25  dB  below  the  mode-2  peak. 

For  the  shallower  cavity  (L/D  = 5.5),  the  energy  is  somewhat 
more  evenly  distributed  between  modes;  broadband  energy  at  the 
leading  edge  is  only  25  dB  below  the  mode  - 2 peak.  At  the 
trailing  edge,  the  mode-2  peak  exceeds  the  broadband  energy  level 
by  about  18  dB. 

6.5.2  Resonant  Frequencies 

As  has  been  discussed  in  previous  sections,  flow-exposed 
cavities  resonate  at  discrete  frequencies,  which  order  them- 
selves along  certain  modes.  In  order  to  determine  resonant 
frequencies  accurately,  the  frequencies  were  counted  by 
electronic  means. 

Figure  52  shows  nondimensional  frequencies,  S = fmL/Uoo. 
as  function  of  freestream  Mach  number  M^.  here,  fij  is  a 
resonant  frequency  with  m = 1,  2,  3 •••  denoting  the  mode  under 
consideration,  L is  the  cavity  length , and  U*  is  the  freestream 
velocity.  Four  sets  of  data  are  contained:  (1)  test  results 

for  2.3  < L/D  <5.5  obtained  at  the  NASA  Lewis  Research  Center 
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FIG.  52.  STROUHAL  FREQUENCIES  OF  CAVITY  MODES  AS  FUNCTION 
OF  MACH  NUMBER. 
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8 x 6 ft  wind  tunnel;  (2)  test  results  from  a previous  research 
effort  obtained  for  4 < L/D  < 7 in  the  MIT  Naval  Supersonic  Wind 
Tunnel  (Heller  et  al.t  1970);  (3)  Air  Force  flight  test  data  for 
L/D  = iJ.O  (Smith  et  al.a  197*0;  and  (O  small-scale  test  data 
obtained  in  the  wall-jet  facility  for  2.0  < L/D  < 5.0.  This 
plot  spans  the  Mach  number  range  from  C.22  to  3.0. 

It  is  evident  that  all  data  points  order  themselves  closely 
among  the  various  mode  curves  that  were  previously  (Heller  et  al., 
1970)  expressed  analytically  as 


if- 

n 

K:  s' 


n — m 

m = ~TT 


M/t  (l+-^j—  • M2)1/2]  + l/kv 


(6.1) 


where  a is  an  empirical 
L/D  = 4,  y is  the  ratio 
empirical  constant  tied 
the  cavity  shear  layer 
effect  of  the  L/D  ratio 
discernable  in  the  2 < 
scatter.  Qualitatively 
ratios  tend  to  increase 


constant,  which  Is  close  to  0.25  for 
of  specific  heats,  and  k is  another 
to  a disturbance  convection  speed  in 
and  found  to  equal  0.57.  No  systematic 
on  nondimensional  frequencies  is  readily 
L/D  < 7 range , because  of  experimental 
, then,  we  may  state  that  the  lower  L/D 
the  nondimensional  frequency. 


At  Mach  numbers  below  0.5,  the  data  points  rise  faster 
(towards  M = 0)  than  predicted  by  Eq . 6.1,  which  restricts 
the  validity  of  the  equation  to  M > 0.5. 

Resonant  frequencies  up  to  the  fifth  mode  and,  possibly, 
beyond  were  observed;  however,  they  had  dramatically  different 
modal  Intensities,  depending  on  both  Mach  number  and  length-to- 
depth  ratios.  Hence,  while  a resonance  can  only  occur  at 
frequencies  determined  by  Eq.  6.1,  this  representation  of  mode 
frequencies  does  not  determine  whether  a resonance  will,  in 
fact,  occur. 

6.5.3  Effect  of  L/D  Ratio  on  Frequencies 

Since  the  general  scatter  of  the  large  body  of  experimental 
data  does  r.ot  allow  an  accurate  evaluation  of  the  mode  frequency 
dependence  on  the  L/D  ratio,  an  experiment  on  the  small-scale 
wall-jet  facility  was  initiated.  In  this  experiment,  the  L/D 
ratio  was  systematically  varied  from  0.3  to  6.67,  thus  encompassing 
the  range  from  very  deep  to  very  shallow  cavities.  The  external 


flow  speed  was  kept  constant  at  2 ll 5 ft/sec.  Figure  53  presents 
the  results  In  terms  of  Strouhal  numbers  vs  L/D  ratio. 

Deep  cavity  nondlmenslonal  mode  frequencies  Increase  rapidly 
with  L/D  up  to  L/D  * 1.  From  dimensional  considerations,  it  is 
clear  that  for  deep  cavities  the  dependence  of  S on  mode  number, 
Mach  number,  and  L/D  must  be  S 1 (m/M)  • (L/D);  i.e.,  for  constant 
m and  M,  the  Strouhal  number  is  directly  proportional  to  L/D.  On 
the  other  hand,  as  shown  by  Heller  et  al.  (1970),  shallow  cavity 
nondlmenslonal  frequencies  are  nearly  independent  of  L/D. 

This  behavior  is  quite  clear  in  Fig.  53.  For  shallow 
cavities  (L/D  > 2),  the  nondlmenslonal  resonant  frequencies 
decrease  as  the  cavity  becomes  shallower.  At  these  relatively 
low  Mach  numbers  (M  = 0.22),  the  increase  of  the  L/D  ratio  from 
2 to  4 results  in  a 20%  decrease  in  Strouhal  number.  However, 
there  is  reason  to  believe  that  at  higher  external  flow  speed 
the  effect  Is  less  dramatic,  and  frequency  variations  for  the 
important  first  three  modes  at  transonic  and  supersonic  flow 
speeds  are  probably  less  than  5%. 

6.5.4  Levels 

The  question  of  oscillation  onset  is  best  answered  by 
Figs.  54  and  55,  where  nondimensional  resonant  levels  are 
presented.  These  figures  pertain  to  levels  that  were  measured 
in  the  leading-edge  bulkhead  region  and  the  trailing-edge  bulk- 
head region,  respectively.  Data  are  presented  for  modes  1,  2, 
and  3 and  for  L/D  ratios  of  2.3,  4.0,  and  5.1.  Experimental 
results  from  the  large-scale  NASA  Lewis  Research  Center  8 x 6 ft 
tunnel  tests,  the  Air  Force  flight  tests,  the  previous  MT171  Naval 
Supersonic  Wind  Tunnel  test,  and  the  small-scale  wall-jet  tests 
are  included.  It  should  be  emphasized  that  this  representation 
is  not  a plot  of  the  overall  cavity  noise  level  vs  Mach  number, 
but  a plot  of  the  level  of  each  individual  mode  vs  Mach  number. 

An  attempt  was  made  to  connect  the  data  points,  which.  In  the 
leading-edge  bulkhead  region,  order  themselves  along  single 
curves  (see  Fig.  54).  The  uncertainty  about  modal  levels  is 
much  greater  in  the  trailing-edge  bulkhead  data  (see  Fig.  55). 
Here,  the  discrepancy  between  flight  test  and  wind  tunnel  test 
is  quite  pronounced. 

For  cavities  of  L/D  > 2,  It  can  be  concluded  that  the 
oscillation  onset  is  gradual,  as  exemplified  in  Pig.  56. 

For  the  specific  case  of  a mode-2  resonance  and  an  L/D  of 
four,  this  figure  shows  the  rapid  but  smooth  level  increase 
from  M=0.2toM=l,  reaching  a relative  maximum  at  M = 1.2 
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and  decreasing  thereafter  more  gradual  L.v . Resonant  mode  behavior 
for  other  conditions,  an  documented  In  Fir,.  54  is  generically 
quite  similar. 

In  the  follow'  ni*  ".erics  of  flg.urer.  (Pig.  57),  the  data 
of  Figs . 54  and  55  are  presented  again  to  allow  direct 
comparison  of  mode  levels  for  various  L/D  ratios, 

6.5.5  Mode  Shapes 

Previous  investigations  (Smith  et  al. , 1974;  Heller  et  al . t 
1970)  demonstrated  that  the  longitudinal  amplitude  distribution 
of  resonant  frequencies  shows  great  similarity  to  the  resonant 
behavior  of  a one-dimensional  oscillator  resonating  in  a length- 
wise direction.  However,  the  resonant  frequencies  do  not 
directly  correspond  to  ideal  "closed-box"  longitudinal  frequencies, 
other  than  at  high  Mach  numbers,  where  the  external  flow  re- 
presents a rather  stiff  boundary. 

Tne  setting  up  of  an  ideal  standing  wave  pattern  between  the 
fore  and  the  aft.  bulkhead  wall  of  a resonating  cavity  lequlres  { 

these  walls  to  be  acoustically  hard.  However,  this  is  not  the  j: 

case,  at  least  as  far  as  the  aft  bulkhead.  The  high  amplitude  ! 

oscillatory  motion  of  the  shear  layer  at  and  near  the  trailing 
edge  affects  the  impedance  of  the  cavity  aft  region.  Thus,  , 

while  standing  wave  patterns  are  definitely  observed,  they  do 
not  coincide  with  a classical  closed-space  standing  wave  pattern. 

Using  the  traversing  rod,  the  mode  shapes  for  the  first 
three  modes  were  traced  across  the  cavity  spar,  from  0.3  < x/L 
' 0.9,  where  x is  the  running  length  coordinate.  These  traces 
were  obtained  only  at  = 0.8  and  M =1.0,  and  for  2.3  < I,/D 
< 5.5*  A typical  trace°°appears  in  Fig . 58.  In  all  cases,  two 
traces,  i.e,,  one  upstream  and  one  downstream,  were  obtained  at 
each  experimental  condition.  The  bandwidth  of  the  trace  is  50  Hs. 

All  mode  shape  traces  are  shown  in  Figs.  59  and  60. 

In  some  of  the  trace  presentations,  data  points  from  the  Air  j 

Force  flight  test  (.Smith  et  al.  , 1974)  are  included.  These  aata 
points  were  obtained  at  Mm  = 0.8  for  altitudes  of  3000  and  20,000 
ft;  thus,  they  are  considered  to  be  comparable  to  the  10,C00-ft  • 

altitude  simulation  in  the  NASA  Lewis  facility  at  Mm  = 0.8.  These 
flight  data  points  were  shifted  so  that  they  coincide  with  tne  < 

mode  traces  from  the  wind  tunnel  tests.  The  points  fall  within  j 

a i4-dD  range  of  the  traces.  It  should  be  realised,  however,  - 

that  the  mode  traces  were  obtained  at  a fixed  distance  from  the 
cavity  mouth  plane,  whereas  tne  flight  test  data  points  were 
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FIG.  57b.  MACH  NUMBER  DEPENDENCE  OF 
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DOWNSTREAM  TRACE — H r- UPSTREAM  TRACE 


MODE  3 


COMPARISON  0!  MODE  SHAPES  FOR  CAVITIES  WITH  VARIOUS  LENGTH-TO-DEPTH 
RATIOS  M = 0.8  (continuous  lines  - present  investigation,  data 
points  -“Air  Force  flight  test  results). 


measured  on  the  floor.  It  is  likely  that  the  acoustic  field 
in  the  cavity  changes  in  a vertical  direction,  and  the  data  at 
various  vertical  positions  cannot  be  compared  directly. 

Inspection  of  the  mode  shape  traces  for  Mro  = 1 allows 
several  general  conclusions: 

• While  all  modes  exhibit  a fairly  well-defined  pressure 
minimum,  the  longitudinal  location  of  the  minimum  does 
not  correspond  to  the  classical  closed-box  modal  patterns. 

• The  pressure  minimum  for  mode-1  is  displaced  from  the 
center  of  the  cavity  towards  the  leading-edge  bulkhead 
for  the  shallow  cavities.  Specifically,  the  minima 
appear  at  the  following  x/L  locations. 


MODE  1 


L/D 

x/L 

2.3 

0.50 

4.0 

0.45 

5.1 

0.37 

Thus,  the  shallower  a cavity,  the  more  the  system  deviates 
from  a classical  one-dimensional  (hard-walled)  oscillator. 
In  all  cases,  the  trailing-edge  maximum  occurs  a small 
distance  upstream  of  the  trailing-edge  bulkhead.  These 
conclusions  are  valid  at  Mach  numbers  of  0.8  and  1.0. 

• The  downstream  pressure  minimum  for  mode-2  is  displa.^d 
In  the  upstream  direction,  while  the  upstream  pressure 
minimum  appears  roughly  at  the  location  for  a classical 
one-dimensional  oscillator  (i.e.,  at  x/L  = 0.25).  How- 
ever, with  increasing  Mach  number  (going  from  0.8  to  1.0), 
the  downstream  minimum  approaches  the  hard-wall  oscillator 
value  of  x/L  = 0.75.  Specifically,  the  minima  appear  at 
the  following  x/L  locations: 
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Furthermore? , the  trai  1 ing-edge  bulkhead  pressure  maxima 
are  displaced  in  t.he  downstream  direction. 
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• From  the  displacements  of  the  pressure  maxima  and  minima, 
we  can  derive  that  the  cavity  acts  acoustically  as  an 
effectively  shorter  cavity.  Qualitatively,  the  lower  the 
Mach  number  of  the?  external  flow,  the  more  the  cavity 
behavior  deviates  from  that  of  an  ideal  one-dimensional 
oscillator . 

• From  the  standing  wave  patterns,  we  can  define  an  effective 
reflection  coefficient  of  the  rear  wall.  This  reflection 
coefficient  does  not.  indicate  in  a physical  sense  a true 
absorptive  property  of  the  rear  bulkhead,  but,  rather,  it 
is  a measure  of  the  relative  magnitude  of  the  pressure 
maximum/minimum  ratio  within  each  modal  pattern. 

The  following  "Reflection  Coefficients"  were  determined, 
averaged  over  all  L/D  ratios  (<1.3  < L/I")  < 5«1): 


M 

Mode  1 

Mode  2 

Mode  3 

0.8 

0.40 

0.55 

0.50 

1.0 

o 

ro 

CO 

0.?6 

0.40 

1 


As  the  reflection  coefficient  rises,  the  trailing-edge  walls 
act  harder  acoustically.  However,  the  qualitative  nature  of 
this  conclusion  should  be  emphasized,  and  no  extrapolation  of 
the  above  numbers  is  attempted  at  lower  and  higher  Mach  numbers. 

6.5.6  Broadband  Noise 


Previously,  it  was  discussed  that  flow-induced  cavity  noise 
is  composed  of  both  discrete  and  random  (i.e..  broadband) 
components.  It  was  noted  that  broadband  noise  is  substantially 
higher  in  the  aft  cavity  region  and  relatively  weak  in  the 
leading-edge  area.  Thus,  an  attempt  was  made  to  separate  t.he 
broadband  from  the  discrete  sound,  and  to  find  some  commonality 
for  all  broadband  spectra  in  the  trailing-edge  region. 


There  is  some  speculation  in  deriving  the  broadband  portion 
from  a spectrum  that  contains  both  components.  Figure  6l 
presents  a subsonic  (M  = 0.8)  and  a supersonic  (M  = 1.5)  cavity 
pressure  spectrum,  as  measured  through  Sensor  7 (trailing-edge  ’ J 
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bulkhead).  These  data  pertain  to  the  deep  cavity  (L/D  = 2.3). 
where  It  was  previously  determined  that  only  modes  1,  2,  and  4 
appear  strongly.  If  we  disregard  those  1/3-octave  bands,  which 
contain  the  appropriate  resonance  frequencies,  we  can  easily 
sketch  a haystack-shaped  "lower  bound"  spectrum  that  is  attri- 
butable to  the  nondiscrete  components.  All  tralling-edge  spectra, 
encompassing  the  range  2.3  < L/D  < 5.5  and  0.8  < M < 2.0,  were 
obtained  in  this  fashion.  It  was  found  that  all  broadband 
spectra  [expressed  in  terms  of  (20  log  Prms/o°°)  vs  Strouhal 
number  (S  = fL/Um)]  fall  within  a range  of  about  ±2.5  dB,  and 
that  there  is  not  any  discernible  dependence  on  Mach  number  in 
this  Mach  number  range  (Fig.  63). 

Broadband  noise  levels  in  the  leading-edge  region  were  found 
to  be  at  least  10  dB  lower  than  at  the  trailing  edge.  A similar 
result  was  reported  by  Smith  et  at.  (1974)  in  their  flight  test 
program. 

6.5.7  Consecutive  Cavities 

A series  of  tests  was  conducted  on  two  consecutive  cavities 
as  shown  in  Fig.  63.  The  cavities  were  formed  by  inserting  a 
1-in.  thick  center  plate  with  the  cavity  floor  in  the  topmost 
position,  i.e.,  6 in.  down  from  the  cavity  mouth.  Hence,  two 
cavities,  which  were  17.5  in.  long,  6 in.  deep, and  9 in.  wide, 
were  formed.  Thus,  the  L/D  ratio  of  each  is  2.9.  Sensors  3 
and  4 (see  Fig.  38)  are  now  located  near  the  leading-  and 
the  traillng-edge  bulkhead,  respectively,  of  the  upstream  cavity; 
Sensors  5 and  6 are  at  corresponding  locations  in  the  downstream 
cavity . 

Experiments  were  conducted  at  Mach  numbers  of  0.8,  0.9,  0.95, 
1.0,  1.05,  1.1,  1.2,  1.5,  and  2.0.  Figure  64  presents  an  ex- 
ample of  1/3-octave  band  spectra  at  two  corresponding  measure- 
ment locations  in  both  cavities  for  subsonic  (M  = 0,8)  and 
supersonic  (M  = 1.5)  Mach  numbers.  It  is  evident  that  sceotra 
at  corresponding  locations  are  similar,  with  the  same  relative 
distribution  of  discrete  and  broadband  energy. 


A detailed  plot  of  the  change  of  nondimensional  mode  level 
with  freestream  Mach  numbers  in  the  range  0.8  < Mro  < 2.0  for 
the  first  three  modes  appears  in  Fig.  65.  This  figure  shows 
the  similar  behavior  of  the  upstream  and  downstream  cavities. 
Levels  in  the  leading-edge  area  are  3 to  9 33  lower  than  at 
the  trailing  edge. 
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f.  The  results  shown  in  Fig.  6 5 can  be  compared  to  Figs.  54 

1 and  55,  where  the  Mach  number  dependence  of  niode  levels  is 

shown  for  the  single  cavity.  For  the  double  cavity,  the  normalized 
levels  seem  to  be  somewhat  lower  than  for  the  single  cavity  case. 

The  significant  difference  between  the  single  and  the  double 
cavity  experiments  is  the  relative  width.  For  the  single  cavity, 
length-to-width  ratio  was  4;  for  the  double  cavity,  it  was  about 
1.9-  Although  the  cavity  oscillation  phenomenon  essentially 
has  a two-dimensional  nature  (as  documented  in  the  independence 
of  resonant  frequencies  from  the  length-to-width  ratio),  the 
levels  are  apparently  affected  by  the  width,  since  the  proximity 
of  the  side  walls  is  likely  to  influence  the  characteristics  of 
the  cavity  internal  recirculating  flow.  The  water  table  experi- 
ments also  indicated  that  consecutive  cavities  oscillate  strongly 
in  phase,  largely  independent  of  the  width  of  the  dividing  wall. 

We  are  uncertain  whether  this  finding  pertains  also  to  the 
conditions  in  air.  However,  the  phase  characteristics  of  con- 
secutive cavities  are  important  in  the  suppression  of  oscillations. 
Figure  65  presents  motion  picture  frames  from  the  water  table 
experiments  for  consecutive  cavities,  indicating  the  in-phase 
oscillation  process. 

6.6  Aerodynamic  Mean  Drag 

Flow-induced  cavity  drag  was  measured  with  the  apparatus 
described  in  Sec.  5.2.  Only  data  for  subsonic  flow  speeds 
were  obtained  because  of  the  limitations  of  the  small-scale 
flow  facility.  The  reason  for  these  measurements  was  twofold: 

First,  no  reliable  data  exist  for  the  drag  of  long  and  shallow 
cavities  in  flow  for  an  approaching  thin  boundary  layer;  second, 
some  information  on  the  effect  on  drag  of  the  oscillation 
suppression  devices  was  desired  to  evaluate  the  drag  penalty 
to  be  expected  under  flight  conditions. 

In  a typical  test  run,  the  flow  speed  over  the  plate 
containing  the  cavity  was  increased  from  zero  to  about  600  ft/sec; 
thereafter,  it  was  reduced  back  to  zero.  In  each  test,  a continuous 
plot  of  ihe  dc  drag  force  vs  the  wall-jet  total  pressure  was 
obtained.  A typical  plot  appears  In  Fig.  67.  Since  both  the 
drag  force,  D,  and  the  total  pressure,  P0,  in  the  flow  are  roughly 
proportional  to  (up  to  about  600  ft/sec),  a straight-line 

plot  is  obtained  in  a D vs  P0  linear-linear  representation,  where  the 
slope  is  a measure  of  the  drag  force  coefficient,  Cp.  Any  change 
in  slope  indicates  a change  in  the  value  of  Cp. 
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TYPICAL  TRACE  OF  FLOW-INDUCED  CAVITY  DRAG 
LAMINAR  INFLOW. 


r-wsMasMM^p^T!^  «w-^ 


The  measured  drag  is  composed  of  that  which  is  induced  by 
the  cavity,  and  th.'t  which  is  due  to  the  surrounding  surface 
being  under  the  highly  turbulent  mixing  region  of  the  wall  jet. 
Since  the  cavity  drag  force  and  the  plate  drag  force  add  to 
each  other,  the  cavi uy  contribution  can  be  obtained  by  sub- 
tracting the  plate  contribution. 

Cavities  of  various  length-to-depth  ratios  were  tested 
under  laminar  inflow  conditions.  Some  test  data  were  obtained 
for  turbulent  inflow,  which  was  achieved  by  tripping  the  nozzle 
boundary  layer  at  the  nozzle  exit. 

A plot  showing  the  drag  coefficient  vs  length-to-depth 
ratio  is  shown  in  Fig.  68.  The  drag  coefficient  of  the 
cavity  is  defined  as 


^total  ~ Opiate 


Acavi.ty  * 


Here,  Dtotai  is  the  measured  drag,  D -,ate  is  that  portion  of 
the  drag  attributable  to  the  surrounding  plate,  Acav:j_ty  is  the 
mouth  area  of  the  cavity,  and  q„  = l/2oU„  is  the  freestream 
dynamic  pressure.  The  data  in  Fig.  68  are  valid  in  a speed 
range  from  0 < Um  < 600  ft/sec.  The  determining  factor  for  the 
drag  coefficient  seems  to  be  the  cavity  L/D  ratio,  which  causes 
drag  maximum  at  an  L/D  value  of  10.  Fluid  dynamic  considera- 
tions indicate  that  the  free  shear  layer  over  short  deep  cavities 
span  the  gap,  resulting  in  a low-energy  recirculating  flow 
pattern.  On  the  other  hand,  long  shallow  cavities  allow  the 
free  shear  flow  to  attach  to  the  floor,  causing  onlv  two  weak 
recirculation  regions  near  the  leading-edge  and  trailing-edge 
bulkheads.  Presumably,  cavit ios,  with  a certain  critical  length- 
to-depth  ratio  (evidently  near  10),  cause  two  strong  recirculating 
regions  near  the  leading-  and  the  trailing-edge  bulkhead  with  a 
relatively  brief  attachment  region  in  between.  This  flow  pattern 
causes  high  drag.  Qualitative  flow  patterns  appear  in  Fig.  68. 

Figure  68  is  qualitative  in  nature;  the  effect  of  cavity 
width  is  neglected,  and  the  ratio  of  boundary  layer  to  depth 
has  not  been  considered.  Although  the  cavity  mouth  area  (length 
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times  width)  was  found  to  be  a rood  normalising  parameter,*  the 
width-to-length  ratio  may  be  ,1ust  as  important  as  the  length- 
to-depth  ratio  in  determini  nr.  the  actual  cavitv  drar.  however, 
Fir,.  68  could  be  used  to  obtain  rough  estimates  of  the  drag 
of  long+  shallow  cavities  at  fairly  low  subsonic  speeds  for 
approaching  thin  laminar  boundary  layers. 

dome  data  points  (open  marks)  in  Fir,.  68  pertain  to 
turbulent  inflow.  They  indicate  that  turbulent  inflow  causes 
less  cavity  drag  than  laminar  inflow. 

6.7  Wake  Flow  Characteristics 

The  flow  downstream  of  a resonating  cavity  is  highly 
disturbed.  As  discussed  previously,  the  shear  layer  above 
the  cavity  oscillates  about  the  cavitv  trailing  edge  and 
periodically  ejects  fluid  mass  into  the  downstream  flow  field. 
Thus,  there  are  both  mean  (or  average)  and  unsteady  velocity 
components  in  the  downstream  flow  field.  The  mean-velocity 
profile,  which  shows  a velocity  deficit  downstream  of  the 
cavity,  provides  information  about  the  steady-state  drag  and 
the  wake  thickness,  i.e.,  the  region  above  the  cavity  within 
which  the  wake  flow  represents  a disturbance  to  the  mean  flow. 

Velocity  profiles  on  the  cavity  centerline  were  measured 
with  the  aft  total  pressure  rake.  Figure  69  shows  results 
for  Moo  = 1 and  M„  = 2;  both  profiles  contain  data  on  cavities 
with  a length-to-depth  ratio  of  2.3  < L/D  < 5.5.  Qualitatively, 
the  Mach-2  flow  profile  indicates  a much  thinner  disturbance 
region  for  the  downstream  flow.  This  is  consistent  with  the  » 
smaller  dimensions  of  the  approaching  boundary  laver  at  Mach  2; 
it  is  also  consistent  with  the  lower  amplitudes  of  the  resonant 
mode  levels  at  Mach  2,  as  documented  in  Figs.  5A  through  56 
A free  shear  layer  above  a cavity  is  known  to  become  more 
stable  (and  stiffer)  with  increasing  supersonic  Mach  numbers. 
Hence,  we  would  expect  the  induced  mean  drag  to  decrease  with 
increasing  Mach  number. 


*MeGregor,  1969,  in  his  thesis  on  cavity  drag  under  thick  shear 
layer  uses  the  cavity  m<j>uth  area  successfully  as  a normalising 
parameter  into  the  drag  Coefficient. 

•j* 

The  term  "long"  indicates  that  the  length  of  the  cavity  is 
greater  than  the  width,  by  a factor  of  2 to  5. 


Smoothed  velocity  profiles  at  the  leading  edge  (forward 
rake)  and  downstream  of  the  trailing  edge  (aft  rake)  indicate 
the  growth  rate  of  the  shear  layer  above  a (resonating)  cavity 
for  Mach-1  flow  (Fig.  70).  The  boundary  layer  thickness  in- 
creases roughly  by  a factor  of  7 along  tne  cavity  length. 

By  using  the  momentum  thicknesses  determined  from  the  up- 
stream and  downstream  boundary  layer  rakes,  it  is  possible  to 
make  an  estimate  of  that  portion  of  the  cavity  drag  which  appears 
as  a momentum  deficit  in  the  boundary  layer.  This  can  only  be 
an  estimate  because  the  downstream  rake  was  located  only  on  the 
cavity  centerline;  therefore,  no  information  is  available  as 
to  the  variation  of  the  downstream  momentum  thickness  over  the 
cavity  width.  Furthermore,  the  downstream  rake  was  too  close 
to  the  cavity  trailing  edge  for  the  pressure  in  the  boundary 
layer  to  reach  equilibrium;  this  possibility  adds  another  un- 
certainty to  the  calculation.  In  addition,  in  order  to  compute 
the  total  drag,  the  wave  drag  component  would  also  have  to  be 
accounted  for,  which  has  not  been  done. 

If  Dd  Is  the  drag  component  due  to  the  momentum  deficit 
appearing  in  the  boundary  layer,  we  have 


m/. 


a,  H 


Dd  = 2(|pU*)  (6M2-«M1)  * W = |pU^(L-W)CD 


(6.3) 


where  SjyQ  and  6^2  are  the  upstream  and  downstream  momentum 
thicknesses,  respectively,  CD  is  a drag  coefficient  referenced 
to  the  cavity  mouth  area,  with  L the  cavity  length,  and  W its 
width . 

Table  II  lists  the  drag  coefficients  at  transonic  and 
supersonic  flow  speeds  for  three  L/D  ratios,  as  derived  from 
the  momentum  thickness  detriment.  The  transonic/supersonl.c 
drag  coefficients  are  lower  than  those  determined  in  the  small- 
scale  experiments  for  subsonic  speeds.  However,  it  should  be 
noted,  that  the  wave  drag  component  is  not  included. 
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TABLE  II.  CAVITY-DRAG  COEFFICIENTS  AT 
TRANSONIC/SUPERSONIC  SPEEDS 


n (in-} 


0.032 

0.034 

0.033 
0.062 
0.058 
0.055 
0 . 066 
0.064 
0.065 
0.078 
0.075 
0.081 
0.072 
0.082 
0.081 


12  (in.) 


0.463 
0.449 
0.399 
0.373 
0.501 
0.572 
0.612 
0.576 
0.429 
0.469 
0.605 
0.618 
0.692 
0 ,665 

0.476 
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6.8  Prediction  Approaches 

The  following  methods  are  proposed  for  the  prediction  for 
resonant  mode  frequencies  and  amplitudes  within  the  cavity  space, 

6.8.1  Frequency  Prediction 

Resonant  frequencies  can  be  accurately  nredicted  with  the 
Rossiter  equation  as  modified  by  Heller  et  al . (1970).  The 
present  research  data,  as  well  as  flight  test  data  reported 
by  Smith  et  al . (1974)  and  data  from  another  large-scale  wind 
tunnel  experiment  conducted  by  Maurer  (1974),  show  conclusively 
that  the  modified  Rossiter  equation  is  the  best  available 
expression  to  compute  resonant  frequencies  in  the  Mach  number 
range  0 , 5 < M < 3 • 0 . 
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Figure  52  shows  Eq.  6.1  in  a graphic  form.  Although  a 
dependence  of  mode  frequency  on  length-to-depth  ratio  exists 
(see  Sec.  6.5.3),  it  is  pointless  to  enter  this  complexity  in 
Eq.  6.1,  since  the  variation  of  mode  frequencies  for  shallow 
cavities  (L/D  > 2)  is  quite  small. 

6.8.2  Level  Prediction:  Resonant  Frequencies 

In  Fig.  57,  the  levels  of  inodes  1,  2,  and  3 are  presented 
for  the  leading-edge  and  the  trailing-edge  regions,  respectively, 
referenced  to  the  freestream  dynamic  pressure  as  a function  of 
freestream  Mach  number.  This  detailed  information,  which  also 
accounts  for  the  length-to-depth  ratio,  represents  a substantial 
improvement  in  the  accuracy  of  prediction  levels. 

6.8.3  Level  Prediction:  Broadband  Noise 

Broadband  noise  within  a Mach  number  range  of  0.8  < M^  < 2.0 
near  the  trailing  edge  can  be  predicted  on  the  basis  of  Fig.  62, 
where  a nondimensional  1/3-octave  band  spectrum  is  presented. 
Broadband  noise  was  found  to  be  a weak  function  of  the  length-to- 
depth  ratio.  Furthermore,  as  a rule  of  thumb,  broadband  noise 
levels  decrease  by  about  10  dB  in  a linear  fashion  towards  the 
leading  edge . 

6.8.4  Mode  Shape  Prediction 

The  complexity  of  mode  shape  prediction  is  documented  in 
Figs.  59  and  60.  These  mode  shapes  were  obtained  within  the 
free  cavity  volume  rather  than  at  the  cavity  floor.  Data  obtained 
along  the  cavity  floor  by  Shaw  et  at.  (1974),  which  used  a cavity 
with  a length-to-depth  ratio  of  4,  and  data  by  Heller  et  at. 

(1970)  are  presented  in  Fig.  71.  These  data  points  can  be 
described  by  an  empirical  equation,  which  is  developed  below. 

The  "standing  wave  pattern"  for  each  longitudinal  mode  in 
the  cavity  can  be  described  by  a cosine  relationship.  It  is 
assumed  that  the  pressure  intensities  increase  exponentially 
from  the  leading  to  the  trailing  edge.  Furthermore,  it  is 
assumed  that  the  physical  location  of  the  trailing-edge  bulk- 
head does  not  correspond  to  the  acoustically  effective  trailing- 
edge  bulkhead.  This  assumption  follows  both  from  physical 
considerations  and  from  the  pseudopiston  analysis  advanced  in 
Sec.  4.  From  the  data  presented  in  Fig.  71,  it  seems  that 
the  acoustically  effective  trailing-edge  bulkhead  is  located 
some  distance  upstream  of  the  actual  trailing-edge  bulkhead. 


k-  4 
(•;/  n 


■ i*‘>  <.*3Jm  V&ffGti 


Thus,  we  define  L*  = mL/(m+a),  where  L*  i s the  acoustically 
effective  cavity  length,  m is  the  mode  number,  and  a is  the 
empirical  constant  in  the  modified  Rossiter  equation.  The 
constant  is  -0,25  for  cavities  of  L/D  = 4. 

Thus,  the  mode  shape  in  the  cavities  should  be  described 
in  logarithmic  form  as 

x 


20  logCp/q^)  - 20  log(p/q  ) = 20  log 

L.E. 


a- 

e k | cos/ — • m * 180 

\ L*  i 


(6.4) 


Here  (p/q^)  is  the  rms  pressure  normalized  with  the  freestream 
dynamic  pressure;  (p/QoJl.E.  is  the  normalized  rms  pressure  at 
the  leading  edge;  a determines  the  exponential  growth  rate,  and 
x/L*  is  the  lengthwise  position  in  the  cavity,  normalized  with 
the  acoustically  effective  cavity  length  L* . 


The  agreement  of  the  proposed  analytical  experiment  with 
this  set  of  data  points  is  considered  rather  good,  considering 
the  vastly  different  experimental  conditions.  In  particular, 
mode-1  mode  shapes  as  measured  in  the  Air  Force  flight  test 
(Smith  et  at.,  1974)  are  extremely  well  described  by  the  proposed 
formulation.  In  a logarithmic  representation,  levels  at  the  nodal 
points  go  to  -<®.  In  this  representation,  a direct  exponential 
growth  was  assumed,  i.e.,  the  exponent  a was  taken  as  unity. 
However,  it  is  possible  that  the  actual  growth  rate  is  a function 
of  Mach  number  and  length-to-depth  ratio;  in  mode  3>  as  depicted 
in  Fig.  71,  the  growth  rate  seems  to  be  slower. 

It  should  be  possible  to  obtain  agreement  of  the  measured 
mode  shapes  as  presented  in  Figs.  59  and  60  by  assuming 
a different  growth  rate  a,  and  an  appropriate  value  of  the  length- 
to-depth  ratio  dependent  quantity  a.  However,  because  of  the 
highly  empirical  nature  of  the  above  approach  and  the  general 
insufficiency  of  data,  no  mode  shape  prediction  scheme  is  offered. 
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SECTION  7 

TEST  RESULTS:  OSCILLATION  SUPPRESSION 

7.1  Conceptual  Consideration 

The  analytical  considerations  presented  in  the  previous 
sections  suggest  several  potentially  successful  concepts  to 
minimize  the  amplitude  of  cavity  pressure  oscillation  or  to 
eliminate  che  occurrence  of  oscillations  altogether. 

Prom  our  understanding  of  the  physical  mechanism  of  pres- 
sure oscillations,  any  of  the  following  methods  should  affect 
the  oscillation  process: 

• Introduction  of  vorticity  into  the  shear  layer. 

• Suppression  of  the  feedback  mechanism. 

• Forcing  of  the  approaching  shear  layer  at  a frequency 
that  is  different  from  the  natural  oscillation  frequency. 

• Change  of  the  phasing  of  the  cavity  internal  pressure 
wave  propagation. 

Implementation  of  these  concepts  requires  geometric  changes; 
to  be  practical,  these  changes  should  not  reduce  the  usable 
cavity  volume.  Furthermore,  any  suppression  device  should  be 
simple  and  not  result  in  an  excessive  weight  or  drag  penalty. 

Many  configurations  were  tested  In  the  small-scale  wall-jet 
and  the  water  table  facilities  before  a few  selected  concepts 
were  evaluated  in  the  large-scale  tunnel. 

7.2  Unsuccessful  Concepts 

Figure  72  presents  several  concepts  that  did  not  reduce 
oscillatory  amplitudes,  or  that  sometimes  even  caused  higher 
levels  than  observed  in  the  basic  cavity.  Concepts  (a)  through 
(d)  were  evaluated  at  subsonic  speeds  (M  < 0.5);  concepts  (e) 
and  (f)  were  evaluated  at  subsonic  and  supersonic  speeds  (0.8 
< M < 1.5) • 

Cavity  Internal  Transverse  Spoilers  (Fig.  72a).  Jt  was 
thought  that  several  transverse  spoilers  would  interrupt  the 
process  of  shear  layer  oscillation  at  the  leading  edge  each  time 
that  the  shear  layer  dives  into  the  cavity  volume.  Apparently, 
however,  the  layer  assumes  a new  streamline  above  the  spoilers, 
and  the  oscillation  process  is  maintained. 
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Internal  Baffles  (Pig-  72b).  Insertion  of  perforated 
baffles  with  an  open  center  portion  was  thought  to  interrupt  the 
process  of  the  feedback  mechanism.  However,  the  prevention  of 
wave  propagation  along  the  cavity  rims  does  not  discourage  oscil- 
lat ions . 


Upstream  Spoiler  Cavity  (Fig.  72a),  Although  the  upstream 
spoiler  cavity  was  strongly  resonating,  and,  therefore,  should 
cause  a highly  disturbed  wake  flow  into  the  main  cavity,  no 
reduction  in  tone  level  was  observed. 

Oblique  Leading-  and  Trailing -Edge  Bulkheads  (Fig.  72d). 
Oblique  leading-  and  trailing-edge  bulkheads  were  thought  to 
offset  the  phasing  of  the  internal  pressure  wave  reflection  at 
the  leading-edge  bulkhead,  and/or  spread  the  generation  process 
of  the  upstream  traveling  pressure  wave  near  the  trailing-edge 
bulkhead  over  a certain  time  period.  Thus,  the  oscillatory 
process  would  be  randomized,  which  should  result  in  less  pro- 
nounced discrete  tones.  However,  no  reduction  in  tone  level 
was  observed. 
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Leading-Edge  Air  Entrainment  (Fig.  72e).  Forced  fluid 
entrainment  at  the  leading  edge  was  thought  to  counteract  the 
motion  of  the  recirculating  trapped  vortex.  However,  no  favorable 
effect  was  observed,  probably  because  only  a small  region  of 
recirculating  flow  occurred  close  to  the  leading-edge  bulkhead. 
This  left  the  main  flow  pattern  essentially  unaffected. 

Rounded  Leading-Edge  Bulkhead  in  Conjunction  with  Guidevane 
(Fig.  72 f) . It  was  suspected  that  the  guidevane  over  a sharp 
leading-edge  corner  (Fig.  72e)  would  not  permit  enough  fluid 
flow  to  enter  the  cavity  volume  to  counteract  the  recirculating 
internal  flow,  which  would,  thus,  discourage  oscillation.  There- 
fore, a more  favorable  flow  passage  in  the  forward  section  was 
constructed  by  rounding  the  leading-edge  bulkhead  with  an 
essentially  flat  guidevane  above  (Fig.  72f).  This  setup 
allowed  variation  of  the  gap  width  by  moving  the  vane  in  a 
vertical  direction,  thus  varying  the  amount  of  fluid  entrain- 
ment . 
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With  the  vane  in  the  lowest  position,  the  oscillation 
amplitudes  were  large;  in  fact,  discrete  tone  levels  exceeded 
those  of  the  basic  rectangular  cavity.  When  the  vane  was  moved 
into  the  freeetream  flow,  the  discrete  frequencies  disappeared, 
but  the  broadband  noise  level  increased  substantially;  thus, 
the  primary  purpose  of  reducing  cavity  noise  levels  was  defeated. 


7.3  Successful  Concepts 
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The  successful  reductions  of  oscillatory  amplitudes  were 
based  either  on  the  stabilization  of  the  shear  layer,  the 
prevention  of  the  periodic  mass-exchange  process,  or  on  a 
combination  of  both.  Shear  layer  stabilization  is  achieved  by 
two  means:  either  through  introduction  of  vorticity  into  the 

shear  layer  through  upstream  spoilers,  or  by  provision  of  an 
inherently  stabilizing  trailing-edge  shape.  Prevention  of  the 
mass-exchange  process  is  achieved  through  a detached  trailing- 
edge  cowl. 

Figure  73  contains  sketches  of  successful  configurations, 
and  Figs.  74  through  76  give  the  dimensional  information. 

These  configurations  are:  leading-edge  spoilers,  trailing-edge 

slant,  leading-edge  spoilers  together  with  trailing-edge  slant, 
and  detached  trailing-edge  cowl. 

These  configurations  were  evaluated  in  the  NASA  Lewis 
8 x 6 ft  tunnel  for  a variety  of  freestream  Mach  numbers  and 
cavity  length-to-depth  ratios.  The  setup  allowed  variation 
of  the  spoiler  angle  of  attack  over  a 360°  range.  The  spoilers 
were  always  counter-rotated  in  order  to  maintain  symmetry  of 
the  shed  vortex  wake  structure.  The  trailing-edge  cowl  (basically 
a thick  airfoil)  could  be  adjusted  in  its  vertical  position. 

Both  the  spoilers  and  the  cowl  could  be  adjusted  during  tunnel 
operation,  so  that  optimization  of  their  positions  was  possible. 
Figures  77  through  83  present  the  test  results. 

The  data  showing  the  effect  of  oscillation  suppression  devices 
have  been  presented  in  1/3-octave  bands.  If  narrowband  analysis 
had  been  employed  In  the  data  reduction,  the  difference  between 
peak  levels  of  unmodified  cavities  and  cavities  with  oscillation 
suppression  devices  would  have  been  even  more  pronounced. 


7.3.1  Shear  Layer  Stabilization 

Figure  77  shows  the  effect  of  deployed  spoilers  on  the 
sound  signature  near  the  leading  edge  (Sensor  3,  see  Fig.  36) 
at  a subsonic  Mach  number  (M  = 0.9)  for  length— to-depth  ratios 
of  2.3,  4.0  and  5,1.  Optimum  reduction  of  the  dominant 
mode  2 tone  at  250  Hz  Is  achieved  for  the  L/D  =2.3  cavity,  when 
the  spoilers  are  at  45°.  The  tone  level  is  reduced  by  20  dB. 

For  the  shallower  cavities,  spoiler  angles  between  45°  and  90°  have 
similar  effects,  and  the  tones  are  effectively  suppressed. 

Spoilers  deployed  at  45°  reduce  tone  levels  also  at  supersonic 
speeds  (M  = 1.5) 5 as  shown  in  Fig.  78. 
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FIG,  76.  DETACHED  Ttf A I L I NG-EPGE  COWL  AND  TRA I L I NG -EDGE  SLANT 
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ONE-THIRD  OCTAVE  BAND  CENTER  FREQUENCY  (Hz) 


FFFECT  OF  UPSTREAM  SPOILERS  AT  45°  ANGLE  OF  ATTACK 
ON  LEADING-EDGE  BULKHEAD  PRESSURE  SIGNAL:  Mm  = 0.9; 

L/D  = 4.0  (continuous  line  refers  to  unmodified 
cavity). 
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63  125  250  500  1000  2000  4000  8000  16,000 

ONE-THIRD  OCTAVE  BAND  CENTER  FREQUENCY  (Hz) 


EFFECT  OF  UPSTREAM  SPOILERS  AT  45°  ANGLE  OF  ATTACK 
ON  LEADING-EDGE  BULKHEAD  PRESSURE  SIGNAL:  = 1.5; 

L/D  = 2.3  (continuous  line  refers  to  .unmodified 
cavity). 


Figures  79  through  8l  show  the  effect  of  a 45°  trailing- 
edge  slant  as  measured  with  Sensor  3 (see  Fig.  36).  Substantial 
tone  reductions  are  achieved  at  subsonic  (M  = 0.8,  Fig.  79), 
transonic  (M  = 1.2,  Fig.  80),  and  supersonic  (M  = 1.5,  Fig. 

8l)  speeds.  However,  effective  tone  level  reduction  for 
the  deepest  cavity  (L/D  = 2.3)  requires  the  addition  of  deployed 
spoilers;  then,  mode  2 tone  levels  are  reduced  by  5 dB  (subsonic 
Fig.  79a)  and  by  15  dB  (transonic.  Fig.  80a,  and  supersonic. 

Fig.  81a). 

Appendix  C presents  the  physics  and  theory  of  flow  stabiliza- 
tion through  a trailing-edge  slant. 

7.3.2  Neutralization  of  Mass-Exchange  Process 

In  the  absence  of  leading-edge  spoilers,  the  performance  of 
the  trailing-edge  slant  is  improved  by  the  addition  of  a detached 
cowl,  as  documented  in  Figs.  82  (M  ■ 0.8)  and  83  (M  = 1.5). 

This  configuration  has  an  impressive  suppression  potential 
especially  for  the  relatively  deep  cavity  (L/D  = 2.3),  as  shown 
in  Fig.  82  and  83.  However,  the  position  of  the  cowl  is 
critical.  For  example,  at  M = 0.8,  the  optimum  cowl  position 
is  2 in.  above  the  surface  plane.  In  fact,  an  additional 
15  dB  were  gained  by  moving  the  cowl  from  the  zero  position  to 
the  +2-in.  position.  At  supersonic  speed,  the  optimum  cowl 
position  seemed  to  be  about  0.5  in.  below  the  surface  plane. 

The  reduction  of  the  mode-2  tone  is  more  than  30  dB . Appendix  D 
discusses  the  physics  of  the  flow  about  a detached  cowl. 

7.4  Consecutive  Cavities 

On  the  basis  of  the  successful  reduction  potential  of 
slanted  trailing  edges  in  conjunction  with  upstream  spoilers, 
a subsonic  flow  experiment  was  conducted  using  the  small- 
scale  wall-jet  facility.  In  this  experiment,  pairs  of  triangular 
spoilers  were  employed  upstream  of  both  cavities,  and  both 
cavity  trailing-edge  bulkheads  were  slanted. 

Figure  84  presents  results,  which  compare  1/10-octave  band 
spectra  on  the  forward  cavity  leading-edge  bulkhead  for  the  basic 
double  cavity  and  the  modified  double  cavity  for  an  external 
flow  speed  of  730  ft/sec.  The  dominant  tone  amolit.ude  is  re- 
duced by  more  than  30  dB.  However,  the  Inflow  into  the  cavitv 
configuration  was  laminar;  therefore,  part  of  the  attenuation 
must  be  attributed  to  tripping  the  boundary  layer,  thus  making 
the  inflow  turbulent,  which  is  an  Inherently  quieter  condition. 

Information  was  not  obtained  within  the  aft  cavity,  but 
pronounced  tones  were  not  observed  in  these  exneriment s , which 
implies  the  effectiveness  of  the  suppression  concept. 
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EFFECT  OF  TRA I L l NG-E DGE  SLANT  WITH  AND  WITHOUT 
UPSTREAM  SPOILERS  ON  LEADING-EDGE  BULKHEAD  PRESSURE 
SIGNAL:  Mo,  = 0.8;  L/D  = 2.3  (continuous  line  refers 

to  unmodified  cavity;  dashed  line  refers  to  trailing- 
edge  slant  only;  dotted  lino  refers  to  combination  of 
trai 1 i ng-edge  slant  and  upstream  spoilers). 
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FIG.  79b.  EFFECT  OF  TRAIL  I NG-EOGE  SLANT  ON  LEADING-EDGE 
BULKHEAD  PRESSURE  SIGNAL:  = 0.8;  L/D  = 40 

(continuous  line  refers  to  unmodified  cavity). 
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FIG.  79C. 


EFFECT  OF  TRAILING-EDGE  SLANT  ON  LEADING-EDGE 
BULKHEAD  PRESSURE  SIGNAL:  Mro  *0.8;  L/D  -5.1 

(continuous  line  refers  to  unmodified  cavity). 
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FIG  80a  EFFECT  OF  TRA I L I NG- EDGE  SLANT  WITH  AND  WITHOUT 

UPSTREAM  SPOILERS  ON  LEADING-EDGE  BULKHEAD  PRESSURE 
SIGNAL:  M<»  = 1.2;  L/D  - 2.3  (continuous  line  refers 

to  unmodified  cavity;  dashed  line  refers  to  trailing- 
edge  slant  only;  dotted  line  refers  to  combination 
of  trail ing-edge  slant  and  upstream  spoilers). 
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FIG.  80b.  EFFECT  OF  TRAIT'TN'G-E'DG'E  'SLANT  ON  'L'FAOI  NG- EDGE  "BULK.H  EAD 
PRESSURE  SIGNAL:  MlX)  = 1.2;  L/D  = 4.0  (continuous 

line  refers  to  unmodified  cavity). 
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' FIG.  '80c.  EFFECT  OF  TRAILING-LDGE  SLANT  ON  LEADING-EDGE 
BULKHEAD  PRESSURE  SIGNAL:  Mro  = 1.2;  L/D  = 5.1 

(continuous  line  refers  to  unmodified  cavity). 
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FIG.  81a.  EFFECT  OF  TRA I L I NG- EDGE  SLANT  WITH  AND  WITHOUT 

UPSTREAM  SPOILERS  ON  LEADING-EDGE  BULKHEAD  PRESSURE 
SIGNAL:  Mra  = 1.5;  L/D  = 2.3  (continuous  line  refers 

to  unmodifTed  cavity;  dashed  line  refers  to  trailing- 
edge  slant  only;  dotted  line  refers  to  combination  of 
trail ing-edge  slant  and  upstream  spoilers). 
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FIG.  81b.  EFFECT  OF  TRA I L I NG- EDGE  SLANT  ON  LEADING-EDGE 
BULKHEAD  PRESSURE  SIGNAL:  Miu  = 1.5;  L/L  = 4.0 

(continuous  line  refers  to  unmodified  cavity). 


RELATIVE  LEVEL  (dB) 


i-  waif  is#  jiYii 


|?>  '*. 


31.5  63  125  2 50  000  1000  2000  4000  8000  16,000 
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FIG.  81c.  EFFECT  OF  TRAIL  I NG-EDGE  SLANT  ON  LEADING-EDGE 
BULKHEAD  PRESSURE  SIGNAL:  = 1.5;  L/D  = 5.1 

(continuous  line  refers  to  un’modi  f i ed  cavity). 
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31.5  63  125  250  500  1000  2000  4000  8000  16,000 

ONE- THIRD  OCTAVE  BAND  CENTER  FREQUENCY  (Hz) 


FIG  82a  EFFECT  OF  TRAIL  ING-EDGE  SLANT  IN  COMBINATION 
WITH  A DETACHED  COWL  ON  LEADING-EDGE  BULKHEAD 
PRESSURE  SIGNAL:  = 0.8;  L/D  = 2.3  (continuous 

line  refers  to  unmodified  cavity;  dashed  line 
refers  to  cowl  trailing  edge  located  at  cavity 
mouth  level;  dotted  line  refers  to  cowl  trailing 
edge  located  2 in.  above  cavity  mouth  level). 
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FIG.  82c, 


EFFECT  OF  TRAILING-EDGE  SLANT  IN  COMBINATION 
WITH  A DETACHED  COWL  ON  LEADING-EDGE  BULKHEAD 
PRESSURE  SIGNAL:  M*,  = 0.8;  L/D  = 5.1  (continuous 

line  refers  to  unmodified  cavity). 
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FIG.  83a, 


EFFECT  OF  TRA I L I NG- EDGE  SLANT  IN  COMBINATION 
WITH  A DETACHED  COWL  ON  LEADING-EDGE  BULKHEAD 
PRESSURE  SIGNAL:  = 1.5;  L/D  = 2.3  (continuous 

line  refers  to  unmodified  cavity). 
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FIG  83b.  FFFECT  OF  TRA I L I NG- EDGE  SLANT  IN  COMBINATION 
WITH  A DETACHED  COWL  ON  LEADING-EDGE  BULKHEAD 
PRESSURE  SIGNAL:  Mro  =1.5;  L/D  = 4.0  (continuous 

line  refers  to  unmodified  cavity). 
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EFFECT  OF  TRA IL I NG-EDGE  SLANT  IN  COMBINATION  WITH 
A DETACHED  COWL  ON  LEADING-EDGE  BULKHEAD  PRESSURE 
SIGNAL:  = 1.5;  L/D  = 5.1  (continuous  line 

refers  to  unmodified  cavity). 
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7.5  Drag  Induced  by  Suppression  Devices 

Using  the  small-scale  wall-jet  facility,  a rectangular 
(basic)  cavity  with  length-to-depth  ratio  of  2 was  measured 
for  drag  due  to  (1)  a trailing-edge  slant,  (2)  two  inclined 
triangular  spoilers  at  angle  of  attack  of  35° , and  (3)  a combina- 
tion of  these  two  measures.  Figure  85  shows  the  drag  as  a 
function  of  flow  speed  for  the  four  conf Lgurations  under  considera- 
tion. These  data,  which  pertain  to  a specific  configuration, 
indicate  that  the  trailing-edge  slant  scarcely  increases  the 
drag.  However,  spoilers  which  have  particular  dimensions  and 
are  deployed  at  a particular  angle  of  attack  raise  the  drag 
by  about  40?;  the  combination  of  trailing-edge  slant  and 
spoilers  raise  it  by  about  60%  above  that  of  the  basic  cavity. 

7.6  Evaluation  of  Suppression  Devices 

The  value  of  an  oscillation  amplitude  reduction  device  is 
measured  in  terms  of  its  effectiveness  over  a broad  range  of 
flight  regimes,  its  ease  and  simplicity  of  implementation  into 
an  existing  or  new  aircraft  design,  and  its  possible  adverse 
effects  on  aircraft  performance. 

7.6.1  Acoustic  Effectiveness 

Obviously,  deeper  cavities  (L/D  £ 2)  require  some  drastic 
measures  to  reduce  the  oscillation  amplitude,  while  shallower 
cavities  (L/D  > 2)  can  be  "quieted"  in  a fairly  simple  and 
straightforward  manner.  The  following  conclusions  can  be 
drawn  from  the  experimental  results,  which  are  valid  for  long 
(L/W  s 4)  and  shallow  (2.3  < L/D  < 5.1)  cavities. 

• At  subsonic  speeds,  two  slightly  separated  flat-plate 
rectangular  spoilers,  which  are  twice  the  height  of 
the  local  boundary  layer  thickness  and  1/10  as  long  as 
the  cavity  length  at  a 45°  angle  with  respect  to  the 
freestream  flow  direction,  effectively  reduce  the 
dominant  discrete-tone  levels  to  the  broadband  noise 
level  (in  1/3-octave  bands)  for  cavities  with  length-to- 
depth  ratios  of  2.3  < L/D  < ®.  The  trailing  edge  of  the 
spoiler  should  be  close  to  the  cavity  leading,  edge.  The 
two  spoilers  should  be  separated  laterallv  by  about  one 
third  of  the  cavity  width  and  placed  symmetrically  with 
respect  to  the  cavity  centerline.  (For  reference,  see 
Fig.  77.) 

• At  supersonic  speeds,  the  same  spoilers  still  show  a 
significant,  although  lesser,  degree  of  level  reduction. 
(For  reference,  see  Fig.  78.) 
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• Dramatic  level  decreases  are  achieved  by  a 45°  slanted 
trailing  edge,  whose  projected  vertical  dimension 
should  be  at  least  twice  the  vertical  dimension  of  the 
naturally  grown  shear  layer  at  the  location  of  the 
trailing  edge.  The  tralling-edge  slant  can  protrude 
above  the  cavity  surface  plane  one  fourth  to  one  third 

of  the  slant  length.  Although  not  specifically  exnlored, 
this  protrusion  may  be  superior  to  a level  termination 
at  the  surface  plane. 

• The  trailing-edge  slant  is  effective  at  subsonic  (M»  = 0.8), 
transonic  (Mm  = 1.2),  and  supersonic  (Moo  = 1.5)  speeds 
since  it  not  only  reduces  the  discrete-tone  levels  but 

also  large  portions  of  the  broadband  noise.  (For  re- 
ference, see  Figs.  79>  80*  and  8l.) 

• For  the  relatively  deeper  cavities  (L/D  = 2.3)  s employment 
of  the  spoilers  at  a deployment  angle  between  30°  and  45° 
brings  further  reduction  in  noise  levels.  In  fact,  spoilers, 
in  conjunction  with  a slanted  trailing  edge,  represent  the 
best  configuration  explored  within  this  investigation;  thus 
they  are  recommended  for  implementation  in  aircraft  design. 

• In  the  absence  of  upstream  spoilers,  the  performance  of 
the  slant  can  be  improved  by  a thick  inclined  airfoil 
(also  referred  to  as  cowl)  some  distance  upstream  of  the 
slanted  trailing  edge,  so  that  a channel  is  formed  be- 
tween the  surface  of  the  slant  and  the  lower  (flat) 
surface  of  the  airfoil.  For  best  results,  the  airfoil 
maximum  thickness  at  the  quarter  chord  should  be  about 
25  to  30&  of  the  length.  The  length  Itself  should  be 

of  the  same  dimensions  as  the  slant  length.  This 
arrangement  achieves  dramatic  reductions,  but  the  critical 
posiuion  of  the  airfoil  depends  in  an  as  yet  not  fullv 
understood  manner  on  the  Mach  number.  (For  reference, 
see  Figs . 82  and  83 . ) 

• The  combination  of  spoilers  and  slanted  trailing  edge 
also  seems  to  be  effective  when  employed  with  two 
consecutive  cavities.  (For  reference,  see  wig.  84.) 

7.6.2  Implementation  into  Aircraft  Design 

The  least  space-consuming  reduction  device  Is  the  double 
spoiler  upstream  of  the  cavity  leading  edge.  The  spoiler  could 
be  retracted  for  cruising  at  the  deployment  angle  of  30°  or  45°, 
and  only  deployed  when  needed.  Alternatively,  the  spoilers  can 
continuously  protrude  above  the  aircraft  surface  at  a 0°  angle; 
this  angle  can  then  be  changed  to  30°  or  45°  when  needed. 
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The  slanted  trai ling-edge  configuration  takes  more  space. 

However,  by  letting  the  slant  extend  above  the  surface,  the 
space  penalty  can  be  minimized.  A substantial  amount  of  near 
trailing-edge  volume  would  be  consumed  by  the  slant  and  detached- 
cowl  configuration.  a 

a 

7.6.3  Effects  on  Aircraft  Performance  ] 

I 

The  effect  on  aircraft  performance  of  any  of  the  proposed 
suppression  devices  would  be  minimal,  since  the  open  cavity 
already  represents  a major  generator  of  steady  and  unsteady 
drag. 


Although  a resonant  cavity  does  not  generate  more  drag 
than  a quiescent  cavity,*  a stabilized  flow  would  represent 
aerodynamically  a cleaner  flight  condition,  even  if  there  was 
some  small  increase  in  drag. 

Prom  this  viewpoint,  the  slanted  trailing  edge  represents 
an  optimum  configuration,  since  it  causes  no  real  increase  in 
cavity  drag  (see  Fig.  85).  Even  in  conjunction  with  the 
deployed  spoilers,  the  drag  Increase  is  less,  than  if  only  the 
spoilers  were  deployed. 


*This  Is  in  contrast  to  McGregor’s  (1969)  findings.  However,  In 
McGregor's  experiments,  the  boundary  layer  thickness  was  about 
the  same  dimension  as  the  cavity  depth  and,  furthermore,  re- 
flecting tunnel  walls  were  used  to  enhance  the  sound.- 
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the  cavity.  A stronr , essentially  plane  pressure  wave  is  pro- 
duced by  this  process  and  propagates  forward  within  the  cavity. 
Since  this  wave  is  supersonic  relative  to  the  external  flow, 
a forward  traveling,  oblique  external  wave  is  also  produced. 

The  internal  wave  reflects  off  the  front  bulkhead  and  returns 
to  the  rear  of  the  cavity.  The  motion  of  the  shear  layer  is 
controlled  by  the  forcing  Imposed  by  these  traveling  pressure 
waves  over  the  length  of  the  cavity  mouth.  The  phase  of  the 
shear  layer  Is  such  that  as  one  pressure  wave  reaches  the 
forward  bulkhead  and  reflects,  another  is  simultaneously  gene- 
rated at  the  rear.  In  this  manner,  the  process  sustains  itself. 


This  explanation  of  the  oscillation  mechanism  is  somewhat 
different  from  that  proposed  by  previous  investigators.  The 
process  is  generally  viewed  as  a feedback  mechanism.  Shear 
layer  interaction  with  the  trailing  edge  produces  a wave  that 
reaches  the  leading  edge  and  triggers  a disturbance  in  the  shear 
layer.  This  disturbance  then  convects  downstream  and  amplifies. 
It  may  take  the  form  of  sinuous  instability,  vortex  shedding,  or 
some  combination.  The  disturbance  interacts  with  the  trailing 
edge  to  generate  another  upstream  wave,  which,  again,  triggers 
the  shear  layer,  and  so  on.  Particularly  severe  oscillations 
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are  believed  to  occur  when  the  frequency  of  this  process  matches 
the  natural  frequency  of  standing  wave  modes  in  the  cavity,  thus 
causing  a resonance  to  occur. 

The  viewpoint  derived  from  the  water  table  tests  differs 
in  the  following  way.  The  upstream  wave,  which  was  generated  by 
the  trailing-edge/shear  layer  interaction,  and  the  corresponding 
downstream  reflected  wave  directly  control  the  shear  layer  phase 
ty  pressure  forcing  over  the  length  of  the  cavity,  not  .lust  by 
triggering  disturbances  at  the  leading  edge.  The  shear  layer 
dynamics  are  not  primarily  the  result  of  a free  instability,  or 
of  at',  instability  triggered  solely  at  the  leading  edge,  but  they 
are  a consequence  of  the  forcing  from  both  upstream  and  down- 
stream traveling  pressure  waves  within  the  cavity.  The  super- 
position of  these  upstream  and  downstream  waves  forms  the 
experimentally  observed  pressure-mode  patterns  in  the  cavity. 
Special  shear  layer  phenomena,  such  as  vortex  shedding,  are 
viewed  as  a manifestation  of  this  pressure  wave  forcing,  rather 
than  as  an  essential  part  of  the  oscillation  mechanism. 

8.3  Analytical  Work 

On  the  basis  of  the  water  table  flow  visualization  experi- 
ments, analytical  models  were  developed  for  the  wave  motion  in 
the  cavity  and  for  the  trailingc-edge  mass  addition  process. 

The  water  table  results  have  shown  that  an  essential 
feature  of  the  phenomena  is  the  periodic  addition  and  removal 
of  mass  at  the  rear  of  the  cavity.  This  is  caused  by  the  inter- 
action of  the  shear  layer  with  the  trailing  edge.  This  mass 
addition  and  removal  produces  an  effect  very  similar  to  replacing 
the  rear  bulkhead  of  the  cavity  with  an  oscillating  piston. 

This  "pseudopiston"  effect  explains  why  pressure  modes  on  the 
cavity  do  not  behave  as  if  the  rear  bulkhead  were  a hard  wall. 

An  analytical  model  was  developed  to  describe  the  pseudopiston 
behavior  in  terms  of  cavity  geometry  and  flow  parameters. 

A separate  analysis  was  developed  for  the  wave  motions  in 
the  cavity.  The  cavity  was  modeled  as  having  an  oscillating 
piston  at  the  rear  bulkhead  to  account  for  the  pseudopiston 
effect.  This  approach  also  treats  the  Interaction  of  the 
wave  modes  with  the  shear  layer  as  an  integral  part  of  the 
oscillation  process.  Solutions  for  upstream  and  downstream 
waves  traveling  between  a shear  layer  and  a solid  boundary  are 
superimposed  to  approximate  an  oscillating  cavity  configuration. 
At  high  Mach  numbers,  this  analysis  predicts  pressure-mode 


Bk 


\ 


Vp* ! ; 


shapes  similar  to  those 
the  qualitative  features 
hers,  the  theory  predict 
and  fails  to  predict  the 
work  shows  are  present, 
caused  by  the  omission  o 
in  the  analytical  model, 
ment  be  the  subject  of  a 
alytical  results,  which 
range,  will  be  available 


observed  experimentally;  in  fact,  all 
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8.4  Experimental  Results 
8.4.1  Basic  Cavity 

Small-scale  and  large-scale  wind  tunnel  experiments  yielded 
detailed  information  on  the  aeroacoustic  behavior  of  shallow 
rectangular  cavities  over  a length-to-depth  ratio  range  of 
2.3  < L/D  < 5.1  ar.d  the  Mach  number  range  of  0.2  to  2.0.  The 
previously  derived  law  of  the  Mach  number  dependence  of  the 
resonant  mode  frequencies  was  substantiated  over  the  Mach  num- 
ber range,  except  at  low  subsonic  Mach  numbers;  in  this  range, 
resonant  frequencies  were  observed  to  be  higher  than  predicted. 
Additional  research  is  necessary  in  this  area  to  resolve  this 
discrepancy . 

Cavity  internal  temperatures  were  measured.  Recovery 
factors  were  found  to  be  closer  to  unity  than  to  zero,  whereby 
shallow  cavities  exhibited  higher  recovery  factors  than  deep 
cavities . 

A weak  dependence  of  the  resonant  frequencies  on  the 
length-to-depth  ratio  was  found  for  cavities  in  the  range 

1.5  < L/D  < 5 at  subsonic  flow  speeds.  Conclusive  data  at  super- 
sonic speeds  could  not  be  obtained.  More  research  in  this 

area  is  necessary  to  incorporate,  the  length-to-depth  ratio 
effect  into  the  present  prediction  schemes.  Detailed  informa- 
tion on  the  Mach  number  dependence  of  resonant  mode  levels  for 
modes  1,  2,  and  3 was  obtained.  Levels  were  observed  to  peak 
at  transonic  speeds.  Levels  were  also  found  to  increase 
generally  from  the  leading-  to  the  trailing-edge  region. 

Cavity  internal  pressure-mode  shanes  were  determined  for 
the  first  three  modes  and  for  the  range  of  2.3  < L/D  < 5.1.  In 
all  cases,  definite  amplitude  minima  were  observed;  however,  they 
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were  qualitatively  displaced  in  the  upstream  direction.  In  all 
cases,  a pressure  maximum  appeared  at  the  leading-edge  bulkhead 
indicating  that  this  surface  acts  as  a "hard-reflecting  wall." 

In  contrast,  pressure  amplitude  maxima  were  observed  a small 
distance  upstream  of  the  trailing-edge  bulkhead;  this  indicates 
that  acoustically  the  trailing-edge  bulkhead  is  undefined 
because  of  the  violent  flow  patterns.  In  this  study,  mode 
shapes  were  obtained  that  differ  from  those  observed  in  Air 
Force  flight  tests,  however,  in  the  latter  experiment,  mode 
shapes  were  determined  on  the  cavity  floor,  while  in  the  former 
experiments,  mode  shapes  were  determined  within  the  cavity  volume. 
Thus,  dependence  of  mode  shape  on  the  vertical  position  in  the 
cavity  seems  to  exist;  this  issue  needs  clarification. 

A nondimens ional  broadband  spectrum  was  derived  from  the 
experimental  data.  All  spectra  fall  within  a ±4-dB  range  for 
cavities  with  a range  of  2.3  < L/D  < 5.1  and  in  a Mach  number 
range  of  0.8  < Mto  < 2.0.  However,  the  nondimensional  peak  fre- 
quency seems  to  be  a function  of  the  thickness  of  the  approach- 
ing boundary  layer. 

Experiments  on  consecutive  geometrically  identical  cavi- 
ties reveal  that  they  are  strongly  coupled  and  resonate  in 
phase.  The  Mach  number  dependence  of  mode  levels  for  consecutive 
cavities  was  found  to  differ  from  that  of  single  cavities.  How- 
ever, this  may  be  because  of  the  differing  length-to-width  ratio 
of  the  single  and  the  double  cavity  system.  While  resonant  fre- 
quencies are  essentially  determined  by  the  cavity  length,  mode 
levels  seem  to  be  affected  by  the  relative  cavity  width.  This 
problem  has  not  been  studied  and  requires  further  research. 

8.4.2  Oscillation  Suppression 

Of  the  various  ways  to  affect  the  oscillation  process, 
introduction  of  vortlcity  into  the  shear  layer  and  the  provision 
of  a slanted  trailing-edge  bulkhead  were  found  to  have  a stabi- 
lizing effect  on  the  external  free  shear  layer.  Oscillation 
amplitudes  can  be  minimized  solely  by  a slanted  trailing-edge 
bulkhead,  over  a Mach  nuniuer  range  of  at  least  0.8  to  2.0,  for 
cavities  with  length-to-depth  ratios  above  4,  Cavities  with  a 
length-to-depth  ratio  below  4 require  the  addition  of  upstream 
vortex  generators  (spoilers),  which  further  reduce  resonant 
amplitudes . 

A physical  explanation  of  the  stabilizing  effect  of  a 
trailing-edge  slant  is  available.  It  would  be  desirable  to 
Investigate  this  concept  further  through  extensive  experimental 
exploration . 


■w"^g'rr^Tr>rg^?,"^Tq?  '•■*-;  \^j~  (s’-iy. 


Steady-state  drag  increase  with  a slanted  trailing  edge  is 
minimal.  The  effective  drag  of  the  spoilers  depends  on  their 
size  and  angle  of  attack.  Drag  increases  up  to  50S  above  that 
of  the  basic  rectangular  cavity  have  been  observed. 

Oscillation  amplitudes  could  not  be  reduced  by  the  forced 
entrainment  of  boundary  layer  fluid  into  the  cavity  at  the 
leading  edge,  by  tilting  the  leading-  and/or  trailing-edge 
bulkheads,  by  implementation  of  an  upstream  "spoiler  cavity," 
or  by  cavity  internal  transverse  spoilers  and  baffles.  J 
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APPENDIX  A 


PSEUDOPISTON  ANALYSIS 


A.l  Overview 

It  is  believed  that  the  cavity  oscillation  phenomenon  is 
driven  by  the  periodic  injection  and  removal  of  mass  at  the 
trailing  edge.  This  process  has  an  effect  that  is  similar 
to  placing  a fictitious  piston  at  the  rear  bulkhead.  The 
oscillation  of  this  pseudopiston  sets  up  an  essentially  one- 
dimensional wave  pattern  in  the  cavity  that,  in  turn,  forces 
the  shear  layer.  The  deflection  of  the  shear  layer  at  the 
trailing  edge  is  responsible  for  uhe  periodic  injection  of 
flow  from  the  shear  layer  and  the  free  stream,  and  for  the 
periodic  removal  of  mass  from  th$  cavity.  Under  proper 
conditions  of  phase  and  amplitude,  the  behavior  of  the  pseudo- 
piston and  the  dynamics  of  the  cavity  are  related  so  that  a 
self-sustaining  feedback  mechanism  occurs.  Of  course,  the 
details  depend  on  freestream  conditions  and  cavity  geometry. 

The  pseudopiston  analysis  relates  the  strength  and  phase 
of  the  fictitious  piston  to  the  phase  and  amplitude  of  the 
shear  layer  at  the  trailing  edge,  and  to  the  freestream  flow 
conditions  and  the  cavity  geometry. 

A. 2 Analytical  Model  for  the  Pseudopiston  Process 

The  actual  process  of  mass  introduction  and  removal  at 
the  trailing  edge  is  complicated  and  messy.  Section  3.2.3  gives 
a general  explanation  of  this  process;  thus,  it  is  sufficient  to 
state  here  that  a given  mass  flow  is  introduced  or  extracted  at 
the  trailing  edge.  The  relation  between  this  mass-flow  rate  and 
the  shear  layer  structure  and  deflection  will  be  discussed  later. 
Our  present  concern  is  the  piston-like  behavior  of  this  mass  flow. 

The  pseudopiston  can  be  modeled  in  several  ways:  e.g.,  as 
a simple  mass  flow  model,  or  as  a thermodynamic  mass  addition 
model . 


If 


u 

r * 


A. 2.1  Simple  Mass  Flow  Model 

In  this  kinematic  approach,  the  motion  of  the  fictitious 
piston  is  related  directly  to  the  mass  flow  rate.  Once  intro- 
duced into  the  cavity,  the  fluid  is  assumed  to  have  the  same 
thermodynamic  properties  as  the  fluid  already  in  the  cavity:  namely, 
freestream  static  pressure  and  stagnation  temperature  (approxi- 
mately). Hence,  the  fluid  density  is  determined  by  P = P^/RTo . 
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Then,  if  mass  flows  are  equated,  v;e  obtain 


m(t)  = p DV  ( t ) , (A .1) 

c p 

or 

v*>  - f%i  • (*-2> 

The  piston  displacement  is  determined  directly  from  the 
specified  mass  addition  or  removal.  Steady  flow  assumptions 
are  used.  However,  no  allowance  is  made  for  the  effect  of  the 
unsteady  expansion  that  the  added  flow  must  undergo  to  achieve 
the  thermodynamic  conditions  in  the  cavity.  Furthermore,  there 
is  no  pressure-matching  condition  across  the  interface  between 
the  cavity  fluid  and  the  newly  added  fluid,  since  these  pressures 
were  assumed  to  be  approximately  equal.  This  model  is  illustrated 
in  Fig.  86. 

A. 2. 2 Thermodynamic  Mass  Addition  Model 

This  more  refined  model  accounts  for  the  important  thermo- 
dynamic and  dynamic  effects  of  the  mass  addition  process  as  it 
relates  to  the  behavior  of  a fictitious  piston,  or  pseudopiston. 
This  process  is  illustrated  in  Fig.  87a. 

Another  model,  which  is  really  equivalent,  is  shown  in 
Fig.  87b.  This  model  is  based  on  the  idea  that  there  is  a 
steady-state  mass  addition,  m , into  the  cavity,  and  that  mass 
addition  and  removal  are  perturbations  around  this  steady-state 
value.  The  effective  mass  addition  or  removal  is  given  by 
/he  = rtir  - rtis , where  fir(t)  is  the  actual  mass  addition  at  the 
trailing  edge;  it  is  related  to  the  instantaneous  shear  layer 
position  in  an  unprescribed  manner.  For  present  purposes,  me 
can  be  considered  a given  function  of  time.  Likewise,  the 
piston  displacement  can  be  considered  as  a specified  function. 
These  parameters  could  be  interrelated  by  matching  pressure 
and  displacement  across  the  piston  and  bv  specifying  the 
relation  between  the  mass  flow  and  shear  layer  displacement. 

The  reservoir  enthalpy  is  given  and  assumed  independent 
cf  time.  Its  value  may  represent  some  average  value  found 
in  the  shear  layer.  With  this  analysis,  we  can  assume  that 
the  enthalpy  of  fluid  injected  into  the  cavity  is  the  free- 
stream  enthalpy,  because  use  of  a recovery  factor  makes  a 
difference  of  only  a few  percent.  This  assumption  was  confirmed 
by  temperature  measurements  (see  Sec.  6.4). 
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FIG.  66. 


SIMPLE  MASS  FLOW  MODEL. 


The  size  of  the  mass  addition  region,  l,  is  small  compared 
co  the  length  of  the  cavity,  L,  which  means  that  the  time  scale 
for  the  oscillation,  l>/&c  , is  long  compared  to  the  time  scale 
for  adjustment  of  conditions  within  the  mass  addition  region, 

J,/^  . In  comparison  to  the  time  scale  of  oscillation,  the  mass 
addition  region  has  considerable  time  to  adjust  to  changes  in 
pressure,  mass  flow,  etc.,  which  are  produced  by  the  oscilla- 
tion cycle.  Therefore,  we  can  assume  that  significant  departures 
from  thermodynamic  equilibrium  do  rot  occur  in  the  mass  addition 
region.  This  Justifies  the  quasi-steady  treatment  of  the  mass 
addition  region  and  the  assumption  of  uniform  properties  through- 
out the  mass  addition  chamber.  Furthermore,  the  pressure  and 
energy  changes  resulting  from  flow  velocity  in  this  region  are 
neglected,  since  the  velocity  is  presumably  small  enough  to  make 
these  effects  insignificant. 


I- 


The  mass  in  the  imaginary  chamber  is 


p(t)DH(t)  = 


p djl 

o 


me(t)dt 


(A.  3) 


where  me(t0)  is  eq •■al  to  zero,  ^or  steady  flow,  iQ  is  not 
necessarily  the  steady-state  value  of  fc(t).  The  reference 
value,  p0,  is  a constant  density  in  the  cavity  corresponding 
to  conditions  at  time,  tQ. 

At  this  time  me(t)  and  i,(t)  can  be  considered  as  given 
functions  of  time.  Therefore,  the  density  in  the  chamber  is 


po  TTtry  + 


X me(t'dt 


CA.il) 


This  equation  determines  the  density  in  terms  of  quantities 
assumed  given.  The  equation  of  state  then  relates  pressure 
and  temperature: 

p(t)  = p(t)RT(t)  , 

p = p RT  , 

, *o  o o * 

and 

D = p RT 

for  steady-flow  conditions. 


(A. 5a) 
(A. 5b) 
(A. 5c) 
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If  we  neglect  the  effect  of  flow  velocity  in  the  chamber,  the 
energy  equation  Is 


AE  = 0 + W + E, 


(A. 6) 


where  AE  is  the  change  in  internal  energy,  o is  the  heat  addition 
(e.g.,  heat  transfer  to  or  from  system),  W is  the  work  done  cn 
the  fluid  volume  in  the  chamber  (e.g.,  by  the  piston),  and  EjvjA  is 
the  energy  addition  due  to  mass  addition. 

Assuming  that  significant  heat  transfer  does  not  occur  in 
the  chamber  during  an  oscillation  cycle,  let  Q *»  0, 

The  energy  equation  expressed  in  detail  is: 

t t 

E(t)-Eo  = - f pDfc(t)dt  +/  me(t)hQdt  . (A. 7) 

to  to 

Note  that  Ep-E(t0)  and  A0=A(tc)  are  reference  values  determined 
at  t = tgj  they  are  not  necessarily  the  corresponding  steady- 
state  values.  Nor  does  the  formulation  given  thus  far  reauire 
that  they  be  steady-state  values.  The  time,  tn,  is  chosen  so 
that  m(tQ)  = 0.  ° 

The  specific  internal  energy  is  e(t).  The  temperature  is 
T(t ) , and 


F(t)  = p(t)Ut)De(tO  = pAD(cvT  X)  , 


(A. 8) 


where  K is  a constant, 


Thus  , 


E(t)-E0  = p4DovT(t)  - P0D*ocv  T(tQ)  + DK( p.t-p^ ) 


s&SsaN 


Substituting  into  the  energy  equation. 


ptDovT  - P0^0DovT(tQ)  + KD(pfc-pofc0) 


,t 


-f  pDfcdt  +1  me(t)hQdt  . 


(A. 10) 


Differentiating 


(PJID)CgvT+K)  + CpilD)ovT  - -pD*  + AehQ  . 


Using  the  density  relation,  (Eq.  A.3')»  we  obtain 


( pDfc ) » m„  ♦ 


Thus, 


(A. 11) 


( A .12) 


me(cvT+K)  + (PqD^o)  + J roedt  cyT  + pD  l » ™ehQ  » 

t 


(A. 13) 


cv  {(poDV  I "1edt  * + ”1ecvT  + pD*  * ™e(ho“K^ 
l to 

where  Tq  is  tne  freestream  stagnation  temperature, 
Next, 

I.  *.  A!  ”edt  \ 


">e°pV  (A'1II) 


p « pRT  = RT  I °o  _0  4. 


(A. 15) 


&T,  I*. 


and,  therefore, 


pDi  -RT  |(Poi0D)  |+  | f medt]  . 


Substituting  into  Eq . A.14  gives 


(A. 16) 


t * * 

(P0^0)  + ft  ®edt]  * + ®ecv  + R(poV}I  + RI  { Q ”edt  T " VpTo  * 


from  which 


(A. 17) 


• . m yT„ 

m s e o 

* + 7 fr~. r + (y  - Df  t — : » 

l(poD*o+/  inedt)  P0D£0+J[  medt 

° ° (A. 18) 


r\  H 

where  R = cp  - cv,  and  y * ; therefore,  — * y - 1. 

The  above  is  a first-order  ordinary  differential  equation 
for  temperature,  T.  The  solution  is 


lY"1(poD1o  + / ™edt) 


(meYTo&Y~;Ldt  + cj  , (A.  19) 


where  C is  an  arbitrary  constant.  Using  Eq . A. 15,  the  pressure 
is,  therefore. 


P = § Y TqJTy  J medt  + C§  JTY  • 


(A. 20) 


The  above  expression  gives  the  pressure  in  terms  of  ?(t)  and 
me(t).  The  second  term,  containing  the  arbitrary  constant, 
corresponds  to  an  adiabatic  compression  without  mass  addition 
or  removal.  The  constant  is  determined  by  the  volume  and 
energy  of  the  gas  being  compressed. 
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Note  that  if  the  arbitrary  constant  is  set  equal  to  zero, 
and  if  it  is  further  assumed  that  me  is  directly  proportional 
to  i,  the  result  is 


5°sTo  • 


(A.  21) 


where  li.  is  C l, 
e s 

Mass  conservation  suggests 


PCD*  => 


(A. 22) 


Thus,  for  this  special  case,  we  have  constant  pressure  p » pcRTQ 
* pc(*pw).  The  flow  into  the  cavity  is  treated  as  being  quasi- 
steady. Therefore,  with  this  mass  conservation  assumption,  we 
have  the  same  result  that  would  be  obtained  using  the  simple 
mass  flow  model  for  the  pseudopiston.  Thus,  the  second  model 
developed  in  this  section  is  consistent  with  the  first  when  we 
make  the  same  assumptions. 

A. 3 Discussion  of  Mass  Flow,  Shear  Layer  Deflection,  and 

Pseudopiston  Motion 

The  analytical  model  for  mass  addition  at  the  rear  of  the 
cavity  relates  the  pressure  of  a fictitious  piston  to  the  piston 
displacement  and  the  mass  flow  rate  behind  the  piston.  The 
time-dependent  mass  flow  rate  depends  on  the  shear  layer  deflec- 
tion process  at  the  trailing  edge.  The  shear  layer  dafleation 
is  itself  dependent  on  the  piston  displacement.  This  dependence 
is  determined  by  solving  the  problem  of  shear  layer  forcing  by 
waves  in  the  cavity  which  are  generated  by  the  piston.  This 
solution  involves  the  dynamics  of  the  cavity  as  a whole,  and  it 
is  discussed  in  Sec.  4. 

From  this  analysis,  a prescribed  piston  displacement,  JL  ( t ) , 
will  determine  the  pressure,  p(t),  on  the  front  face  of  the 
piston  and  the  shear  layer  amplitude  near  the  trailing  edme, 
r'T,E.(t).  the  mass  flow  relation,  me  [n(t)],  were  also  known, 
it 'would  be  possible  to  match  displacements  and  pressures  with 
the  analytical  model  for  the  pseudopiston  and  determine  the 
oscillatory  behavior  in  this  manner.  The  actual  process  of 
mass  addition  at  the  trailing  edge  is  extremely  complex,  and 
perhaps  it  cannot  be  modeled  analytically  to  the  required  accuracy, 
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Nevort  lie  loss  , the  present 
.-.eneral  dependence  of  tin.1 
and  flow  parameter:’,. 


,o  present  analysis  Is  sufficient,  to  indicate  the 


mass  addition  procesr 


a re o me try 


In  relatinr,  !>(<■)  to  rw.i.-  U5 
two  problems:  the  phase  relation, 


and  P ( t, ) , there  are  essentially 
and  the  functional  dependence 


with  time.  It  seems  reasonable  to  believe  that  m (t)  and  n l - } 
could  be  directly  in  phase  with  a negative  deflection  producinr 
a positive  mass  flow.  Furthermore , from  a nseudopiston  viewpoint, 
the  mass  flow  causes  the  displacement  so  that,  we  mirht  expect 


be  in  phase  wr 


Figure  88  illustrates  this  rrunhically . 
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FIG.  88.  COMPAkISON  OF  THE  TIME  BEHAVIOR  OF  THE  SHEAR  LAVER 
DISPLACEMENT,  THE  MASS  FLOW  RATE,  ANO  THF.  PSEUDO- 
PISTON MOTION, 

Although  the  mass  addition  process  is  complex  and  involves 
the  dynamics  of  a stagnation  point  within  an  unsteady  shear 
layer,  the  assumption  of  an  in-phase  relation  between  n and  rvp.E. 
seems  reasonable.  The  relation  between  n-p.E.  and  ) depends 
on  the  dynamics  of  the  whole  cavity,  and  it  is  not  necessary  that 
_ and  2. ( t ) be  exactly  in  phase.  There  could  he  a nhase  shift, 
because  of  the  pressure-match! nr,  relation  across  the  psoudooiston , 
When  waves  are  forced  by  a piston  in  a rectangular  enclosure,  the 
pressure  and  displacement  at  the  piston  are  in  phase  and  the 
pressure  and  velocity  are  out  of  phase,  so  that  no  net  work  is 
done.  In  the  cavity  problem,  enerrv  is  radiated  from  and  adaed 


to  the  system.  Presumably,  this  Is  manifested  as  a phase  shift 
between  pressure  and  velocity  that  is  imposed  by  the  effects 
of  shear  layer  motion. 

The  greatest  source  of  uncertainty,  however,  involves  the 
functional  relation  between  me  and  tvj>,  E . • J f shear  layer 
deflection  is  small  compared  to  shear  layer  thickness,  then  a 
linear  relation  between  mass  flow  and  displacement  is  likely. 

If  displacements  of  the  order  of  the  shear  layer  thickness  are 
expected,  then  the  functional  dependence  is  more  uncertain  and 
is  probably  complicated. 

The  expression  for  pressure,  determined  from  the  previous 
section,  is 

P - ~ Y T0JfY  / *Y-1me  dt  + . (A. 23) 

This  expression  is  nonlinear  in  £(t).  If  displacements  of  the 
chamber  are  small  compared  to  its  steady-state  size,  i.e., 
[x(t)-£g]  / [ 5. ( t ) ] <<  1,  then  the  expression  can  be  linearized. 
This  condition  corresponds  to  shear  layer  displacements  that  are 
small  compared  to  shear  layer  thickness.  If,  on  the  other  hand, 
shear  layer  displacements  are  of  the  order  of  its  thickness,  the 
expression  for  p cannot  be  linearized  with  respect  to  5..  In 
this  case,  an  alternative  for  analytical ^purposes  Is  to  assume 
a mass  flow  dependence  on  displacement,  me[£(t)],  which  rives  a 
tractable  functional  form  for  the  pressure  In  terms  of  displace- 
ment. Such  an  approach  is  somewhat  arbitrary  because  It  contains 
the  assumption  that  an  extremely  complex  process  will  adjust  its 
behavior  to  some  simple  form.  However,  it  receives  support 
because  the  rest  of  the  cavity  would  like  to  behave  In  a simple 
modal  way.  An  example  of  this  approach  was  previously  illustrated 
by  the  assumption  me  « i,  for  the  quasi-steady  flow  model. 
Similarly,  more  refined  forms  can  be  assumed  to  rive  the  pres- 
sure in  a desired  functional  form. 

Now  let  us  return  to  the  possibility  of  linearizing  the 
pressure  expression  when  shear  layer  displacements  are  small. 

Let  nT  £ be  small  and  (£-Jls)/S,  <<  1.  Write  H = and 

hm  tt  ■ en,»  where  e is  the  maximum  of  the  small  amplitude. 

Let  me  ■ p £US£erl>  where  the  PsAuC£  is  the  product  of  density 

in  the  shear  layer  and  a characteristic  shear  layer  velocity. 

Then, 
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(A. 24) 


P - Y | V~7/  ^'lAedt  + C | *~Y  * 


Substituting, 


■>  ‘ srar2  psiu8i/niat  ♦ c§  ls"V(1  + + <*-25) 


(1  + tip*  ■ eY  Jtn(X  + £f,p  eYle>ll  + 1 + YEti  +0[c!J  , 


P = c 5 *~Y+ 


uoo  /* n,dt  + C§4  “yyA,I  + 


(A. 26) 


DT~  ''ai.valj  "l 


D s ’"1 


Y*n  + o[e2]  . 


Choose  the  arbitrary  constant  C so  that 


(A. 27) 


m 


i*n» 

Vt*t& 

>-  - «u 


C § rY  = Pc  85  P"  =>  C = Pc  | . 


(A. 28) 


Therefore , 


P = Pc  + e 


ryto 


psAUsjt/\dt  + PCY^)+  OCe2]  • 


(A. 29) 


Note  that  to  order  e the  pressure  changes  from  mass  addition  and 
piston  displacement  are  uncoupled  and  simply  added.  Cross  terms 
will  appear  at  0[e2].  As  mentioned  before,  the  relation  between 
hi  and  ?. j results  from  determining  the  dynamics  of  the  whole 
cavity  forced  by  a piston,  which  is  an  essentially  separate 
problem. 

In  practice,  there  is  really  no  way  to  determine  theoretically 
&s,  which  is  the  spanwise  size  of  the  mass  injection  and  mixing 
region  at  the  rear  bulkhead.  The  product,  «,gD,  is  the  volume  of 

this  region.  Roughly,  Jl  = 8D , where  8 is  some  constant  <0[l]. 
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For  shallow  cavities,  It  Is  possible  that  6 may  be  relatively 
independent  of  overall  cavity  geometry  and  external  flow 
conditions . 
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Since 


= p RT  , let  us  write 
c Kc  o’ 


- 

/ p a , 

“ 

_x_ 

gD2 

(57) 

“s*  J 

f n dt  + yA, 
1 1 

(A. 30) 


Assume  the  velocity  in  the  shear  layer  is  ^uU  , where  a ~0[1]. 

This  should  be  the  velocity  on  the  dividing  streamline  of  the 
shear  layer  near  the  trailing  edge,  i.e.,  the  streamline  dividing 
mass  entrained  from  above  and  the  mass  entrained  from  below.  If 
the  stagnation  enthalpy  is  unchanged  through  the  shear  layer, 
then 

ho  = 2pTo  = h + = V + |v2  , (A. 31) 


and 


RT 


s£ 

si 


RTo(l-V2/2cpTo) 


P. 

RT 

o 


Therefore , 


ps&  „ 1 a 1 

pc  l-Va/2cpTo  l-a2U2/8cpTQ 


(A. 32) 


(A. 33) 


(A. 34) 


194 


Since  aw  * /yRTo  , R = cp-cv,  and  y = cp/cv,  this  becomes: 


si 


1 + M2 
2 » 


i + f1  - x)3^ 


M' 


(A. 35) 


Finally, 


Thus, 


= M a 
00  00  00 


(A. 36) 


p-p 

ep.  • X 


c _ 


Y-,/a  (HTo) 1/2 


^l+l^i  M, 


l+d-Tf)  ^ Ma 


M. 


«(f)/eni 


dt  + yea, • 


(A. 37) 


This  result  gives  the  perturbation  pressure  of  the  pseudopiston 
in  terms  of  effective  piston  displacement  and  shear  layer 
displacement  based  on  a linearization  of  the  pressure  relation 
derived  earlier.  The  linearized  model  implies  a specific  set 
of  assumptions  about  the  working  of  the  pseudopiston  mechanism. 
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TABLES  OF  ROOTS 
Root  A 

(to  obtain  A*  replace  i with  -i  everywhere) 


M 


0.1  K 
±T 
+B 
C 


±T 

+B 

C 


S = 1 .0 


S = 0.5 


S = 0.1 


P 


P 


-10.035  + 110.000 
-10.010  + 110.025 
9.975  + 110.060 
- 0.050  - 1 0.050 

- 2.13-4  + 1 1.984 

- 2.012  + 1 2.105 
1.862  + 1 2.251 

- 0.251  - 1 0.234 


-5.017  + 15.000 
-5.005  + 15.012 
4.987  + 15.029 
-0.050  - 10.050 

-1.170  + 10.948 
-1.107  + 11.002 
0.886  + 11.214 
-0.258  - 10.209 


-1.103  + 10.958 
-1.100  + 10.960 
0.956  + 11.105 
-0.052  - 10.045 

-0.319  + 10.281 
-0.310  + 10.289 
0.254  + 10.322 
-0.177  - 10.156 


1.0 


1.5 


3.0 


K 

±T 

±B 


K 

±T 

±B 


- 1.348  + 1 0.919 

- 1.107  + 1 1.H9 
0.677  + 1 1.487 

- 0.506  - 1 0.345 

- 1.142  + 1 0.587 

- 0.807  + 1 0.831 

0.181  + 1 1.228 

- 0.693  - 1 0.356 


K 

+T 

±B 

C 


- 1.060  + 1 0.234 

- 0.534  + 1 0.465 

- 1.113  + 1 0.728 

- 0.900  - 1 0.199 


-0.768  + 10.504 
-0.658  + 10.589 
0.341  + 10.811 
-0.455  - 10.299 

-0.642  + 10.338 
-0.493  + 10.441 
0.052  + 10.654 

-0.610  - 10. 321 

-0.587  + 10.166 
-0.370  + 10.264 
-0.732  + 10.496 
-0.789  - 10.223 


-0.202  + 10.159 
-0.187  + 10.172 
0.087  + 10.200 
-0.306  - 10.241 

-0.172  + 10.103 
-0.151  + 10.117 
-0.023  + 10.156 
-0.428  - 10. 256 

-0.161  + 10.066 
-0.134  + 10.079 
-0.274  + 10. 190 
-0.532  - 10. 218 


P 


m 


I 


NOTES:  T = aD  and  T2  = K2-S2,  where  K = kD  and  S = — and  a - a . 

0 00 

B = 3D  and  B2  = Q-K2,  where  Q = (MK  + aS)2  . 

c s 

c = -2., where  C = £ is  the  phase  speed.  Thus  Cp  = ^ . 

P ac  ? 

Dispersion  Relation:  Q2T2  sinh2T  + S (Q-K  ) cosh  T - 0 . 
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APPENDIX  C 

AERODYNAMIC  FLOW  STABILIZATION  THROUGH  A SLANTED  TRAILING  EDGE 

This  section  will  present  the  physical  and  theoretical 
basis  for  slanting  the  upper  portion  of  the  cavity  trailing- 
edge  bulkhead  as  a means  of  reducing  pressure  oscillations 
in  the  cavity.  Although  this  is  a simple  modification,  a 
considerable  reduction  of  levels  has  been  observed  experimentally. 
Slanting  the  edge  seems  to  be  more  effective  than  other  modifica- 
tions of  roughly  the  same  size,  such  as  rounding  the  trailing 
edge.  By  studying  the  reasons  for  the  success  of  this  modifica- 
tion, we  can  reach  a deeper  understanding  of  the  shear  layer 
impingement  process  at  the  trailing  edge. 

First,  it  is  appropriate  to  recall  the  discussion  in 
Sec.  3.2  entitled  Steady-Flow  Considerations.  It  was  argued 
that  the  shear  layer  must  have  a stagnation  point  at  the 
trailing  edge,  and  that  part  of  the  shear  laver  flow  must  be 
returned  to  the  cavity  (to  replenish  entrained  fluid)  and  the 
rest  must  go  downstream  with  the  mean  flow.  This  view  of  the 
steady-flow  configuration  was  then  used  to  explain  the  mass 
addition  and  removal  processes  at  the  trailing  edge,  when  the 
shear  layer  behaves  in  an  unsteady  manner. 

The  trailing-edge  stagnation  point  must  lie  within  the 
shear  layer  itself,  and  it  is  expected  that  the  stagnation 
streamline  will  be  near  the  center  of  the  shear  layer,  where 
the  vorticity  is  greatest.  Therefore,  it  is  of  interest  to 
examine  in  some  detail  the  stagnation  flow  when  the  flow  has 
vorticity.  The  physics  of  this  process  is  examined  below  in  a 
simplified  form.  The  effects  that  have  been  neglected  can  be 
expected  to  alter  the  flow  details,  but  they  will  not  affect 
the  conclusion. 

Consider  the  nature  of  the  local  region  around  the  stagna- 
tion point,  i.e.,  a region  of  dimensions  less  than  the  shear 
layer  thickness.  It  is  appropriate  to  assume  a two-dimensional 
incompressible  flow  in  this  region.  Shear  layer  solutions, 
either  laminar  or  turbulent,  show  that  the  shear  is  nearly 
constant  (constant  vorticity)  In  the  middle  region  of  the 
layer.  Therefore,  an  irrotational  stagnation-point  flow  and 
a simple  shear  flow  will  be  combined  to  find  the  actual  stream- 
line pattern  near  the  stagnation  point.  For  the  present  case, 
the  stream  functions  can  be  simply  added  together.  Generally, 
stream  functions  for  rotational  and  Irrotational  flows  cannot 
be  simply  combined,  because  the  equations  of  motion  are  nonlinear 
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due  to  the  presence  of  the  convection  terms  when  the  flow  has 
vorticity.  However,  this  simple  addition  can  be  done  when 
certain  conditions  are  satisfied.  The  flow  must  be  two- 
dimensional,  the  vorticity  must  be  constant  everywhere,  and 
the  rotational  and  irrotation  parts  must  independently  satisfy 
the  same  boundary  geometry.  In  the  present  case,  these 
conditions  are  satisfied  as  long  as  the  region  of  validity  is 
restricted  to  the  immediate  neighborhood  of  the  stagnation 
point  and  the  orientation  of  the  simple  shear  flow  is  correctly 
chosen. 

The  stream  function  of  a stagnation  point  in  a constant 
vorticity  flow  is 


axy  + \ by*, 


(C.l) 


where  the  first  term  is  the  stagnation  point,  and  the  second 
is  the  shear  flow.  The  corresponding  velocity  components  are 


If  * ax  + by, 


(C.2) 


v - - It 

ax 


(c.3) 


The  constant  a is  related  to  the  strength  of  the  stagnation- 
point  flow,  while  the  constant  b is  related  to  the  strength 
of  the  shear.  The  vorticity  is 


„ 3u  K 

c • - ay  b 


(C.4) 


The  streamline  pattern  is  shown  in  Fig.  89,  where  the  values 
of  the  constants  have  been  absorbed  into  the  coordinates. 
Nevertheless,  Fig.  89  gives  a realistic  picture  of  the  actual 
flow  pattern.  The  important  conclusion  is  that  the  effect  of 
shear  is  to  change  the  direction  of  the  stagnation  streamline 
(and,  presumably,  the  centerline  of  the  shear  laver),  so  that 
the  flow  must  impinge  on  the  wall  (cavity  rear  bulkhead)  at  an 
oblique  angle.  It  cannot  impinge  normal  to  the  wall,  as  in  the 
case  of  no  shear. 


i The  fact  that  the  flow  impinges  at  an  angle  has  a simple 

; physical  explanation.  Shear  flow  has  a velocity  gradient; 

i.e.,  velocities  are  higher  hn  the  freestream  side  of  a shear 
f layer  than  on  the  cavity  side.  As  the  flow  impinges  on  a wall, 

I the  streamlines  must  curve,  which  introduces  centrifugal  pres- 

sure forces.  If  the  flow  were  symmetrical  about  the  stagnation 
streamline,  the  opposing  pressures  would  not  balance  because 
! the  velocity  is  higher  on  one  side.  However,  if  the  streamline 

radii  of  curvature  are  made  larger  on  the  high  velocity  side 
and  smaller  on  the  low  velocity  side,  then  a balance  can  be 
► achieved. 

[ ; This  conclusion  has  definite  implications  for  the  phenomena 

of  cavity  oscillations.  In  order  to  suppress  oscillations, 

1 a steady  flow  must  be  physically  possible.  Figure  90a  depicts 

t the  behavior  of  the  shear  layer  in  a conventional  cavity  configura- 

tion. Note  that  the  information  about  the  angle  of  the  stagna- 
tion streamline  near  the  rear  bulkhead  was  incorporated  into 
Fig.  90a,  since  this  must  be  true  in  a steady  flow.  In  order 
to  meet  this  condition,  a considerable  curvature  of  the  shear 
layer  over  the  cavity  mouth  is  required.  This  curvature  will 
cause  the  freestream  flow  to  produce  static  pressure  variations 
over  the  cavity  mouth  that  cannot  be  balanced  by  the  pressure 
within  the  cavity,  which  is  uniform  to  first  aoproximation  (re- 
circulation velocities  are  much  less  than  the  freestream  velocity). 
Hence,  an  unsteady  motion  of  the  shear  layer  will  result. 
Alternatively,  if  the  shear  layer  is  assumed  to  be  straight  so 
that  unbalanced  pressures  do  not  occur  over  the  cavity  mouth, 
the  impingement  angle  for  steady  flow  cannot  be  satisfied  and 
the  flow  in  the  immediate  vicinity  of  the  rear  stagnation  point 
must  be  unsteady.  This  type  of  reasoning  suggests  that  steady- 
flow  configurations  over  conventional  unmodified  cavities  are 
difficult,  if  not  impossible,  to  achieve.  When  this  basic 
tendency  towards  unsteadiness  is  combined  with  the  ability  of 
cavities  to  sustain  a modal  unsteady  pressure  field,  which  can 
be  coupled  to  shear  layer  motion,  it  is  easy  to  see  why  the 
suppression  of  pressure  oscillations  is  so  difficult. 


When  the  portion  of  the  trailing-edge  bulkhead  that  | 
interacts  with  the  shear  layer  is  slanted,  a steady-flow  solu-  | 
tion  is  much  easier  to  achieve.  If  the  slant  angle  Is  properly  * 
chosen,  the  shear  layer  can  remain  essentially  straight  over  I 
the  cavity  mouth  and  satisfy  the  proper  impingement  angle  as  ' $ 
well.  This  is  illustrated  in  Fig.  90b,  which  represents  a ■ | 
physically  possible  steady-flow  solution.  However,  once  the  J 
major  conditions  are  satisfied,  details  of  the  real  flow  5 
structure  may  still  preclude  ever  achieving  a completely  steady  | 
state.  However,  levels  of  unsteadiness  may  still  be  considerably  5 
reduced,  as  the  experimental  results  indicate.  fj 
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FIG.  90a.  FLOW  OVER  A CONVENTIONAL  CAVITY  (this  configuration 
cannot  be  a steady  flow). 
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Beyond  the  fact,  that  a steady  flow  must  be  possible,  the 
flow  must,  also  be  stable  to  disturbances.  Small  vertical 
displacements  of  the  shear  layer  do  not  strongly  affect  the 
impingement  process,  since  the  slanted  edge  has  constant  slope. 
This  insensitivity  should  also  considerably  diminish  the  mass 
addition  and  removal  processes  that  occur  when  the  shear  layer 
Is  deflected.  The  constant  slope  also  allows  the  shear  layer 
to  find  its  own  natural  vertical  level  along  the  slanted  edge. 
This  level  depends  in  a complex  manner  on  the  pressure  and 
mass  flow  balances  that  must  be  maintained  in  the  cavity.  In 
contrast,  a rounded  edge  has  the  correct  slope  for  steady  im- 
pingement at  only  one  vertical  height,  which  may  not  be  the 
natural  position  for  the  shear  layer.  Furthermore,  the  geometric 
slope  change,  which  is  caused  by  a rounded  edge  when  the  shear 
layer  is  displaced  vertically,  is  such  that  it  adversely  affects 
the  stability. 

Actual  analytical  determination  of  the  optimum  slant  angle 
based  on  the  mean-flow  nroperties  and  characteristics  of  the 
shear  layer  is  a complicated  problem  in  fluid  mechanics.  A 
i jugh  order  of  magnitude  calculation  was  made  which  shows  that 
the  tangent  of  the  slant  angle  is  of  order  unity.  This  means 
that  typically,  the  slant  angle  will  be  a number  like  *45°  (as 
opposed  to  numbers  like  10°  or  80°).  In  fact,  the  current 
experiments  suggest  that  this  number  may  be  near  ^5° , at  least 
for  subsonic  speeds.  The  length  of  the  slanted  portion  must 
be  large  enough  to  accept  the  impingement  of  most  of  the  shear 
layer  thickness  and  to  allow  the  shear  layer  to  find  its  natural 
vertical  level.  Unfortunately,  this  factor  may  mean  a fairly 
large  slanted  region  at  the  trailing-edge  bulkhead,  particularly 
when  the  shear  layer  is  turbulent.  However,  reductions  in  size 
may  be  possible  by  carefully  testing  the  specific  configurations. 


APPENDIX  D 


FLOW  STABILIZATION  THROUGH  DETACHED  TRAILING-EDGE  COWL 

" ! ‘ ‘ ' ‘ ' 1 

* 

A detached  cowl  placed  a small  distance  upstream  cf  the 
tra iling-edge  bulkhead  of  a rectangular  cavity  has  a stabilizing 
effect  on  the  otherwise  oscillating  shear  layer.  To  be  effec- 
tive, this  cowl  must  have  a curved  surface  that  faces  the  flow. 

The  water  table  flow  visualization  studies  reveal  that  a half- 
cylindrical  cowl  stabilizes  the  flow  under  simulated  supersonic 
flow  conditions  (Fig.  91)*  For  the  tests  in  air,  an  airfcll 
shape  was  used,  as  shown  in  Fig.  76. 

We  believe  that  stabilization  is  achieved  by  accelerating 
the  shear  layer  over  the  curved  surface;  thus,  an  underpressure 
at  the  slot  formed  between  the  cowl  trailing  edge  and  the  cavity 
bulkhead  is  created.  This  underpressure  tends  to  suck  fluid 
mass  out  of  the  cavity  each  time  the  shear  layer  attempt:  to  in- 
ject mass  into  the  cavity;  thus,  a cancellation  effect  is  created. 
In  this  manner,  the  effective  strength  of  the  trailing-edge  mass 
addition  and  removal  process  is  greatly  reduced. 


FIG.  91.  WATER  TABLE 
THROUGH  DET 
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ular cavities,  exposed  to  high-speed  external  flow,  and  (2)  devise  and  eval- 
uate devices  that  would  ei'ther  substantially  reduce  the  amplitude  of  such 
pressure  fluctuations,  or  totally  suppress  the  occurrence  of  pressure 
fluctuations. 
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During  this  investigation,  the  understanding  of  the  complex  Interaction 
of  the  external  shear  layer  and  the  cavity  internal  fluid  medium,  which  con- 
stitutes the  generating  mechanism  of  high-intensity  periodic  pressure  fluc- 
tuations was  substantially  furthered.  Both  the  analysis  and  investigation  of 
the  physical  mechanisms  were  aided  through  extensive  shallow-water  flow- 
simulation  experiments. 

Large-scale  experiments,  using  the  NASA-Lewis  Research  Center  8 X 6 ft 
Supersonic  Wind  Tunnel  were  conducted  to  substantiate  some  of  the  analytical 
predictions  and  to  study  in  detail  the  aeroacoustic  behavior  of  cavities  in 
the  length-to-depth  ratio  range,  L/D,  of  2.3  to  5.5,  and  the  Mach  number  range 
of  0.8  to  2.0.  Detailed  information  was  obtained  on  the  normalized  levels  of 
the  first  three  resonant  modes  in  the  cavity  for  a range  of  cavity  length-to- 
depth  ratios  and  freestream  Mach  numbers. 

Several  concepts  for  pressure  oscillation  suppression  were  developed  and 
evaluated  in  wind  tunnel  experiments.  The  most  promising  concept  utilizes  a 
slanted  trailing  edge,  which  stabilizes  the  free  shear  flow  above  the  cavity, 
thus  effectively  suppressing  discrete-tone  generation. 
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4.0,  and  either  5.1  or  5.5*  Not  all  of  these  data  were  reduced; 
however,  all  data,  which  were  obtained  through  Sensor  3,  were 
reduced  and  are  available. 

Signals  from  Sensor  3 were  found  to  be  close  to  those  from 
Sensor  2 (on  the  leading-edge  bulkhead);  thus.  Sensor  3 signals 
can  be  considered  representative  for  the  entire  cavity .Internal 
leading-edge' region.  Figure  48  shows  1/3-octave  band  pressure 
level  spectra,  obtained  through  Sensor  3 for  seven  freestroam 
Mach  numbers  and  three  L/D  ratios.  Levels  are  given  in  dB  re 
20yN/nr  . Each  spectrum  is  composed  of  several  discrete 
tones  and  a broadband  noise  floor.  The  discrete  tones  correspond 
to  the  various  mode  frequencies;  in  general,  few  modes  with 
relatively  high  Intensity  appear  in  the  deep  cavity,  and  more 
modes  with  more  evenly  distributed  intensity  occur  in  the  shallower 
cavities.  In  a 1/3-octave  band  representation,  a discrete  tone 
determines  the  band  level.  Hence,  direct  comoarison  of  levels 
in  each  spectrum  for  corresponding  modes  is  possible. 

Following  the  L/D  =2.3  column,  in  'Mg.  48  we  can  pursue 
the  level  and  frequency  change  of  the  discrete  peak  at  250-  to 
315-Hz  bands.  This  peak  corresponds  to  a mode-2  oscillation  (see 
Sec.  6.7).  Levels  reach  a relative  maximum  at  transonic  "inch 
numbers.  With  increasing  Mach  numbers,  the  mode-1  oscillation 
in  the  160-Hz  band  becomes  significant.  Clearly,  most  of  the 
energy  for  such  relatively  dee;'  cavities  is  in  the  discrete 
modes,.  For  the  shallower  cavities,  discrete  tones  are  less 
pronounced,  and  relatively  more  energy  is  :ln  the  broadband  portion. 
This  figure  is  intended  to  give  an  overview  of  the  spectral  changes 
with  the  two  most  Important  parameter::-,  i.e.,  Mac!)  number  and 
lengtn-to-depth  ratio.  The  levels  would  he  those  encountered  by 
a cavity  in  an  aircraft  flying  at  the  Mach  number-dependent 
altitudes  as  shown  in  Fig.  28.  Each  altitude  corresponds  to  a 
certain  frecstream  dynamic  pressure. 

Since  levels  in  the  leading-edge  bulkhead  region  are  known 
to  scale  with  freestroam  dynamic  pressure,  (Smith  et  al . , 

1974;  Heller  et  al. , 1970),  the  ordinate  scale  of  the  Fig.  48 
spectra  can  be  converted  intq  20  log  (prr,K/q<0)  by  subtracting 
the  numbers  given  in  Table  1 from  the  flu'ctuatlng-pressure 
levels . 


